NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



THE MATERIALS PROCESSING RESEARCH BASE 

of the 

MATERIALS PROCESSING CENTER 


submitted by 

MERTON C. FLEMINGS, DIRECTOR 
MATERIALS PROCESSING CENTER 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
CAMBRIDGE, MASSACHUSETTS 02139 


submitted to 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
400 MARYLAND AVENUE S.W. 

WASHINGTON, D.C. 20548 


ANNUAL REPORT FOR FY 1981 
PREPARED UNDER GRANT NO. NSG 7645 


(NASA-CK-169096) TAB MATERIALS PROCESSING N82-27J82 

RESEARCH BASE OF lilE MATERIALS PROCESSING THRU 

CENTER Annual Report for EY iy81 N 82-4 74 Ob 

(Massachusetts Inst, of Tech.) 295 y Unelas 

HC A 13/M F A0 1 CSCL 110 03/2 3 2 8272 




Si ilfT;':*" 'I I3h\i. KJ i, 1 :* 
Si- », ■w|v. abv \.f V A# KvriJ&i lifcJ 


Massachusetts Institute of Technology 
Materials Processing Center 
Cambridge, Massachusetts 02139 
Tel, (617) 153-3233 

June 28, 1982 



Dir. Louis Testardi, Manager 
Materials Processing In Space Office 
NASA Headquarters, EN-1 
600 Independence Avenue 
Washington, D.C. 20546 

Dear Lou, 

Enclosed is the Fiscal Year 1981 Annual Report of the Materials Pro- 
cessing Center to the National Aeronautics and Space Administration for 
Grant No. NSG-7645 which is entitled "The Materials Processing Research 
Base of the Materials Processing Center," Three (3) copies are enclosed 
and two (2) copies have also been sent to the NASA Scientific and Techni- 
cal Information Facility. 

The Annual Report is presented in four (4) sections beginning with 
0 description of the Materials Processing Center and summary of structure, 
goals, and current activities. Section II reports on the initiation of 
the Materials Processing Center/Industry Collegium, which now has more 
than thirty (30) industry members. Section III reports on the specific 
research activities within the thrust areas of nucleatlon and rapid solidi- 
fication, fluid flow in crystallization processes, and adaptive materials 
processing which are supported by NASA Grant NSG-7645. Section IV outlines 
the overall materials processing activities at Massachusetts Institute of 
Technology, many of which benefit either directly or indirectly from the 
NASA supported activities, and vice versa. 

The Materials Processing Center continues to enjoy much success and 
has a current total research volume in excess of $3 million per year. 
Substantial progress has been made in broadening the scope of interaction 
with industry with the formation of the above mentioned Collegium. This 
Collegium not only provides a forum for direct interaction with industry, 
but will also provide fellowships, seed funds for research projects, etc. 

Our success to date is directly attributable to the founding and continuing 
support of NASA, for which we are very grateful. We look forward to con- 
tinued collaboration with NASA, 


Sincerely yours, 



Merton C. Flemings 

Professor Merton C, Flemings, Director 
Room S-407 
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I. INTRODUCTION AND SUMMARY 
THE MATER IALS PROCESSING CENTER 


This is the second annual report of the Materials Processing 
Center, describing activities during calendar year 1981. 

Our activities grew greatly during this second full year of 
the Center's existence. Total funding of the Center grew to 
an annual rate of $3.5 Million, and we are particularly 
pleased that of this total, approximately $800,000 is 
funding by industry. 

A generous gift from Toyota Motor Company to MIT established 
the endowed Toyota Professorship of Materials Processing, 
which has been awarded to the Center Director, Professor M. 
C. Flemings. The Associate Director of the Center, 
Professor H. K. Bowen, has been awarded the endowed Ford 
Professorship of Engineering. In addition, Norton Company 
has provided funds for three years to permit hiring an 
Assistant Professor of Materials Processing, and a search 
for this new faculty member is underway. 

A major step of the Center has been to institute the 
Materials Processing Center/Industry Collegium, which is 
described in detail in a separate section of this report. 
The central aim of this Collegium is to provide a mechanism 
for improved interaction of MPC faculty, staff, and students 
with industry. This Collegium is now underway and is 
operating in close cooperation with MIT's Industrial Liaison 
Program. 

The central thrust of activities of the Materials Processing 
Center remains the transition or "action" stages of the 
well-known "Materials Cycle" (see Fig. 1). These comprise 
extraction, processing into bulk materials, processing into 
engineering materials, fabrication, recycling, and disposal. 
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Figure l» Materials Processing steps, as transition stages In the 
Materials Cycle. (Taken from the "Materials Cycle" of 
the Cosmot Report, Materials and Man's Noodt ) 


We continue to view Materials Processing as the engineering 
field that seeks to control structure, shape, and properties 
of materials, and to do so in a cost effective way with 
acceptable social costs. It lies at the heart of the 
broader field of Materials Science and Engineering, linking 
basic science to societal need and experience (see Fig. 2). 
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KNOWLEDGE KNOWLEDGE 


Figure Z< Materials Processing as part of Materials Science' and Engineering 
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CURRENT CENTER ACTIVITIES 

The Materials Processing Center is organized to enable it to 
conduct, or facilitate conduction of, three specific types 
of research programs: (1) discipline oriented research 

problems of the type ordinarily conducted by an individual 
faculty or staff member, and his associates and students, 

(2) technologically related problems that are so broad in 
scope that they require effective focused cooperation of a 
team of faculty or suaff members, and (3) development of 
processing "science" or a "materials processing base" that 
cuts across traditional disciplinary and materials boun- 
daries. 


The current annual operating budget for the Center is $3.5 
million dollars, of which approximately 0 comprises 

the central "nucleus" of a broad based grant from the 
National Aeronautics and Space Administration. Table I 
lists the programs of the Center. These cover a wide range 
of activities that fall within four central thrust areas of 
the Center — areas which the Center views as being at the 
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heart of its research charter* These are* 

• Process innovation and development of new materials 
through processing. 

• Mathematical and experimental modeling of processes 

• Computer aided processing, 

© Economic issues relating to materials processing. 

The research programs of the Center/ listed in Table I 
involve 19 faculty members/ 58 graduate students/ 21 and 
staff members. Industrial interaction is a keystone of the 
Center and a number of programs within the Center involve 
such interaction. Approximately $800/000 of the Center's 
funding is provided by industry. 

INDUSTRIAL ADVISORY BOARD AND THE CENTER INAUGURATION 

Our Industrial Advisory Board, comprised of industrial 
leaders in materials related industries, continues to be 
essential in its role of advising the Center. The members 
of this Board and their affiliates are listed in Table II. 
The Board held its second meeting on March 5, 1981 and sub- 
mitted a report of its conclusions and recommendations to 
the Dean of Engineering. An annual meeting is planned. 
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TAB LK 11 


* Advis ory B oard Cor 

Materials >roocs slinf Con ter , M. l.T, 


nr. Turner All toy , Jr.*' 

Ho s o a r o h Sclent i, s t 

Dow Chemical Central uo no arch 

Dr. Donald J . Biiekwedo 
Vico President & Director 
of Research 

Bethlehem Steel Corporation 

Dr. Kenneth J. Brondyke 
Director 

A! coa Labor a tor i e s 
Aluminum Company ol America 

Dr. Harris M. Bur to 

Chiei ol Metal k Ceramics niv. 

Wright Patterson Air to too Base 


Dr. John R, Car rut hers 
Laboratory Project Manager 
Solid State Laboratory 
Hewlett Packard Labs 

Lr. Paul A. Floury 

Director/ Materials Research Lab. 

Be 1 1 Lo born tor ies 

Dr. Rodney F., Hannoman 
Vice President, Corporate Research 
k Development 
Reynolds Metals Division 

Mr. Frank B. Herlihy 
Vice President, Group Executive 
- Hydraulics 
Ab e x Co v p o rati o n 


Mr. Winston R. Hindle, Jr. 
Vice President, Corporate 
Operations 

Digital Equipment Corporation 

Mr, John R. Hutchins, III 
Vice President & Director 
of Research & Development 
Corning Glass Works 

Dr* Anthony D, Kurtz 
President 

Kulite Semiconductor 
Products Inc. 

Dr. Horace N. Lander 
Senior Vice President, 
Research & Development 
AMAX Molybdenum Division 

nr. George Mayer 
Director, Metallurgy 
8, Materials 
Science Division 
Army Research Office 


Dr, Robert Mohr ab inn 
Chief, Metallurgy Division 
National Bureau of Standards 
U.S. Department of Commerce 

Mr. John H, Mori son 
Chairman of the Board 
Hi tchinor Manufacturing 
Company, Inc. 


* We regret the passing ol. a highly valued member of our Advisory 
Board, Dr, Turner Allrey, Jr, 
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Dr, Richard K. l’itlor 
Dr. Vice President-Technical 
Director 
He search Center 
Allegheny hud ium steel 
corporation 

Mr. Hichard P, polich 
Chairman 

Tall ix Xncorpora ted 

Mr, James F. Bochin 

vice president 

Turbine Components Division 

Equipment Group 

TRW Xncorpora ted 

Dr. Goo mo 8. Reichenbaeh 
V i c e P re s ident , Do n d e d 
Abrasives 
Norton Company 

Dr, Edward I. Snlkovitz 
Director , Research Programs 
Cilice of Naval Research 

Dr. Maurice B. Shank 
Director, Engineering Tech. 
Commercial Products Division 
Pratt & Whitney Aircraft 
Group 

Un i ted Terhnolog ies 
Corporal ion 


Mr. Charles U. Smith, Jr, 
Chairman of the Board 
S I PCO Ind us tries , 

Incorporated 

Dr, M, A. Steinberg 
Director of Technology 
Applications 

Lockheed Aircraft Corporation 


Dr. Louis Testardi 
Manager, Materials 

Processing in Space Of Lice 
NASA Headquarters 

Mr. Thomas R. Wiltse 

General Manager, Central Foundry 

General Motors Corporation 


Dr. Edward S, Wright 

Director 

AMMRC 



I I . THE INDUSTRIAL COLLEGIUM 


INTRODUCTION 


The Materials Processing Center/Industry Collegium was 
officially begun on November 1, 1981 after extensive 
discussions with MPC Advisory Board Members, faculty and 
staff of the Center, and the MIT Administration, The 
central aim of this Collegium is to foster mutually 
productive relationships between industry and Center 
faculty, staff, and students. Response of industry has been 
strongly positive and we have now ahead of us much to do to 
make the Collegium as effective as it can be, 

A full description of the Collegium and its plan of 
operation has been prepared in the form of a brochure given 
to those companies interested in joining. That brochure is 
reprinted here in its entirety and comprises the -remainder 
of this Chapter of the report. 
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An Invitation to join 

the Materials Processing Center/Industry Collegium 
at the Massachusetts Institute of Technology 
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The Materials Processing Center 'Industry Collegium is a partnership 
established to promote innovative research and development programs, 
mtormation exchange, and personnel exchange Its aim s to assist 
industry to protit trom Ml T leseaich in materials processing, and tor MIT 
to benefit through improved industrial interaction 

through special publications, seminars, workshops, and personnel 
exchanges, the Collegium provides a forum for faculty, staff, students, and 
industrial representatives to exchange ideas, develop innovative research 
topics, and carry out cooperative or sponsored research and development 
programs oriented towards specific goals 
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The Materials Processing Center was established in the School of 
Engineering at MIT in 1979 Tne Ce r ter aims to work in concert with 
industry and government on broad technological problems relating to 
materials production and forming including aspects of those problems 
which relate to economic or societa' costs and benefits Some current 
central activities of the materials P'ocessing center are the following 

• Electronic Matena's Procesvng Research 

• Welding and Joining 

• Powder Piocessing/Rap'd Sc idi* ca'ion 

• Materials Systems 

• Ceramic Processing and E"g neeung Research 

• Electroprocessing 

• Polymer Processing 

• Mrunematica! and Physica Modeling 

• Corrosion Research 

• Mechanical Prope' ties ana Met i Forming 

• Composites and Non-Dest r ..ct've Evaluation 

• Solidification Processing 

• Computer Aided Processes 

• Chem . ai Metallurgy 

The foregoing activities o' "e Center involve over 20 f acuity 
members. 50 research staff, and *-*0 graduate students. 

An importont par! of the KVe'als Processing Center is ris 
Collegium by wnich complies i ' >ours can stay abreast of new 
technology m the Materials Proc* ng a 'oa, influence the course of 
processing research at MIT meet >.cu ty evolved in piocessi'ig and 
follow tne professional developed of unoe*giaduate and graduate 
students Close contact between representatives and students is 

expected to be a particularly frui" asooct of Center activities 


v VYiiiTE PHOTOGRAPH 
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The Center plays a special role in acting as host and intermediary 
between MIT researchers rnd members of your company Short visits for 
research discussions are readily arranged and the Center also 
encourages and facilitates personnel exchanges of all types including 

• Student cooperative work/study programs 

• Faculty summer work periods in industry 

• Visiting scientist appointments at MIT by industrial colleagues 
on extended leave 

• Shorter visits for special purposes by industrial colleagues 

In addition, to keep members abreast of current research oppor- 
tunities, the Materials Processing Center makes available periodically 
abstracts of current materials processing research, and provides a 
comprehensive annual summary of materials processing research at MIT. 


Collegium members receive Opportunity Briefs not less than twice a year 
Each Opportunity Brief outlines an area of ongoing research in Materials 
Processing at MIT which has long or short range business potential 

Following the publication of an Opportunity Brief, interested Col- 
legium members attend a Workshop where they meet oirectly with faculty 
and students whose work has been summarized in the Brief With partici- 
pation fmm industry, university ana government agencies, Workshops 
allow for a fuller explo'ation of the research and its potential applications. 

In addition, Workshops give membeis an opportunity to suggest directions 
for future research and to lay the groundwork for future participation in 
research programs. Some examples of planned Workshops include 

• Materials Cost and Availability 1980-2000 

• Mathematical and Physical Modeling of Materials Processes 

• High Energy Source Processing of Materials 

• Innovations n Welding ana Joining 

• Bapid Solidification Proceso.ng 

• Strategies for Systems Management of Materials Resources 

t 




17 



We cordially invite your company to deccrne a member of the Materials 
Processing Center/Industry Collegium The Collegium membership tee is 
$10,000 per year Special arrangements are possible for companies 
which are also members of the MIT Industrial Liaison Program, or are 
otherwise supporting programs related to trie Materials Processing 
Center 

To | 0 iri the Collegium, or for more information, please write or call 

Matonals Processing Center 
Room 4 415 

Massachusetts Institute of Technology 
Cambridge Ma 02139 

Professor Merton C Flemings. Director. (617)253-3233 
Prolessor H Kent Bowen Associate Director, (617)253-6892 
Dr George B Kenney. Assistant Director (617)253-3244 
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I.TX, MATERIALS PROCESSING RESEARCH BASE 

Research activities conducted under the Materials Processing 
Research Base during this second year are in three areas: 

(A) nucloation and rapid solidification, (B) fluid flow in 
crystallization processes , and (C) adaptive materials pro- 
cessing. The individual programs are listed below, and 
detailed research summaries follow. 

A. NUCLEATIQN AND RAPID SOLIDIFICATIO N 

» l i mn hi i »» —■■■ m . n »i....'i ■■■ . mu m i— —«— i^wiiwi 


1. S, M. Allen, "Processing of Sendust-Type Soft Ferro- 
magnetic Alloys." 

2. M, C. Flemings and J. Szekely, "Convection in Grain 
Refining ." 

3. N. J, Grant, "The Structure and Properties of Rapidly 
Solidified High Alloy Aluminum Materials," 

4. F. J. McGarry, "Rapid Solidification of Polymers." 

5. G. J. Yurek, "Development of Rapidly Solidified Oxida- 
tion Resistant Alloys." 

6. J. S. Haggerty, "Laser Materials Processing Facility," 

B. FLUID FLOW IN CRYSTALLIZATION PROCESSES 

1. R. A. Brown, "Fluid Flow in Crystal Growth: Analysis of 

the Floating Zone Process." 

2. D. Roy lance, "Numerical Modeling and Optimization for 
Polymer Melt Processing Operations," 

3. D, R, Sadoway, "Studies of Materials Electroprocessing in 
Molten Salts." 

4. A. F, Witt, "Heat Flow Control and Segregation in Direc- 
tional Solidification." 

C. ADAPTIVE MATERIALS PROCESSING 

1. T. W. Eagur, "Adaptive Control of Welding Processes." 
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A. NUCLEATION AND RAPID SOLIDIFICATION 

PROJECT Al : PROCESSING OF SENDUST-TYPE SOFT" FERROMAGNETIC 

ALLOYS 


Principal Investigators Prof. S. M. Allen 
Personnel: Mr. E. P. Kvam 

RESEARCH ABSTRACT 


The specific processing techniques which will optimize the 
magnetic properties of an Fe-Si-Al-Ni Sendust-type alloy are 
under investigation. The alloy being studied has magnetic 
properties comparable to Permalloys but with far less nickel 
content and hence less cost. Our objective is to develop a 
thorough understanding of the ways in which solidification 
processing, deformation processing, and heat treatment can 
be best utilized to produce material in suitable shapes, 
with excellent magnetic characteristics. 

RESEARCH SUMMARY 

I. INTRODUCTION 

The Senduat composition, Fe-9.6 w/o Si-5.4 w/o M, has long 
been recognized as an alloy with exceptional characteristics 
as a soft ferromagnetic material. These are: High 
permeability, low coercivity, high electrical resistance, 
and reasonably good saturation magnetization. These 
excellent properties result primarily from the fact that the 
composition has near-zero magnetostriction as well as 
near-zero magnetic anisotropy. Unfortunately, the material 
is also very brittle and applications in the past have been 
limited mainly to cast shapes and to forms made from 
compressed powder. 
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Two recent studios give promising results with respect to 
processing-related improvements Cor these alloys. The 
first is a report by Tsuya et al. (I EBB Trans, on 

Magnetics, 15, 1149 (1979)1 of the use of a melt-spinning 
rapid-solidification technique to produce Scnduat alloy 
ribbons of fine grain size. As cast, these ribbons were 
ductile onouqh to be rolled to 25% of their original 
thickness. After annealing treatments to optimize magnetic 
properties, however, the material was embrittled. The 
second study, by Yamamoto and Utsushikawa [Trans. Japan 
Inst. Metals, 19, 3261 (1970)], focussed on alloy 
modification through 3 w/o additions of nickel to an Fe-6 
w/o Si-4 w/o A i alloy. The resulting material was called 
Hupei' Scnduat because many of its magnetic properties were 
superior to those of the ternary Senduat alloy. Although 
the mechanical properties of the alloy were not studied, 
there is reason to expect that it will have better ductility 
than Senduat. Nickel additions to soft ferromagnets are 
generally believed to contribute to both magenetic and 
mechanical softness. Also, Super Senduat has approximately 
one-third less aluminum and silicon, each of which 
embrittles iron. 

I I . RESEARCH IN PRO GRESS 

Senduat alloy buttons have been prepared by vacuum arc 
melting. These ingots are then remelted and cast onto the 
inside of a rapidly rotating drum to produce thin ribbons 
approximately 2 mm wide. The microstructure of the as-cast 
ribbons is being studied to determine grain size, grain 
shape distribution through the ribbon, and crystal 
structure. Ductility of the as-cast ribbon material is also 
being characterized. 

The magnetic behavior of the ribbon material is under 
investigation. Initially, we are studying the influence of 
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annealing treatments on coercivity. As cast, the ribbons 
have a coercive force of approximately 2 Oersteds* 

Annealing at 900°C, followed by slow cooling, reduces the 
coercive force to approximately 0.7 Oersteds. Further 
reductions in coercive force should be possible through 
thermal treatments, and these are being pursued. At these 
low coercivities , the surface condition of the ribbons may 
also be important. As cast, the ribbons are rough on the 
broad face away from the quenching wheel. Deformation of 
the as-cast ribbons by rolling will be tried as a method to 
smooth the surfaces. 

Phase transformations in Super Sendust are being studied via 

transmission electron microscopy analysis of samples 

annealed over a range of temperatures. An order-disorder 

transition has been located at about 675°C. The high 

temperature phase has the B2(CsCl) structure, while the low 

temperature structure is DO (Fe Si or Fe A1 ) which is 

3 3 3 

reported to be the desired structure for the best magnetic 
properties in the ternary Senduet alloy. The primary 
defects in these phases are antiphase boundaries (APBs). If 
antiphase domains are fine enough, APBs may interfere with 
magnetic domain wall motion. Our transmission electron 
microscopy results enable us to estimate antiphase domain 
coarsening kinetics at different temperatures, and this 
information should prove useful in optimizing annealing 
schedules for reducing coercive force. 

In summary, our aim is to use a rapid-solidification 
processing technique to produce microcrystalline Super 
Senduet alloy ribbons. Such material, as cast, has a useful 
shape and because of its small grain size and homogeneity, 
should have reasonably good ductility. Without additional 
processing, however, it has relatively poor magnetic 
properties. A combination of deformation processing and 
thermal annealing treatments is being designed to optimize 
magnetic properties. 
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PROJECT A2 : CONVECTION IN GRAIN REFINING 

Co-Principal Investigators: Prof. M. C. Flemings 

Prof. J. Szekely 


Personnel: Dr. 

Dr . 
Dr . 
Dr . 
Dr . 
Dr . 
Dr . 
Mr . 
Mr’. 


R. Abbaschian 

C. W, Chang 

M. Choudhary 

N. El-Kaddah 
J. McKelliget 
Y. Shiohara 
G. Oreper 

G. M. Chu 

D. Maclsaac 


RESEARCH SUMMARY 


The scientific aim of this program is to obtain a better 
understanding of the relationship between fluid flow pheno- 
mena, nucleation, and grain refinement in solidifying metals 
both in the presence and in the absence of a gravitational 
field. One ultimate technical aim is to determine ways to 
achieve significant grain size reductions in hard-to-process 
melts; another is to understand the effects of undercooling 
on structure in solidification processes, including rapid 
solidification processing; a third is to better understand 
how to control this undercooling to improve structures of 
solidified melts. The project has been divided .into two sub- 
tasks. The actual grain refining and supercooling experi- 
ments are being carried out by Professor M. C. Flemings' 
group, together with thermal modeling of the solidification 
process, while the main thrust of the work of Professor J. 
Szekely' s group has been to study the heat and fluid flow 
phenomena in the levitated meta] droplets. 
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The main accomplishments in the heat flow, fluid flow, 

modelling area are; 

1. A computational capability has been developed to determine 
the electromagnetic force field, the fluid flow field and 
the temperature field in induction stirred systems, 
including contained cylindrical melts and levitated 
spherical melts. 

2. Calculations were carried out for a variety of conditions, 
including heat and fluid flow in a metal held in an 
inductively stirred cylindrical crucible and levitation 
melted specimens both on the ground and in a zero gravity 
environment. 

3. Calculations have shown that the fluid flow field is 
markedly different for ground base and for zero gravity 
conditions. The theoretical predictions for the velocity 
fields and the temperature fields were found to be in 
good agreement with experimental measurements reported in 
the literature. Experimental measurements are being 
conducted to study the velocity fields in a low melting 
alloy system using hot wire anemometry. This is a 
pioneering undertaking which shows considerable promise 
in testing the model predictions. 

4. The techniques developed for solving MHD type problems in 
molten metal and glass systems and the results generated 
are thought to have made an important contribution to 
this overall field. 

The more important and significant results of the undercooling 

and structures research are; 

1. Nickel base alloy samples of approximately 1 g have been 
successfully levitated in inert atmospheres and under- 
cooled by amounts up to 270°C, and a wide range of grain 
sizes and solidification structures obtained, depending 
on amount of undercooling and cooling rate. 

2. Two important innovative techniques have been developed 
to obtain much larger amounts of undercooling in high 
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temperature (iron# nickel, and cobalt base) alloys. In 
one of these, the metal is melted and then "emulsified" 
(stirred into fine droplets) in a molten oxide or salt. 

In the seond, small pre-alloyed metal droplets are 
interspersed at room temperature with finely crushed 
oxide or salt. The admixture is then melted. 

3. Extremely large undercoolings have been obtained in the 
above two methods because of the fine particle size and 
cleansing action of the slag. For example, in 316 
stainless steel, undercoolings up to 480°C have been 
obtained. This is approximately 300°C greater than any 
undercooling that has been previously recorded for any 
iron base alloy, and we do not think we have yet reached 
the maximum number. At the higher undercoolings, 
metstable BCC steel is obtained in place of the equili- 
brium FCC structure. These results are of fundamental 
importance in understanding structures produced by 
"Rapid Solidification Processing." They also point the 
way directly to new engineering processes for producing 
controlled solidification structures. The results have 
special implications for space processing as discussed 
below. 

4. Emphasis of the experimental work at the present time is 
on increasing the amounts of undercooling obtainable, 
and therefore the types of structures obtainable through 
(a) use of alternate emulsification media, (b) increas- 
ing rate of heat extraction, and (c) process variations. 
One factor that limits applicability of the process in 
earth processing is the need, at 1 g, to have the slag 
viscosity high enough to prevent settling and coalescence 
of the metal particles. 
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1 . THE MAT H EMATICAL AND PHYSICAL MODELLING OF FLUID FLOW IN 
CONTAINED AND IN CONTAINERLESS MELTS 


In containerless processing applications the metallic speci- 
mens are positioned with the aid of an electromagnetic force 
field. In ground based applications of this positioning 
technique (viz levitation melting) the force field has to be 
very strong, while in space processing applications a much 
weaker force field should be sufficient. 

In both these cases the applied force field will generate 
an electromagnetically driven flow, the quantitative 
assessment of which is thought to be crucial for both the 
interpretation of measurements and for the rational design 
of in-flight experiments. 

This work has the following components: 

(i) Calculation of the electromagnetic force field, 

(ii) Computation of the fluid flow fields resulting from 
these electromagnetic force fields, 

(iii) The verification of the calculations through experi- 
mental measurements, 

(iv) Interfacing of these results with the grain refining 
studies . 


( i ) Calculation of the Electromagnetic Force Field 

The calculation of the electromagnetic force field was per- 
formed by using the concept of the mutual inductances. This 
work was undertaken by the Space Science Center of the 
General Electric Corporation, who delivered a flexible com- 
puter program to us ^ . This program or some modification 
thereof was being utilized to compute the electromagnetic 
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force field for spherical specimens processed in a con- 
tainerless mode and for cylindrical melts held in crucibles, 

( i i ) Computation of the Fluid Flow Field 


A formulation has been developed to represent turbulent, 
recirculating flow as driven by the combination of buoyancy 
and electromagnetic forces in levitated specimens and also 
in molten metals held in cylindrical containers. These 
latter calculations were performed because extensive 
measurements have become available concerning such systems, 
against which the predictions may be tested. 

The general vectorial form of the governing equations is 
given as; 


V’U " o 


( 1 ) 


equation of continuity 


(t) 


- Vpyii - Zp - ft,' + P P 


( 2 ) 


where V’t is the turbulent stress tensor and F is the body 

b 

force field, which in general will include both the buoyancy 
and the electromagnetic forces. 


The thermal energy balance equation takes the following form; 

UVT « VkVT (3) 

where k is the effective thermal conductivity, 
e 

The system of equations (1-3) was solved numerically, uti- 
lizing the k - e model to represent the components of the 
turbulent stress tensor. 
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In essence three sets of calculations have been performed; 

(a) The heat and fluid flow field was modelled in a levi- 
tation melted specimen under ground based conditions. 

(b) Calculations were carried out modelling heat flow and 
fluid flow in an electromagnetically positioned metal 
sphere under zero gravity conditions, and finally, 

(c) Calculations were carried out to represent the beha- 
vior of a molten mercury pool which was agitated by a 
current passed through induction coils surrounding 
the melt. 

Some calculated results are given in the following and some 
additional computed results will be given subsequently in 
conjunction with measured data. 

Figs. 1 and 2 sho^ the computed lift force, acting on a 1 g 
levitated metal droplet, under ground based conditions, for 
coil currents of 200A and 300A respectively^ 2 '^ . Also shown by 
the horizontal line is the weight of the levitated specimen. 

It is seen that the lift force would not be adequate to 
levitate for a coil current of 200A, but would be sufficient 
for a coil current of 300A. This finding agrees with ex- 
perimental measurements. Fig. 3 shows the computed velocity 
field and the computed temperature field for a 1 g levitated 
specimen, under ground based conditions. It is seen that 
the linear fluid velocity is quite high, of the order of 0.2 
- 0.3 m/s and that the temperature in the sphere is not 
completely uniform ' '• These predicted linear melt velocities 
were again found to be in good qualitative agreement with 
observations. Figs. 4 and 5 show the computed velocity and 
the temperature fields and the computed turbulence charac- 
teristics for a levitation melted Beryllium sphere, under 
zero gravity conditions. The input parameters for the 
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modelling calculations correspond to those used experimen- 
tally in the SPAR 1 runs. It is noted that the fluid flow 
field differs quite appreciably from that computed for the 
ground based conditions, exhibiting four rather than two 
circulating loops. The theoretically predicted temperature 
appears to be in good agreement with that measured in the 
experiments , 


( ill ) Verification of the Predictions, by a Comparison with 
Experiments 


In this regard, three sets of experimental measurements have 
to be considered: 

(a) Mass transfer in a levitation melted sphere, 

(b) Experimental measurements in an induction stirred 
vessel, and 

(c) Experimental measurements in a molten metal pool agi- 
tated by the direct passage of a current between two 
electrodes . 

The former two refer to experiments conducted elsewhere, 
while the latter group involves experimentation in our 
laboratory. 


Mass Transfer in a Levitation Melted Sphere 


Experimental measurements were conducted on the rate at 
which a levitation melted iron sphere picks up carbon, when 
contacted with a C0 2 /C0 mixture. The overall rate of this 
process depends on both the rate of mass transfer in the gas 
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phase and on the fluid flow and circulation in the metal 
sphere. Fig. 6 shows a comparison between the experimental 
measurements and the theoretical predictions!. It is seen 
that the agreement is quite good, providing direct confir- 
mation of the appropriateness of the mathematical model. 


Fluid Flow in an Induction Stirred Vessel 


These experiments were conducted in the Engineering 
Department at Cambridge University, by Dr. Hunt and Mr. D. 
Moore, who kindly provided us with their data. 

Fig. 7 shows a sketch of the experimental arrangement, which 
indicates a vessel holding molten mercury, which is being 
agitated by induction coils. Fig. 8 shows the computed 
radial body force field, while the corresponding axial com- 
ponent is shown in Fig. 9. The experimentally measured velo- 
city field is shown in Fig. 10, while the corresponding 
theoretically predicted values are given in Fig, 11. The 
very good agreement is readily apparent. Fig. 12 shows a 
more detailed comparison of the theoretically predicted 
(broken lines) and the experimentally measured velocity 
field, again giving very reasonable agreement considering 
the fact that the theoretical predictions were based on 
first principles. Finally, Fig, 13 shows a comparison bet- 
ween the experimentally measured and the theoretically pre- 
dicted values of the turbulence intensity, again giving very 
good agreement. 


Experimental Verification of the Model in this Laboratory, 


In order to provide a direct experimental verification of 
the predicted electromagnetically driven flow fields 
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obtained in this laboratory under controlled conditions, an 
apparatus has been constructed in which molten Woods Metal 
is being agitated by the passago of a current between two 
submerged electrodes, as sketched in Pig* 14. The resultant 
velocity field will be measured by using a hot film anemo- 
meter, the output of which will be analyzed by a microcom- 
puter. The apparatus has been constructed and the 
preliminary results are being obtained at present. 


( iv ) Interfacing the Modelling Work with the Grain 
Refining Experiments . 


The results of the mathematical modelling efforts do provide 
a sound basis for helping with the interpretation of the 
grain refining experiments, conducted by Professor Flemings’ 
group. Work has been initiated in this area, particularly 
by modelling the solidification of supercooled melts. 

The work done to date may thus be summarized by stating that 
major advances have been made regarding the modelling of 
turbulent electroinagnetically driven recirculating flows in 
both spherical and cylindrical geometries. A capability has 
been developed to predict the velocity fields in electro- 
magnetically agitated melts, in both the spherical and the 
cylindrical geometry. 

The model has been verified by a direct comparison with 
experimental measurements conducted in other laboratories, 
for ground based conditions and work is in progress to per- 
form additional experiments in this laboratory with this 
objective. Predictions have also been made regarding the 
heat and fluid flow in a zero gravity environment, 
corresponding to one of the previous space experiments. The 
predicted temperature fields were in good agreement with the 


Lifting Force F^xIG^CN) 
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Figure Is The computed lift force acting on a 1 g iron 
droplet for a coil current of 200 A, 
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Figure 2 % The computed lift force acting on a 1 g 
droplet for a coil current of 300 A. 
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Figure 4: The computed velocity field and temperature 

field for a Beryllium sphere which is induction 
heated under zero gravity conditions. 
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Fiquro bs Comparison of the experimentally measured and 
the theoretically predicted carburization 
curves for a levitated iron droplet. 
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Figure 7: Sketch of the experimental arrangements for measuring 

the velocity fields in an induction stirred mercury 
pool under ground based conditions. 
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Figure 8: Comparison of the experimentally measured (discrete 

data points) and the theoretically predicted radial 
body force field for the mercury system. 
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Figure 10: 


The experimentally measured velocity field in the 
ground based mercury system. 
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Figure 12: Comparison of the experimentally measured and 

the theoretically predicted axial velocity 
component at one practical vertical position 
within the mercury system. 
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Figure 13: Comparison of the experimentally 

measured and the theoretically pre- 
dicted turbulence intensity in the 
mercury system. 
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Figure 14: 


Schematic sketch of the apparatus designed 
for measuring the electromagnetically driven 
velocity field in the Woods Metal system. 
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measurements. Work is currently in progress in order to 
relate these findings to the grain refining experiments and 
to the planning of an experimental program to be conducted 
in a zero gravity environment. 
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2 . EXPERIMENTAL PROGRAM ON NUCLEATION AND STRUCTURE IN 
UNDERCCOLED MELTS 


INTRODUCTION 


The central focus of this work has been on undercooling and 
structural refinements in droplets of molten metal levitated 
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in an induction field and/or by dispersion in a fluid 

carrier. Initial work, summarized in the last annual report 

( 1 ) 

and in a technical paper was on metal base alloys. More 

recent work has been on nickel base and lower melting point 

(2,3) 

alloys levitated in molten carr.ier fluids , Ma^or 
thrusts of current work are: (1) dispersion of molten alloy 

droplets in a high temperature fluid following the proce- 
dures developed by Perepezko and co-workers for lower 
melting point alloys, (2) obtaining a similar dispersion by 
room temperature mechanical mixing of particles of the metal 
and solidified liquid carrier, and (3) solidification of 
single relatively large droplets in a transparent fluid 
carrier, enabling high-speed temperature measurement of the 
recalescence and subsequent cooling behavior. 


During the last year, it has become evident that work on 
undercooling of metal has broad practical implications, both 
for space and terrestrial processing, and particularly in 
the important technological area of rapid solidification 
processing . 


As the solidification rate of a metal alloy is increased 
progressively from those rates encountered in conventional 
casting processees, the structure (dendrite arm spacing) is 
first observed to become finer, then supersaturation of the 
primary phase may be observed, metastable phases may be 
encountered, and ultimately the alloy may solidify as a 
metallic "glass." In all of these instances, substantial 
undercooling must exist in the liquid in front of the 
growing solid front, and in many instances, much of the 
undercooling may be developed in the liquid before any soli- 
dification begins. For example, the relatively homogeneous 
structures, sometimes termed "microcystalline , " found in 
atomized ultra-fine liquid droplets probably result from 
undercooling in the droplet before. The metastable BCC 

phase in atomized Fe-Ni alloys, and stainless steels, 
certainly results from undercooling before nucleation. 


(7,8) 
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Thus, the question of melt undercooling before nucleation 
and its influence on solidification behavior and solidifica- 
tion structure is fundamentally important with direct prac- 
tical implications. And, as has been pointed out by 
Perepezko , this undercooling can be studied not only in 
usual "rapid solidification processes" but also, and more 
directly, in any of several processes designed to produce 
significant undercoolings at relatively low cooling rates 
or, more exactly, at relatively low rates of heat extraction 
to the surroundings. Examples of such processes are under- 
cooling of bulk samples and undercooling of dispersed 
droplets. Of course, even though the rate of heat extrac- 
tion to the surroundings in these processes may be very 
slow, solidification rate immediately after nucleation is 
not. In the Case of small dispersed droplets, the cooling 
rate after recalescence may also be quite rapid because of 
the temperature difference between the recalesced droplets 
and the surrounding fluid. 


3. EXPERIMENTAL AND ASSOCIATED ANALYTICAL WORK. 


Experimental work has been on a low melting alloy, Sn-25% 
Pb, and on high temperature alloys, Ni-Sn and IN-100. 


Tin- 2 5% Lead Alloy 


The apparatus used for the low melting point alloy is shown 

in Fig. 1; it is modelled after the apparatus of Perepezko * 

’ (9,10-13) . . 

and co-workers . Droplet dispersions are produced m 

the apparatus by emulsifying a mixture of the liquid alloy 

and an organic carrier fluid (Polyphenylether ) , at a high 

speed (3800 rpm). The atmosphere in the furnace is argon. 

The temperature is kept at 25°C above the melting point 

during shearing, which continues for 40 minutes, producing 
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fine droplets dispersed in the oil. The emulsified samples 
are then transferred to sealed aluminum pans for analysis in 
a differential thermal analyzer (DTA) or a differential 
scanning calorimeter (DSC) . The emulsified samples are 
later ultrason.ically cleaned for observation of cross- 
sectional microstructures by scanning electron microscopy 
(SEM) . 

Our current work comprises primarily study of metal struc- 
tures as a function of solidification variables (under- 
cooling, cooling rate). Working with droplets in the size 
range 2 to 30 microns, we can readily obtain undercoolings 
up to about 89 °C and find the measured undercooling is a 
function of the cooling rate in the differential scanning 
calorimeter, Pig. 2. Pig. 3 is a typical example of struc- 
tures obtained. The dark areas are tin-rich and light areas 
lead-rich . 


Sn-25% Pb is a two-phase alloy at equilibrium - it would 
have approximately 71 wt% tin-rich phase and 29 wt% lead- 
rich phase after equilibrium solidification. Note in Fig. 3 
that these phases are distributed differently in two 
distinctly different two-phase regions. There is a small 
circular Pb-rich region (hereafter termed V^) that is 
clearly lead rich because of the large volume fraction of 
the light colored phase. There is also a larger Sn-rich 
region, V^, which has predominantly the darker tin rich 
phase. The volume fraction of decreases with increasing 
cooling rate before nucleation (and therefore with 
increasing undercooling, Figs. 4 - 6). Structure fineness 
also increases with increasing undercooling; i.e., spacing 
of B particles in the a matrix decreases, Fig. 4. Data in 
these figures are for droplets in the size range of 10-15 ym. 


In this work we have not found, nor have we been searching 
for, metastable phases such as those found by Perepezko and 
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co-workers. Our central aim has been to study the effects 
of undercooling and process variables on structure fineness, 
solute redistribution and, where present, supersaturation. 
Towards this end, we have begun preliminary study of the 
heat flow during solidification of these droplets, and of 
several models of solidification behavior, as discussed 
below. 


Heat Flow in Solidification of Undercooled Droplets 

There have been three studies made to date of heat flow 
, . . ... . (4,14,15) 

behavior m solidification of undercooled droplets. 

These studies have considered droplets that are undercooled 
and solidify in a gaseous medium (i.e., atomized droplets). 
All assume plane front solidification, and consider resist- 
ance to heat flow within the solid as well as at the 
droplet-gas interface. However, the Biot numbers achievable 
in practice in usual droplet atomization (Bi<<0.1) mean that 
heat flow can be considered Newtonian for purposes of calcu- 
lating overall solidification times. Even in Newtonian 
cooling, significant thermal gradiants can occur on the 
solid droplet and so calculations of G, thermal gradient; 
G/R, thermal gradiant divided by growth rate, etc., require 
consideration of heat flow in the solid. 


AH 


Ps v df s - IPs v f s c p + l^s’PL vc pl dT + qMt (1) 


In our preliminary work on heat flow analysis in droplets 
emulsified in a liquid^ we are using the basic methodology 
of Levi and Mehrabian , modified for the different 
surrounding medium and are beginning to explore application 
of these ideas to alloy (and dendritic) solidif icatition. 
Our basic model is one in which droplets are emulsified in 
oil and cooled uniformly to their nucleat.ion temperature, 
T^j. Upon nucleation, recalescence occurs creating a tem- 
perature difference between the droplet exterior and its 
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surroundings so that heat flow occurs both outwards to the 
surroundings and inwards to the liquid core. We assume the 
liquid droplets are sufficiently far apart that their ther- 
mal fields do not overlap. Our work to date has been using 
a Newtonian model so the heat balance relating heat flow 
from the droplet to fraction solidified is: 

AHp Vdf = [p Vf CS+ (1-f )p r VC L ]dT + oAdt , 
ss s s p sLp 

where: AH is the latent heat of fusion 

V is the volume of metal droplet 
f s is fraction solid 

p s , and p r are densities of solid and liquid 
Cp and Cp are the specific heats of solid and liquid 
q is the heat flux to the oil from a metal droplet 
A is the surface area of a metal droplet 
T is temperature and 
t is time 

Heat flow to the surrounding still liquid medium is given 
not by the convective heat flow correlation used when soli- 
dification occurs in a flowing gas but by the unsteady state 
form of the second law of heat conduction: 

3T 1_ 3__ 2 3T 

8 r “ a 2 3r (r 3R ) (2) 

r 

where a is thermal diffusivity of the surrounding medium and 

r is radius. In the example to be given, the solid is 

assumed to grow so the remaining liquid is spherical. In 

addition to the usual boundary conditions, the only other 

needed expression is a relation of growth velocity to melt 

undercooling. In addition to exploring kinetic coefficients 

(16) 

of the types proposed by Turnbull , and Cahn and co- 
workers , ( ) and employed by Levi and Mehrabian^^, we 
are examing the behavior of the growth front assuming the 
solid is growing dendr itically with a narrow-liquid-solid 
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range, and a dendrite tip velocity versus undercooling of 
the form 

R = aAT n 

where a and n are constants with n « 2. Relations of this 
form have been obtained by many investigators, and specifi- 
cally by Schaeffer and Glicksman . ) 

Combined Solidification and Heat Plow Model 


To relate the preceding thermal equations to solidification 
variables of an alloy requires some assumptions as to soli- 
dification behavior, and we are beginning to examine some 
different possibilities. One example of these is illu- 
strated in Fig. 7. Here, solidification is assumed to be 
partitionless between and T Q so solid composition is C Q . 
The "T 0 " temperature is calculated from the equal molar free 
energies of the liquid and the solid using the regular solu- 
tion approximation. The dependence of the temperature on 
the molar free energy is calculated from the Gibbs-Helmholtz 
relationship. Above T Q , the solid composition at the inter- 
face is the maximum thermodynamically possible solubility. 
This composition is calculated from the intersection of the 
molar free energy curve for the solid and the tangent line 
to the molar free energy curve for the liquid at the com- 
position as described by Cahn.^ 1 ^) 

Combining mathematical description of this solidification 
model with the preceding heat flow analysis enables pre- 
dicting details of heat flow behavior and some sample com- 
puter calculations for Sn-16 at% Pb (Sn-25 wt% Pb) are given 
in Fig. 8. The heat extraction during recalescence signi- 
ficantly affects thermal behavior and maximum recalescence. 
For example, a sample undercooled 60°C reaches a cooling 
rate of 105°C/sec just after its maximum recalescence tern- 
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pcrature, as it cools through the eutectic temperature of 
183°C (at this time, it is still 0.5 liquid). 

Solidification models such as this, to the extent they are 
proven valid, will ultimately permit prediction of effects 
of solidification variables on structure fineness, and will 
also permit prediction of amounts, compositions, and distri- 
butions of phases present, including any supersaturation. 

We are concentrating our efforts in this direction. Clearly, 
the problem is a complex one, and we intend to closely 
couple any computational work we do in this regard with 
experimental studies. 

High Temperature Alloys 

Our current experimental work on high temperature alloys is 
aimed at developing methods for undercooling droplets of 
iron and nickel base alloys. Three independent and related 
techniques are being employed. These are (1) emulsification 
of droplets in a high temperature liquid carrier (i.e., a 
technique similar to the technique of Verepestko and co- 
workers used for lower melting point alloys), (2) dispersion 
of fine droplets in a liquid carrier by intimate room tem- 
perature premixing of finely divided particles of the metal 
with particles of the solidified liquid carrier, and (3) 
undercooling of a single large metal droplet in a 
transparent slag to permit rapid temperature measurement of 
the recalescence and subsequent solidification. Exploratory 
work performed to date has been on Ni-Sn alloys and on the 
superalloy IN-100, but the techniques being developed are 
intended to be applicable to a broad range of iron and 
nickel base alloys. 

Droplet Dispersal by Emulsification 


The emulsification apparatus for high temperature alloys 
which has been constructed is shown schematically in Fig. 9. 
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Crucible and stirrer are made of aluminum oxide and all runs 
are conducted under argon or nitrogen atmosphere. Rotation 
speeds of up to 2000 rpm are used and droplet sizes in the 
range of 2-10 pm are readily obtained by separation of the 
smaller size droplets from the broader distribution produced 
by emulsification. Slags of several glasses including Pyrex 
and have been used with success, and experiments are 

also planned with salts. 

A high temperature differential thermal analysis ( DTA ) unit 
has been acquired and installed and is being utilized to 
study these materials. The unit is a Perkin-Elmer BTA-1700 
System which can be used at temperatures up to 1550°C, Fig, 
10a shows a typical DTA curve obtained for droplet samples 
of Ni-34% Sn in the 2-10 tarn range, which were emulsified in 
Pyrex. The curve shows undercooling of 219 °C, 

Fig. 10b shows an example of the structure obtained in the 
undercooled Ni-34% Sn dispersed in Pyrex. The electron 
micrograph is taken after DTA. Interphase spacing is the 
order of one micron, which is far finer than that of simi- 
larly undercooled bulk specimens. 

Range of droplet diameters obtained after DTA is typically 
2-100 ym. Current work in this part of the program is on 
obtaining more uniform and finer droplet dispersions through 
varying processing parameters including the liquid carrier, 
and on studying the undercoolings and structures obtained. 

Droplet Dispersal by Powder Processing 

Particles of IN-100, less than 20 ym in diameter , prepared by 
gas atomization have been mechanically dispered among fine 
particles of a solidified carrier fluid for analysis in the 
DTA. Heating of this powder mixture forms the metal droplet 
dispersion in a liquid carrier. Fig. 11 shows a DTA curve 
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for IN-1Q0 droplets which indicates undercooling of 300°G. 
Examination oi the droplets after DTA shows that some have 
combined to form larger droplets while others remained iso- 
lated # Examples of the microstructures in Fig. lib show 
that the larger droplets (20-* 100 tm) are two-phase while 
the smaller droplets <10 urn show no structure. This 
variation in microstructure is perhaps related to the two 
peaks on the DTA cooling curve which correspond to different 
degrees of undercooling. Similar single phase super- 
saturated micros true tur os have also been observed in under- 
cooled gas atomized droplets, but usually only in much 
smaller particles . ) 


barge Dr 0 £l e t Undercool tng 

One limitation of working with the dispersed droplets 
described in the previous sections, is that it is not 
possible to measure the thermal behavior of individual 
droplets. That is, it is not possible to pin-point the 
precise undercooling of any given droplet nor is it possible 
to measure its thermal history during and after recale- 
scence. Therefore, it is not experimentally possible to 
verify calculated curves such as those of Figure 8. 

We have therefore initiated a series of experiments on 
undercooling of individual droplets of high melting point 
alloys, in the size range of about 5 mm. The experimental 
apparatus is shown in Fig. 12. As shown, a small droplet 
is placed inside a Pyrex or quartz tube and surrounded by 
crushed glass. The sample is then induction melted and the 
surrounding glass is also melted. Vigorous stirring of the 
droplet by the magnetic field in the presence of the glass 
slag results, pre-sumably , in absorption of nucleating 
agents, leading to very large degrees of undercooling even 
in the relatively large droplet. Most importantly, thermal 
behavior of the droplet surface during recales-cence can be 
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directly observed. (This experiment can be similarly done 
by levitation melting without the surrounding glass.) 

The apparatus used to measure the rapid recalescence and 
subsequent cooling involves a sensing device, a data 
acquisition and storage device, and a data manipulation 
device. The sensing device is a Capintec Ratioscope III, 
two-color pyrometer sensing head with a 135 mm close-up lens 
kit to observe the small sample. This device provides two 
sig-nals from silicon photodiodes filtered at different 
wavelengths which can then be ratioed and calibrated at the 
observed metal temperature. 

The data acquisition and storage device is a Nicolet 
Explorer III Digital Oscilloscope which converts the analog 
output of the sensing head to digital bits of information 
which can more easily be converted into temperature values. 
The digital oscilloscope is interfaced to an HP-85 personal 
computer for data manipulation and interpretation. 

This measurement system records the thermal information 
directly from the photodiode sensors into the digital 
oscilloscope. In this manner, the thermal and time data are 
not delayed or dampened by amplification systems, and thus 
the fast recalescent profile can be properly recorded. 
Factors which could limit the ability to obtain accurate 
thermal data at the high speeds involved are: (1) the 

response time of the silicon photodiode sensors (which is 
the order of microseconds), and (2) the recording time for 
the digital oscilloscope, which is capable of 0.5 microse- 
cond per bit of information recording time. Initial experi- 
ments indicate that these response delays are inconsequen- 
tial compared to the recalescent rise time of an undercooled 
sample (which, for samples of this size, is on the order of 
milliseconds as will be seen below) . 
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In conducting thermal experiments, the metal sample is 
usually surrounded by a transparent viscous medium ouch as 
crushed Pyrex glass. The sample is placed in an induction 
molting coil within the viewing range of Gapintec R0S-1II 
two-color sensing head. The induction melter consists of a 
copper coil having two reversed turns which is coupled to a 
400 KC, 10 KW high frequency generator. The sample is stop- 
pered, evacuated and backfilled with argon to provide a 
clean environment around the sample. 

The sample is heated into the molten state and allowed to 
cool and solidify, After several thermal cycles, the sample 
comes into intimate contact with the glass and achieves 
better undercooling, A thermal cycle is then run, and the 
two-color measurements of nueleation and recalescence are 
digitized and stored within the Nicolet digital oscillo- 
scope. Conversion of the two-color measurements to 
temperature is being performed by a linear relationship of. 
the two-color ratio. In current experiments we are checking 
the validity and reproducibility of this linear relationship 
using thermal artifacts including the solidus temperature 
and the B-B 1 solid-solid transformation (in the Ni-Sn alloy) 
for calibration. 

One experimental thermal profile obtained to date is shown 
in Fig. 13. Recalescence time is approximately 30 millise- 
conds for these 5 mm diameter droplets. This curve is 
representative of early data acquired with newly obtained 
equipment and much experimental work is necessary to improve 
quality, prove reliability and usefulness, and determine the 
minimum size metal droplet on which the experiment can be 
successfully carried out. The aim of this research is to 
compare experimental curves such as Fig. 13 with calculated 
profiles such as those presented in Fig. 8, or with calcu- 
lated profiles based on other solidification models such as 
a pure metal model comparable to that of Levi and 
Mehrabian. ( 4 ) 
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Figure 1: Schematic diagram of the experimental setup for 

emulsifying droplets of molten low melting 

point metal in oil. • I 
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Figure 2: Effect of cooling rate in the differential 

scanning calorimeter on undercooling of 
Fn~25 wt% Pb droplets. (Data for 10~15pm 
diameter droplets.) 
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Figure 3: The cross sectional microstructure of a droplet 

of Sn-25 wt% Pb solidified with a DSC jooling 
rate of 10°C/min, and with 80°C of undercooling. 
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Figure 7: One possible solidification path for an alloy of 

composition C 0 . Solidification is partitionless 
from T n to T 0 . Between T 0 and Tr the solid that 
forms is of the maximum composition thermo- 
dynamically possible. Thereafter the solid 
forming is of the composition of the equilibrium 
solidus „ 
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Figure 8: Calculated temperature profile for a Sn-25 wt% Pb 

(Sn-16 at% Pb) alloy nucleated at three different 
temperatures, T n . 
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Figure 9: Schematic diagram of apparatus from emulsifying 

high melting point alloys. 
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(b) 

Figure 10: Differential thermal analysis (DTM and micro 

structure of Ni-34% Sn emulsified in Pyrex. 
DTA curve shows sample undercooled 219°C. 
Second peak is solid state transformation. 
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Figure 11: Differential thermal analysis ( DTA) and micro- 

structure of IN-100 alloy mechanically dispersed 
in Pyrex. DTA shows two peaks corresponding to 
undercoolings of 190°C and 300°C respectively. 
Photomicrographs ^at 1000X) show the laraer 
droplets have a J wo phase structure while smaller 
droplets arc apparently sinqlc phase. 
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PROJECT A3 : THE STRUCTURE AND PROPERTIES OF RAPIDLY 

SOLIDIFIED HIGH ALLOY ALUMINUM MATERIALS 


Principal Investigator: Prof. N. J. Grant 

Personnel: Mr, W. Wang 

RESEARCH ABSTRACT 

A series of 2024 type aluminum alloys modified by additions 
of 1 to 2% Li were studied to determine the role of the 
Cu:Li and the (Cu + Mg):Li ratios on resultant strength, 
ductility, notch-tens ile behavior, and crack propagation 
rates. U1 trasonically gas atomized powders, with quench 
rates of 10^ to 10^ °/s, atomized in an argon atmosphere, 
produced yields of powder such that almost 100% was finer 
than 250 microns. The powders are free of gases and 
porosity, are quite spherical, have few satellites (adhering 
fine powder particles) and are of uniform microstructure. 
Strength properties are such that yield strength is 20% 
greater than for lithium-free 2024 ingot alloy, tensile 
strength is 10% greater than that of 2024 ingot material, 
and ductilities are comparable. In terms of specific 
strength and specific modulus, these RS 2024-Li alloys are 
significantly better than IM 2024. 

RESEARCH SUMMARY 

An examination of the ternary diagram Al-Cu-Li shows an 
exceedingly complex system with numerous ternary 
intermetall ic compounds in addition to the expected 
CuMgAl 2 and A1 gLi (delta prime) phases. Small changes in 
the Cu:Li and (Cu + Mg):Li ratios can lead to the appearance 
of three ternary phases involving A1 , Cu, and Li. The 
presence of Li in the ternary phases decreases the amount of 
Li available for formation of AlgLi, one of the important 
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ageing phases in this system. The presence of 1.6% Mg, 0.5 
Mn, and up to 0.25% Cd further complicates the prediction of 
which phases will form and in what quantities. 

Accordingly, Cu was varied from 3.96 to 5.92%, Mg from 1.57 
to 1,92%, and Li from 1.25 to 2.17%. Both Mn and Cd were 
omitted from one of the alloys. See Table 1 for final alloy 
compositions . 

The high Cu alloy (5.92% Cu, 1.29% Li, l./sV* :»3g, 0.36% Mn, 
and 0.13% Cd ) had very poor ductility. Strength values were 
not importantly different than those for. the lower copper 
content, but ductility was decidedly smaller. The presence 
of primary CUAI 2 was readily evident; unfortunately this 
primary phase coarsened during solution heat treatment and 
probably accounted in part for the poor ductility values. 
Measured values of notch toughness and da/dN values were the 
poorest of the three alloys studied. 

To illustrate the room temperature tensile data, Table 2 
shows for Alloy 3 values of YS , UTS, elongation, and 
reduction of area for various heat treatments; in addition 
YS and UTS values are listed after correction for the 
decrease in density due to the presence of 0.96% Li. 

The decrease in density of Alloy 2 (1.63 wt pet Li) is 5% 
while the increase in Youngs modulus is 12.5% for a net gain 
of 18.3% in specific modulus. 

Further, Table 3 shows representative values of the Notch 
Yield Ratio, the Notch Tensile Ratio and Fracture Toughness 
(these values are more nearly equal to Kq rather than 
Kjjr. due to specimen geometry). 

Overall, these Li-modified 2024 alloys demonstrate the 
importance of the Cu:Li ratio on mechanical properties, with 
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Li having a major positive effect on both specific modulus 
and specific strength. The high fracture toughness values 
and the excellent notched yield and tensile properties 
indicate that alloys based on the 2024 compositions can be 
expected to possess excellent commercial potential. 
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PROJECT A4 : RAPID QUENCHING EFFECTS IN GLASSY POLYMERS 


Principal Investigator: Prof, P, J. McGarry 

Personnel: Dr. J. F, Mandell 

Mr. A. Agrawal 
Mr. H, D. Lee 

RESEARCH ABSTRACT 


Using a specially constructed microbalance for hydrostatic 
weighing, density changes in PVC thin film (with no 
additives, 30~100 ym thick), due to rapid quenching 
(~ 500 °C/sec ) through the glass transition temperature, have 
been observed. The more severe the quench, the greater is 
the free volume content. Isobaric volume recovery of PVC 
has also been studied by volume dilatometry. Both show 
aging of relaxing molecular rearrangements taking place an a 
linear function of logarithmic aging time at room 
temperature. Distribution of retardation times and Primak's 
distributed activation energy spectra have been applied to 
the volume recovery data. The concomitant changes in 
mechanical properties of PVC after quenching have been 
monitored by tensile creep and stress-strain to failure. 

All reflect the presence of excess free volume content, due 
to rapid quenching. 

I. INTRODUCTION 


Metallic glasses prepared by quenching molten alloys at the 

6 

rate of approximately 10 °C/sec have been intensively 
{ 1 — 6 ) 

studied, ' Such materials are characterized by amorphous 
st'.te, ductility, lower modulus of elasticity, high flow 
stress, superior fatigue resistance, higher electrical 
resistivity, superior corrosion resistance, etc. 
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Compared to metals, polymers exhibit not only relatively 
poor thermal conductivity but also a narrower temperature 
range for quenching. Accordingly, similarly high quenching 
rates are not possible for polymers. Using thin film 
specimens of glassy polymers and selecting adequate quench 
media, however, relatively fast cooling rates in the range 
of 500-700°C/sec can be achieved. According to the concept 
of fictive temperature (7/8) it is believed that the fictive 

temperature, where glass formation begins, will be greater 
at increasing cooling rates. (7-9) In this connection, 
excess free volume that is required for segmental motion 
will be frozen in. Subsequently, when such a quenched 
glassy polymer is aged at a temperature below T g , 
densification and change in the population of the various 
conformational isomers of the polymer chains will ensue in 
order to reach an equilibrium state. The reduction in the 
excess free volume, thus, results in a significant observed 
embrittlement of amorphous polymers. 

Theoretical and experimental results reported by Williams, 
Landel and Ferry[13] , Turnbull and Cohen^ 4 ) and Kovacs^ 1 ®) 
suggest that molecular transport mobility depends on free 
volume. The relationship between free volume, volume 
recovery, molecular mobility, and mechanical behavior can be 
illustrated by the following: 


FREE VOLUME MOLECULAR MOB ILITY-O MECHANICAL BEHAVIOR 
VOLUME RECOVERY 


Various studies on the subject have not adequately defined 
the effects of free volume on mechanical properties, 
especially for rapid quenching effects. More specifically, 
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the role of excess volume In the ductility of glassy 
polymers must be defined in order to gain further insight 
into the modes of molecular motion that contribute to the 
observed ductile behavior, Based on this understanding, it 
is hoped to predict the changes in material properties and 
further to utilize the benefits of rapid quenching effects 
for glassy polymers as engineering materials. Thus, study 
of rapid quenching effects on the mechanical properties of 
amorphous polymers has been undertaken, 

I I . MATE RIALS AND EXPERIMENTAL^ METHODS 

In order to achieve the highest cooling rates, thin film 
specimens of polyvinyl chloride (Geon 121) which contains 
no additives, are prepared either by solvent-casting in 
T.H.F. or by melt-molding. Typically, film thickness is in 
the range of 30pm to 130 pm. During annealing and quenching, 
the specimens were wrapped in aluminum foil (0,02 mm thick) 
and copper sheet (0.13 mm thick), in order to prevent 
surface contamination and to maintain good shape for 
mechanical tests. The temperature for annealing was in the 
range of 90°C to 160°C at atmospheric pressure. The time 
for annealing was 3-6 minutes. Quenching has usually been 
accomplished by plunging into ice water. 

Among the methods available for density measurements, 
hydrostatic weighing is the best approach for fast 
measurement and accuracy. The use of thin film specimens 
has two serious disadvantages: the weight is only 100 mg 

and the large surface area causes serious static electricity 
problems. Accordingly, to monitor density changes to 
+0.0001 gms/cm , not only a very high sensitivity balance 
but a big chamber to minimize static between the specimen 
and chamber wall are required. 

Following the ideas of Madorsky : ^ ^ ) a microbalance was 
constructed to perform the hydrostatic weighing. This 
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apparatus employs a quartz spring or a fine tungsten wire in 
the form of helical spring. In order to reduce the 
variability in the downward surface tension force, the 
suspension wire which has a diameter of 7 Mm was treated by 
electrical heating under l~10x!0 torr vacuum, yielding a 
uniformly granular coating of chromium oxide on the surface 
of the nichrome wire. Further, a precision thermoregulator 
was made to control the temperature of the distilled water 
within +0.005°C. Consequently, using a cathetometer which 
reads to 10 pm, the balance can detect weight differences of 
5x10 gms.,^giving the desired density sensitivity of 
0.0001 gm/cm or better. 


The volume recovery behavior under various temperatures is 

also studied by a glass dilatometer described in detail by 

Bekkedahl . ( 17 ) This consists of a reservoir, containing 

the specimen of known mass, joined to an accurately 

calibrated capillary of known cross-sectional area. The 

technique of filling mercury has been improved by employing 

a specially designed mercury reservoir. The specimen is 

then subjected to different thermal treatments by immersing 

the dilatometer in a glycerine bath, which is 

thermostatically controlled to within ±0,05°C in the range 

of 80°C to 130°C , and the water bath, which is controlled by 

the thermoregulator within ±0.005°C in the range of 20° to 

82°C, Finally, using a cathetometer reading to 50m, this 

— 5 3 

set-up can detect volume changes to 1x10 cm . 


Because of the thin film nature of the specimens used, it 
was necessary to construct a special apparatus to perform 
creep measurements. This required careful temperature 
control and optical determination of strain because a 
mechanical extensometer would disturb the response of the 
specimen. Stress-strain curves were determined using 
standard procedures. 


Ill, RESULTS 


Typical densif ication results determined by hydrostatic 
weighing are given in Pig, 1, These data show the effect of 
quenching from various initial temperatures to ice water. 
More of the excess volume is retained when quenched from 
higher temperatures. The densif ication due to physical 
aging takes place almost as a linear function of logarithmic 
aging time in each case, 

Dilatometry is used to provide volume relaxation data over a 

broader temperature range. To construct isochronal volume 

versus temperature curves it was necessary to carry out 

volume recovery at various final or aging temperatures. 

This is dilaometrically performed by quenching a PVC 

specimen from 99°C, where PVC is characterized by an 

equilibrium structure, to various temperatures below T , The 

volume recovery behavior in the temperature range of 2B°C to 

65°C, which is far below the glass transition temperature of 

PVC, is shown in Pig. 2. Similar results of volume recovery 

behavior near the glass transition temperature of PVC were 

also carried out. The relative excess volume 

(V-V )/V represents a dimensionless measure of the 
00 00 

departure from equilibrium. The straight lines of Pig. 2 
plotted vs. log time indicate a tendency for the rate of 
contraction to become slower as contraction proceeds with 
time. Furthermore, the relative excess volume retained at 
various temperatures is proportional to the temperature 
difference. 

Crossplots of the data yield two ischronal branches below 

the equilibrium state, at 1 minute and 200 hours after 

quenching from equilibrium at 99°C to various temperatures 

T, shown in Fig. 3. Both branches asymptotically approach 

two parallel straight lines with a slope of a = (d v/d 

T) » 7.35x10 cm /°C. These lines intersect the 
t 
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equilibrium volume lino nt Q.U1°0 and 0 0 , 9 0 c , which can bo 

associated with l ho gl ami transition temperatures T at 1 

g 

min and 200 h , respect ively. 

droop experiments nhow the oitoot oi excess volume cm the 

mechanical proportion. The ol l eel ot quenching from various 

initial temperatures in shown in Fig . 4. The quench from 

lf>)H' to 0°d produces the leant dense structure and the 

hiqhent creep valuer, whereas the quench from 9(>°C to Q°C 

qivon the most dense structure leading to the lowest creep. 

In essence, quench inq from various temperatures above T to 

q 

ice water clearly produces different levels of instantaneous 
creep as well as subsequent, creep at the same aq inq time 
prior to loadinq. At higher stresses, the effect of excess 
volume has been found to enhance the ductility while 
rod ue inq the yield stress. More high stress effects will bo 
studied in the future because of their implications for cold 
forming opera t ions , 

iv. com, w IONS 


Hydrostatic weighing and volume dilatometry results clearly 

demonstrate the more severe the quench, the greater is the 

departure from equilibrium volume, or the greater is the free 

volume content. Further, aginq of relaxing molecular 

rearrangements, based on the studies of phenomenological 

events, takes place as a linear function of logarithmic 

aging time at temperatures far below T . 

g 

Studies of mechanical properties monitored by tensile creep 
reveal the dependence of the creep behavior on the trapped 
free volume content. The phenomena of aging and quenching 
are essentially the r mo re v e r s idle. Further, a linear 
relationship of log (shift factor) vs, lug (aging time) is 
observed from Investigating the creep behavior of PVC film 
specimensi at various aging intervals after quenGhing, This 
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observation strongly corroborates the linear behavior of 
diminishing free volume. In addition, reduction in ductility 
is observed as aging proceeds. These observations indicate 
that the volume changes in the glassy matrix of PVC due to 
rapid quenching and subsequent aging are relatively small but 
they directly and strongly affect the mechanical properties. 
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PROJECT A5 : DEVELOPMENT OF RAPIDLY SOLIDIFIED OXIDATION 

RESISTANT ALLOYS 


Principal Investigator; Prof. G. J. Yurek 
Personnel; Mr. D. Cocke 

RESEARCH ABSTRACT 


The objective of this program of research is to develop 
iron-base alloys that contain low concentrations of expensive 
and/or strategic metals, but that are very resistant to 
oxidation at elevated temperatures. Rapid solidification 
processing is being employed to make alloys that are uniform 
in composition, that comprise very fine grains ( <1 0 urn) and 
that contain a uniform dispersion of very fine (<0.5 pm) 
nonmotallic inclusions. Each of these factors contribute to 
the oxidation behavior of alloys. The work in this program 
is directed to determining the compositional and processing 
variables that yield the optimum oxidation resistance. 

I. INTRODUCTION 


The oxidation of alloys in systems that operate at high 
temperatures limits their useful life by decreasing their 
load-bearing cross sections, by degrading their mechanical 
properties or by changing their shapes. The resistance of 
an alloy to oxidation at elevated temperatures is usually 
dependent on the formation of a thin scale of Cr 2 0 3 t 
A1 2 0 3 , or Si0 2 on its surface by reaction of CR, A1 , or Si 
dissolved in the alloy with oxygen or oxygen-bearing 
molecules in the gas phase. Continuous, pore-free layers of 
these oxides inhibit oxidation because diffusion of 
reactants through them is slow. Good resistance to 
oxidation in actual practice requires that protective oxide 
scales maintain good adherence to the alloy substrate during 
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observation strongly corroborates the linear behavior of 
diminishing free volume. In addition, reduction in ductility 
is observed as aging proceeds. These observations indicate 
that the volume changes in the glassy matrix of PVC due to 
rapid quenching and subsequent aging are relatively small but 
they directly and 3trongly affect the mechanical properties. 
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the application of external stresses and during thermal 
cycling. If a protective oxide scale spalls from the 
surface of an alloy, the alloy must be able to reform the 
protective oxide rapidly to avoid extensive oxidation of the 
bare alloy surface. 

The adherence of an oxide scale to an alloy may be Increased 
considerably by the presence of a fine, uniform dispersion 
of intermetallic or nonmetallic precipitates in the alloy. 

In general, the external protective scale tends to grow 
inward around the precipitates, thereby yielding intrusions 
of the scale into the alloy, which are effective in 
mechanically keying the scale to the alloy substrate. 

The selective oxidation of an alloying element to form an 
external protective scale requires the formation of a 
critical volume fraction of the oxide near the alloy surface 
to cause a transition from internal to external oxidation of 
the alloying element. A higher concentration and a higher 
diffusivity of the alloying element favors external 
oxidation as does a lower solubility and diffusivity of 
oxygen in the alloy. Grain boundaries and dislocations that 
intersect the surface of an alloy promote the formation of a 
continuous layer of the most thermodynamically stable oxide 
( Cr 2°3' a - l 2°3' or s i°2) on surface. 

Very fine-grained alloys may be obtained by rapid 
solidification processing. The sizes (<10 pm) remain stable 
at elevated temperatures as long as the alloy contains a 
dispersion of very fine stable nonmetallic inclusions to 
inhibit grain growth. The combination of a stable, fine 
alloy grain size and a fine, uniform dispersion of 
nonmetallic inclusions should provide an excellent alloy 
microstructure for oxidation resistance, provided the alloy 
composition is adjusted properly. 
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The purpose of this progra , is to use the principles 
outlined above to employ rapid solidification processing to 
develop iron-base alloys that are very resistant to 
oxidation at elevated temperatures and that contain very low 
concentrations of expensive and/or strategic metals. 

II. RESULTS 

The focus of the initial part of this program is the 
processing and testing of very fine-grained RS Fe-Al alloys 
containing 4-10 wt%Al and about 1 vol% MnS. The A1 is 
present in these alloys to form an external scale of 
Al^O^ that will inhibit oxidation of the Fe. The purpose of 
the MnS which is supposed to be present as a uniform 
dispersion of very fine MnS inclusions, is to maintain the 
fine grain size of the alloy at elevated temperatures. 
Manganese sulfide was chosen because it is a very stable 
compound and because the diffusivity of sulfur in iron-base 
alloys is low. These two factors inhibit coarsening of the 
fine inclusions in the RS alloy at high temperatures, 
thereby maintaining their effectiveness in pinning grain 
boundaries. 

The first rapidly solidified alloy produced in this program 
was an Fe-6wt%Al alloy. A master alloy was made by melting 
the elements in an Al 0^ crucible under an argon atmosphere 
in an induction furnace. Manganese and MnS were added to 
the master alloy in quantities that would yield a final 
composition of ~6wt%Al, ~1.6wt%Mn, and ~0.3wt%S. The final 
composition of the RS alloy was, however, 5.5wt%Al, 1.5wt%Mn 
and 0.075%S, thereby indicating that considerable sulfur was 
lost during the melting steps. 

About three pounds of the RS alloy were produced by 
remelting the master alloy in a calcia-stabilized zirconia 
crucible under an argon atmosphere and allowing a thin 
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stream of the liquid alloy to fall between two molybdenum rolls 
that were rotating eountor-ourrontly to each other, The rolls 
wore not cooled by external moans. The RS alloy so produced 
was in the form of flakes that were <vQ,5 to 2.0 mils thick, 

‘vo. 5" wide and ^2 inches long. The flakes were cut into smaller 
pieces before being compacted in a mild steel can, which was 
subsequently evacuated and welded shut boforo being hot extruded. 
The hot extrusion was carried out at 1100°G with a 20 ?! reduction 
ratio. The steel can was machined from the surface of the ex- 
truded rod to yield a bar of fully dense RS alloy that was 
*v0 . 5 inches in diameter. 

The as-cast master alloy and the RS alloy were characterized 
using optical /,nd scanning-transmission electron microscropy 
(STEM), The as-cast alloy comprised large grains vo.l to 1,0cm 
in diameter, The grain size of the RS alloy flakes depended on 
the thickness of the flakes, which reflects a variation of 
grain size with cooling rate? i.e., the thinner flakes com- 
prised fairly uniform grains that were *vl to Gym in diameter, 
while the thicker flakes contained both the fine grains and 
larger grains that were v20 to 40ym in diameter. 

The extruded RS alloy contained a mixture of fine grains (2- 
,10ym) and larger grains (20-50ym) . An etched transverse cross 
section through the extruded bar had a marbled appearance 
caused by selective etching of the alloy around oxide films 
that had existed on the surfaces of the flakes before hot ex- 
trusion. The presence of the oxide film is attributed to the 
rather slow rate of cooling of the flakes after exiting the 
Mo rolls in the rapid solidification apparatus? i.e., the 
flakes were hot for a sufficiently long time to allow oxidation 
by the small amount of oxygen in the argon atmosphere. The 
RS alloy also contained a number of inclusions of Mo (<V100- 
500ym in diamter) , which apparently were picked up from the 
Mo rolls. 
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Examination of the extruded RS alloy with a STEM revealed 
the presence of very fine nonmetallic inclusions (~0.05 to 
0.1 pm in diameter) , which were identified by EDX analysis 
to be mostly MnS. A number of the precipitates were 
identified as ZrO^ and (Mn,Zr)S. The source of the Zr was 
apparently the calcia-stabilized zirconia crucible used 
during the remolting process. Many of the precipitates were 
located along alloy grain boundaries, but once again, a 
considerable range of grain sizes was observed. 

A thermogravlmetric apparatus was constructed to conduct 

both isothermal and cyclic oxidation tests of the rapidly 

solidified alloys. Some preliminary isothermal oxidation 

tests were conducted with the alloys described above. The 

tests were run for lOOh at 900°C in pure oxygen at 1 atm 

pressure. Both the as-cast and the RS alloy exhibited a 

rapid period of oxidation in the early stages of exposure, 

followed by a period of slow oxidation that is attributed to 

the formation of an A1 0 surface film. The transition to 

2 3 

passive behavior occurred more rapidly with the RS alloy, 
which is expected because of its finer grain size and its 
more uniform composition. The RS alloy also exhibited, 
however, accelerated oxidation around the inclusions of Mo. 
These results, although encouraging, are only preliminary in 
nature and are not considered to be a test of the goals of 
this program because the microstructure of the RS alloy did 
not meet fully the specifications of the program. 

The problems encountered in the production of an RS Fe-Al 
alloy using uncooled, twin Mo rolls to achieve rapid 
solidification have prompted us to seek an alternative rapid 
solidification process. The current plan is to produce RS 
Pe-Al alloys using melt spinning with a water-cooled copper 
wheel. Such an apparatus is available in Prof. N. Grant's 
laboratory. We have designed and constructed a melting uni : 
that will be employed with the spinner to produce continuous 
ribbons of the RS Pe-Al alloys. 
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An A1 0 crucible will bG used in the melting unit, which 
2 3 

will eliminate contamination of the melt from the crucible. 
The melt will be forced out of the crucible through a small 
hole in its bottom by pressurizing the crucible with 
purified argon. The atmosphere above the crucible will, 
ideally, bG saturated with respect to sulfur, which will 
help to minimize losses of sulfur. 

The melt spinning apparatus produces a ribbon (~0.2 mil 
thick, ~12. 0 mils wide) that io both thinner and more 
uniform in thickness than the flakes produced by roll 
quenching, and the cooling rates are higher. Thus, the 
grain size of the RS alloy should be finer and more uniform. 
Also, because of the greater cooling rate, subsequent 
oxidation of the alloy directly after solidification should 
bo less, which will help to eliminate the presence of prior 
foil boundaries in the extruded product. 
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PROJECT A6: LASER MATERIALS PROCESSING FACILITY 

Principal Investigator: Dr. J. S. Haggerty 

RESEARCH SUMMARY 


The ceramics group within the Materials Processing Center and 
the Advanced Energy Materials group in the Energy Laboratory 
have jointly supported construction of a laser facility which 
is specifically designed for materials processing research. 

It gives MIT research capabilities which are unique in aca- 
demic institutions. 

The facility uses a custom designed 1500 watt C0 2 laser 
(Photon Sources Model 1033) with controls which permit CW, 
repetitive pulse, and shaped power-time cycles. The controls 
are designed to be interfaced with numerically controlled 
positioning equipment. The laser is designed to emit two 
600-650 watt Tni mode beams, two 600-650 watt Tm n , mode 
beams, one 1200 watt Tm 00 mode beam or one 1500 watt Tm^^ 
mode beam. In all configurations, the laser can be tuned 
to emit less than 20 watts. This power range will satisfy 
requirements from smaller diameter crystal growth up to 
modestly deep penetration welds in steel. The beams go from 
the laser to the experiment stations through conduits posit- 
ioned at ceiling level. Removable mirrors are employed to 
direct the beam(s) to the different experiment stations. 

The first experiment station consists of a machine designed 
for floating zone crystal growth. Pour beams are directed 
onto the heated region. We will also install anvils for 
quenching molten drops to study new glass and nonequilibrium 
crystalline compositions. The laser heat source permits high 
purities, freedom to select ambient atmosphere, and extremely 
high melting points. This chamber and optics can also be 
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used for powder synthesis at rates which are substantially 
higher than is possible with the optics and laser (150 watts 
maximum) we are presently using. 

The second experimental station will consist of a large 
controlled atmosphere chamber with provisions for manipulat- 
ing sample position. It will be used for a variety of 
cutting, welding, alloying, glazing, and heat treatment 
experiments, it can also be used for laser chemistry 
experiments . 

The crystal growth station has been operational since October 
1980. Since then it has been used to produce single crystals 
for several researchers at MIT and has undergone minor modi- 
fications to improve operational characteristics. Examples 
of crystals which have been grown include: •Ba x Sr 1 __ JC Nb 2 0, 
LaBg, LiF 2 , Fe 2 0 3 , Fe 3 C> 4 , FeO, Y 2 C> 3 stabilized Zr0 2 , high 
purity A1 2 0 3 / Ti0 2 do P ed Al 2 C> 3 , and Sc 2 0 3 stabilized Ta 2 C> 5 . 
These crystals were grown for Professors Bowen, Cannon, 
Tuller, Coble, Yurek, and Linz. The equipment has also been 
used for preliminary cutting and shaping experiments with 
Si 3 N 4 to displace diamond shaping processes. In addition to 
the research already completed, this laser heated crystal 
growth facility is the basis for several research proposals 
that have been submitted. 

The funds for the laser heated crystal growth facility came 
from a combination of four sources, including this NASA 
grant. 
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B. FLUID FLOW IN CRYSTALLIZATION PROCESSES 


Project Bl: FLUID FLOW IN CRYSTAL GROWTH: ANALYSIS OF THE 

FLOATING ZONE PROCESS 


Principal Investigator: Professor R, A. Brown 


Personnel: Mr. 

Mr . 
Mr . 
Mr . 
Mr . 


C . J . Chang 
H.M. Ettouney 
G»M, Harriott 
L.H. Ungar 
Y. Yamaguchi 


SUMMARY 


This research program is aimed at a fundamental under= 
standing of the interactions of heat, mass, and momentum 
transport in crystal growth from the molt. Emphasis has 
been on studies of the small-scale floating zone process and 
on a prototype of the vertical Bridgman growth system. In 
both systems detailed numerical calculations are being used 
to disect the interplay between fluid convection and dopant 
segregation. These calculations are based on newly 
developed finite-element techniques that make feasible the 
complete solution of solidification problems which include 
convection. Other numerical methods have been developed for 
solving thermal models of crystal growth processes with 
melt/gas menisci and for simulating the microscale instabil- 
ities in solidification interfaces. 

1. SCOPE OF RESEARCH PROGRAM 


Our research program is a comprehensive theoretical and 
computational study ultimately directed towards a 
fundamental understanding of the interactions of heat, mass, 
and momentum transport in crystal growth from the melt. 
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Although special emphasis has been placed on the floating 
zone process, the techniques for analysis and the physical 
insights developed in this research have broad application 
to other problems in melt crystal growth and studies are 
underway into Bridgman Growth, Edge-Defined Film-Fed Growth, 
and into the microscale instabilities in directional 
solidification. 

In each research project, a combination of analytical 
modelling and computer-aided calculation is being used to 
develop both a qualitative understanding of the relevant 
physics and quantitative models of the selected crystal 
growth processes. To do this, state-of-the-art numerical 
methods have beer developed for handling solidification 
problems. These methods give us the capability for accurate 
calculation of the fluid flow and mass transfer near 
solidification interfaces. 

The five research projects that are currently underway are: 

1) Thermal modelling of crystal growth processes with melt/ 
gas interfaces. (H.M. Ettouney) 

2) The coupling between rotationally-driven convection, 
surface-tension-driven convection, and solute transfer 
in the floating zone process at low gravity. (G.M. 
Harriott and L.H. Ungar) 

3) Fundamental studies of natural convection in low Prandtl 
number melts; effect of a magnetic field on the 
intensity of convection. (Y. Yamaguchi and C.J. Chang) 

4) Effect of buoyancy-driven convection on melt/solid 
interface shape and radial segregation in vertical 
Bridgman growth. (C.J. Chang) 

5) Microscale modelling of melt/solid interface dynamics. 
(L.H. Ungar) 


RESEARCH SUMMARY 


2.1 THERMAL MODELLING OF CRYSTAL GROWTH PROCESSES 

In principle, the macroscopic shape of the interface 
separating melt and solid during solidif ic-'tion depends on 
the transfer of heat, mass, and momentum. The conservation 
statements for the transport of these quantities through 
both melt and solid phases and across phase boundaries are 
mathemically coupled to the location of the phase boundary. 

In turn, the location of the melt/solid interface is set by 
local statements of conservation of energy and mass along 
the phase boundary, each written in terms of fluxes 
extrapolated from the bulk phases, and by assumptions of 
equilibrium for concentration and temperature. The first 
goal of this portion of our research has been to develop 
numerical methods for solving solidification problems 
accurately and efficiently. This goal has been accomplished 
for steady solidification problems and the finite element 
methods that are preferred for these calculations are 
outlined in a manuscript ( 3 ^ based on models that neglect 
fluid flow in the melt. The reasoning behind the use of 
these simple models is two-fold. First, they are good 
descriptions of heat and mass transport for many crystal 
growth geometries, especially in a low-gravity environment. 
More importantly, novel numerical methods are most 
efficiently developed for these less complicated models and 
then extended to account for bulk convection in the melt. 
This strategy has proven successful in the development of 
the finite element techniques outlined in Section 2.4 for 
treating natural convection near a melt/solid interface. 

The Isotherm-Newton formulation, the most efficient finite 
element technique developed by Ettouney^ 2 ^, has been 
extended to calculate simultaneously the temperature and 
concentration fields in a binary melt and to account for the 
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dependence of the melting temperature on concentration. 

Here the conservation equations for mass and energy are 
solved simultaneously with the condition for thermal 
equilibrium at the interface. Although this technique 
requires solution of a large set of nonlinear algebraic 
equations, it is not intrinsically limited to binary systems 
with small slopes of the solidus curve. These calculations 
have uncovered a primary difficulty in calculations of 
combined heat and mass transport in crystal growth systems. 
For liquid semiconductors, the region of variation of 
concentration in the melt is confined to a small layer near 
the melt/solid interface, whereas temperature variations are 
significant throughout both melt and solid. Concentration 
is best approximated by a finite-element basis that is 
concentrated near the solidification interface; however, 
approximations to the temperature distribution must be 
uniform. We have developed methods for simultaneous 
calculation of temperature and concentration on two 
different finite-element grids. These techniques establish 
the methodology necessary for detailed simulation of 
thermosolutal convection near a solidification interface. 

Ettouney has extended the finite-element algorithms to 
calculate the shapes of melt/gas menisci and melt/solid 
interfaces in crystal growth systems where both types of 
boundaries are present. The menisci are located oy solution 
of the equation of capillary statics and by satisfying the 
critical angle condition [Surek and Chalmers, J. Crystal 
Growth , 2 9, 1-11 (1975)] for steady crystal growth. The 
output of these calculations is the shape of the molten 
region along with the thickness of the grown crystal as a 
function of heat transfer conditions and growth rate. 
Applications to Edge-Defined Film-Fed Growth have defined 
limits in terms of growth rate and die temperature for 
steady operation. Similar studies of the floating zone 
geometry are underway. 
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2.2 FLUID MECHANICS OF A ROTATING LIQUID "ONE 

Research has Hocused on undorskandi.no the fluid mechanics of 
a liquid zone trapped between two cylindrical and 
differentially rotated solid rods and is aimed at describing 
the influence on the stability of the molten zone and on 
mass transfer of fluid flows driven by rotation of the food 
and crystal rods. G.M. Harriott has developed asymptotic 
methods for calculating these flows and the deflection of 
the melt/gas meniscus when the Reynolds number Re^iRH/v is 
small, where v is the kinematic viscosity of the melt, R is 
the radii of the rods, and H is the height of the molt; see 
Fig. 1. 

The details of this analysis are contained in a manuscript^. 
One major result shown in Fig. 2 is the cellular streture of 
the Clows produced by rotating the top surface at rate sR 
and the bottom surface at rate u. Depending on the relative 
rotations measured by s, and the aspect ratio of the zone 
AsR/n , there may be either one or two toroidal cells in the 
zone. Sample streamlines for a zone with aspect ratio of 
one are shown in Fig. 3. The value s-~l corresponds to 
exact counter-rotation. Here the stagnation streamline that 
separates the two cells is exactly a plane at the middle of 
the zone. For any other configuration this streamline is 
not a plane but a curve; see Figs. 3a-3f. 

Harriott's analytical results have been compared to detailed 
finite element calculations that are possible for a wide 
range of Reynolds number. The asymptotic results agree 
within ten percent with the numerical results up to Ro-50 
and hence an? very useful, estimates of the fluid flow in a 
micro zone . These flows have been used as the basis for 
calculations of radial dopant segregation under steady state 
cond i t ions . 
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Figure 2. Variation of cellular structure of flow as a function of 
aspect ratio A and relative rotation of the end rods s. 



- 2/3 



Figure 3. Streamlines predicted by analytical solution for a rotating liquid 
zone as a function of relative rotation rate s for A=1.0. 
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Numerical investigatlono of the flow at higher Reynolds 
number have uncovered multiple steady flows beyond a 
critical value that depends on aspect ratio and rotation 
ratios. For exact counter-rotation the onset of these 
multiple flows occurs at a bifurcation point in Re; for A 
equal to one, the critical value of Reynolds number Is 
approximately 109.5, As shown in Fig. 4, new families of 
flows with cells of unequal size start at this Reynolds 
number and evolve toward higher values of Re. The 
difference between these asymmetric flows and the original 
symmetric flow is brought out by the streamline plots in 
Fig. 5. Mathematical arguments indicate the stability of 
the fluid flow changes at the bifurcation point so that only 
the asymmetric flow would be observed experimentally. This 
prediction is being tested. 

In another study^), L.H. Ungar has mapped the stability of 
a rigidly rotating liquid zone held together bv surface 
tension for a wide range of the volume of molt, 
gravitational acceleration, and rotation rates. This study 
is built on earlier efforts at modelling the stability of 
floating zones (Brown and Scriven, Philo s. Trans, R. Soc.^ 
Land., 297, 51-79; Cor tell, Hardy, and Cordes, Colloid _ 
interface .Sc i §0 , } 26-136) and is based on the application 
of bifurcation theory for analyzing the existence of 
multiple steady solutions and their stability. The 
mathematical insights that were gained in studying this 
relatively simple problem have paid great dividends in the 
studies of the rotating flows and buoyancy-driven 
convection? see Sections 2.3 and 2.4 

2.3 FUNDAMENTAL STUDIES OF BUOYANCY- DRIVEN CONVECTION 

The complicated buoyancy-driven flow patterns in vertical 
Bridgman and floating zone growth systems have emphasized 
the need for understanding the transitions from axisymmotrie 
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Figure 4. Multiple solutions calculated by finite element analysis for s=-l 
and a b 1 .0, 


SYMMETRIC 
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Figure 5, Comparison of streamlines for flows in families of symmetric and 

asymmetric solutions. Cases (a-b) are for Re=125 and (c-d) for Re=200. 
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laminar flow to time-par iodic Clow in more simply defined 
heat transfer geometries. Our research has focused on 
axisymmetric flows in a vertical cylinder heated from below 
with either adiabatic or perfectly conducting sidewalls. 
This configuration was chosen for its geometrical 
similarities to the vertical Bridgman system and because 
results for the critical Rayleigh number for the onset of 
convection are also known (Charlson and Sani, In t. J. He at 
Mas_s Transfer 13 , 1479-1496 (1970 ) and JEnt, J. Heat Mass 
Transfer 14 , (1971). Experimental investigations for this 
same geometry are also possible as has been demonstrated by 
Olson and Rosenberger (J. Fl uid Mech . / 9 _ 2 , 609-6 29 (1980)). 


We have combined finite element methods with newly developed 

computer-aided methods for detecting and tracking multiple 

solutions to map out the steady axisymmetric flows that 

satisfy the Boussinesq equations for natural convection. A 

typical set of results is shown on Fig. 6 for a cylinder 

with height L twice its radius R and a conducting sidewall; 

here the Mussel t number for the top surface N of the 

u 

cylinder is plotted against the Rayleigh number 
3 

Ra= 3 g(T -T )L /av for each family of flow fields identified 
in our calculations. 

T 

The rest state corresponds to N =1 and is the stable state 

. . u 4 

up to the first critical value Ra=1.2xl0 where two families 
of singlecelled flows are stable and equally probable in 
experiments. Our predicted locations for the critical point 
is within three percent of the value calculated by Charlson 
and San.i using linear stability analysis. 


For a Prandtl number of one, the single-celled flows 
increase in intensity with increasing Ra until both families 
turn toward decreasing values of Rayleigh number at a limit 
point near Ra=6xl0 . This limit po’nt marks the loss of 
existence of a single-celled flow and may well mark the 
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Figure 6. Flow families for buoyancy-driven convection in a 
vertical cylinder heated from below and with a 
conducting sidewal 1 . 
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upperbound for steady stable ax.lsymmetr.ic flows. As 

depleted on Fig, 6, the two families of single-celled flows 

join families of flows with two rolls stacked axially in the 

cylinder. These two cell flows branch from the static state 

4 

at a second crxtxcal Rayleigh number Ra=3.1xl0 . As 
explained fully in a recent manuscript^), the flow 
structure of Fig. 6 is generic for Prandtl numbers between 
100 and 0,01 and for either insulated or conducting 
sidewalls . 

The detailed knowledge of the flow structure portrayed by 
Fig. 6 is a unique approach to understanding convection in 
melt crystal growth and is the basis for assessing the onset 
of fully three-dimensional steady and time-periodic 
convection, as well as determining precisely the effects of 
a magnetic field on the flow patterns in crystal melts. 


2 . 4 BUOYA NCY-DRIVEN CONVEC TIO N IN VERTI CAL-BRIDGMAN GROWTH 

Calculations of natural convection in melt crystal growth 
have been extended to account for the dependence of the 
shape of the melt/solid interface on the convective flow 
pattern in the melt. C.J. Chang has done this by combining 
the Isotherm-Newton technique developed by Ettouney with 
finite element techniques for calculating buoyancy-driven 
convection. Chang's algorithm is the most efficient one yet 
developed for solving solidification problems involving 
fluid flow and will oe widely used. The method has been 
reported in a recent publication^ for the solution of a 
prototype solidification problem where a melt and solid are 
held in a cylindrical ampoule in a vertically destabilized 
linear temperature field. Just as for a vertical cylinder 
filled with melt, multiple steady flows and melt/solid 
interface shapes are possible for this two-phase system. 

The structure of these flows is currently being determined. 

Pages 101-108 are blank 
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More important applications of Chang's algorithm are to 
studies of convection and mass transfer in the thermally 
stabilized growth experiments of C. Wang in the Materials 
Science Department at MIT, to the Bridgman experiments 
proposed for the Solidification Experiment System of NASA, 
and to a similar system being designed at MIT* A prototype 
for a Bridgman system is shown in Fig. 7 and sample 
isotherms and streamlines^for Prandtl number of 0.01 and 
Rayleigh numbers up to 10 are shown in Fig. 8 for a 
stationary ampoule. 

The me-lt/solid interfaces in these calculations are not 
flat, but slightly curved because of radial temperature 
gradients at the edges of the gradient region and because of 
convective heat transfer. Interfaces are shown in Fig* 9 
that correspond to the flow fields given in Fig. 8. These 
curved interfaces together with the convection in the melt 
cause radial segregation of the dopant in the crystal. 


Chang has calculated the solute distribution in a 

steady-state analog of Bridgman growth, using the velocity 

fields and interface shapes predicted from his convection 

calculations. Fig. 10 shows the distribution of solute 

along the solidification interface for Rayleigh numbers up 
4 

to 10 for the galliumdoped germanium system of Wang. The 
slight amount of radial segregation present in the absence 
of natural convection (Ra=0) is due to the curvature of the 


melt/solid interface as studied by Coriell and Sekerka [ J. 

Crystal Growth , £6, 479 (1979)]. As the convection level is 

increased, the delivery of solute to the interface becomes 

3 

more uneven radially until at Ra=10 the radial segregation 
of gallium is maximized. Further increases in the Rayleigh 
number causes the beginning of a solute boundary layer and 
decreases the amount of radial segregation . 
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figure 7, Prototype of vertical Bridgman system* 
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Temperature fields and streamlines for prototype of Bridgman solidification 
system with Pr=0.01 and Rayleigh numbers of 1, 1Q3, and 105. In each case 
the right side of the vertical slice through the melt is shown. 
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Figure 9. Melt/solid Inter- 
face shapes cal- 
culated with flow 
fields in Figure 

7 - 2 
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Figure 10. Radial concentra- 
tion of dopant in 
the melt to the 
melt/solid inter- 
face , 



RADIAL COORDINATE, r 




114 


2,5 MICROSCALE MODELLING OP SOLIDIFICATION DYNAMICS 


The Isotherm-Newton method tor solving steady solidification 
problems has been extended to account for the variation of 
melting point temperature with the shape of the solidifica- 
tion interface, where the condition for interfftcial 
equilibrium depends on the shape of the interface through 
its mean curvature. This added complication is system- 
atically handled because we write an explicit fbxbcdUelement 
representation for the location of the interface. Ungar is 
combining this new finite element formulation with methods 
for calculating the stability and multiplicity of steady 
solutions (see Ungar and Brown to analyze the evolution 
of microscale variation of aft interface that is constitu- 
tionally unstable. So far, the linearized results of 
Mullins and Sekerka [ J. Applied Physics , 35 , 444-451 (1964)] 
have been reproduced and steady interfaces with finite 
amplitude cellular morphologies have been calculated. 

Analytical methods developed in [7] have been applied to 
study an interface with a grain boundary. Contrary to 
previous analyses [Coriell and Sekerka, J_. Crystal Growth, 
19, 90-104 (1972) and 19, 285-293 (1973)] the presence of 
the grain boundary causes large changes in the amount of 
constitutional undercooling needed for an instability of a 
given wavelength. Ungar' s analytical results are valid only 
for grain angles near 90° (as measured through the solid). 
Numerical calculations are being developed that are valid 
for all angles. 

PRESENTATIONS 
■ " * — — * 

1. C.J. Chang and R. A. Brown, "Radial Segregation Induced 
by Natural Convection and Melt/Solid Interface Shape in 
Melt Crystal Growth." Fifth American Conference in 
Crystal Growth, San Diego, July 1981. 
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2, C.J, Chang and R.A. Brown, "Finite Element Analysis of 
Convection and Mass Transfer in Directional 
Solidification," Invited paper at the 18th Annual 
Meeting of the Society of Engineering Science, 
Providence, September 1981, 

3. C.J. Chang and R.A. Brown, "Finite Element Methods for 
Buoyancy-Driven Convection Near a Melt/Solid Interface." 
Second National Symposium on Numerical Methods in Heat 
Transfer, Maryland, September 1981. 

4. G. Harriott and R.A. Brown, "The Fluid Mechanics of a 
Differentially Rotated Captive Drop." 74th National 
Meeting of the AIChE , New Orleans, November 1981. 

5, Y. Yamaguchi, C.J, Chang and R.A, Brown, "Multiple 
Buoyancy-Driven Flows in a Vertical Cylinder." 34th 
Meeting of the Fluid Dynamics Division of American 
Physical Society, Monterey, November 1981. Abstract in 
Bull. Arner. Phys. Soc. , 26, 1286 (1981). 
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PROJECT B2 5 NUMERICAL MODELING AND OPTIMIZATION OF POLYMER 
MELT PROCESSING OPERATIONS 

Principal Investigator; Prof. David Roylanee 

Personnel: Mr. Brent Collins 

Mr. Craig Douglas 

RESEARCH AB5TRACT 

This project has been concerned primarily with the 
application of f inite~element computer analyses to polymer 
flows of the type encountered in melt processing operations. 
The finite element method is well suited to a wide variety 
of problem types, and is able to incorporate irregular 
boundary conditions and complicated fluid properties more 
conveniently than alternative methods. Our code as now 
developed is able to predict values of fluid velocity, 
pressure, shear stress, and temperature at any point within 
the flow field. As such, it is of value in diagnosing such 
processing problems as regions of fluid stagnation at which 
thermal degradation may occur, or regions of excessive shear 
deformation which lead to thermomechanical damage. It is 
further able to generate predictions of the forces which 
must be applied to the melt to achieve the desired flow, and 
this information is of value in designing processing 
equipment of optimal efficiency and minimum energy 
consumption. Ultimately, the project seeks to develop a 
versatile and easily implemented analytical tool which may 
be used routinely by processing engineers, as well as by 
researchers in rheology seeking solutions to presently 
intractable problems. 
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RESEARCH SUMMARY 

At the beginning of the current year, our finite element 
code had been developed so as to handle a wide variety of 
creeping flow problems, and some preliminary work on coupled 
temperature effects had been completed. The capability for 
coupled viscous flow-heat transfer problems was developed 
further during the current year, so that now many practical 
problems in nonisothermal polymer melt processing can be 
modeled numerically. Excellent agreement has been obtained 
between computed velocity and temperature fields and those 
available from theoretical closed-form solutions of selected 
problems. Finite element solutions have also been obtained 
in several cases in which closed-form solutions are not 
possible, such as the die entry problem, in all of these 
calculations, our code takes the momentum convection to be 
negligible (polymer melts have low Reynolds numbers), but it 
does have the capability to treat thermal convection. 
Convective terms are well known to introduce instability in 
computed results, and our code is able to avoid these 
problems either by successive grid refinement within the 
usual Galerkin method of finite element analysis, or by an 
"'optimal upwinding" scheme in which a greater weight is 
attached to the upwind side of each element. 

During the second half of the current year,, attention turned 
to transient problems, in which temporal derivatives appear 
explicitly in the governing equations, A capability for 
such problems has been implemented by means of a versatile 
and unconditionally stable two-point scheme, known as the 
theta-method. We have been able to predict transient 
flow-heat transfer solutions with this algorithm which agree 
well with closed-form results, but it must be mentioned that 
selection of the proper spatial and time steps is highly 
important; even a versatile algorithm may fail unless the 
problem's parameters are carefully chosen. 
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The code is presently written oo that the various material 
properties - the viscosity, conductivity, density, etc., - 
may be computed as functions of the solution variables 
(velocity gradient, temperature, pressure etc.). Only 
little work has been done in computing these nonlinear 
problems, however, and a more detailed exploration of such 
effects will constitute a major portion of the effort during 
the coming year. Nonlinear problems can be coded and 
attacked by a number of iterative schemes, but in general no 
single scheme is guaranteed of success,, One expects that 
numerical convergence will become more difficult as the 
strength of the nonlinearity increases, and it is likely 
that the very strong (usually exponential) dependence of the 
material parameters on the solution variables in polymer 
flows will prove a most difficult problem. 

In addition to the further exploration of nonlinear 
problems, the code will be developed so as to include the 
capability for handling reactive flows such as occur in 
reaction injection molding or epoxy dry-blend compounding. 
These will be both nonlinear and transient, so that many 
serious problems in convergence are expected. However, we 
are confident that many of these problems can be 
circumvented, and that a novel and versatile analytical tool 
will result. 

PUBLICATIONS ; 

1. Roylance, D. , "Finite Element Modeling of Nonisothermal 
Polymer Flows," to appear in the ACS Symposium Series, 
"Computer Applications in Coatings and Plastics," 

2. Collins, B . , "A Finite Element Model of a Nhite-Metzner 
Viscoelastic Polymer Extrudate," Technical Report 
T-733, Charles Stark Draper Laboratory, February 1981. 
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PROJECT B3 : STUDIES OF METALS ELECTROPROCESSING IN MOLTEN 

SALTS 

Principal Investigator: Prof. D. R. Sadoway 

Personnel: Mr. A. Abdelmassih 

Mr. P. T. Rogers 

RESEARCH ABSTRACT 

Molten salts are very important solvent systems in which the 
electrodeposition of a wide variety of metals may be 
conducted. However, solid electrodeposits from molten salts 
are typically incoherent, powdery, and/or dendritic. To 
understand better the electrodeposition process, fluid flow 
patterns in the electrolyte are being observed in order to 
determine how mass transport affects the morphology of the 
metal deposit. Studies are being conducted on the same 
metal, both in aqueous electrolytes in which coherent solid 
electrodeposits are produced, as well as in transparent 
molten salt electrolytes. Process variables such as current 
density and composition of the electrolyte are adjusted to 
change the morphology of the electrodeposit and, thus, to 
permit the study of the nature of electrolyte flow in 
relation to the quality of the electrodeposit. The results 
of this work will be helpful in electrochemical cell design 
and will serve as a data base for subsequent mathematical 
modeling of fluid flow in such systems. 

RESEARCH SUMMARY 


It is generally true that electrowon metal is of better 
quality than the thermal reduction product. However, solid 
electrodeposits from molten salts are typically incoherent, 
powdery, and/or dendritic. 
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The novelty of the present study is its consideration of the 
effects of fluid flow of the electrolyte in seeking to 
explain the poor quality of solid electrodeposits in molten 
salts. The plan is to electrodeposit the same metal in 
separate experiments from aqueous solutions and molten salt 
electrolytes. Transparent cells allow observation of 
electrolyte circulation patterns by a laser Schlieren 
optical technique. 

Zinc has been chosen as the metal to study. It can be 
deposited both from aqueous chloride and molten chloride 
electrolytes. Furthermore, both types of electrolyte 
systems are transparent to visible light, an important 
property which permits the use of noninvasive optical 
techniques of flow measurement. 

As this project is the first of its kind at the Materials 
Processing Center, much effort was spent in building up a 
laboratory for electroprocessing. This year the facilities 
for flow visualization by laser Schlieren photography were 
added. These include a pneumatically supported optical 
bench measuring 4x6 feet. 

Phase one of the zinc plating study was devoted to the 

electrorefining of zinc in aqueous acid chloride solutions. 

Working at room temperature without the need for furnaces 

and controlled atmosphere cells allowed one to concentrate 

on perfecting the optics of the system. Following a survey 

of the literature on acid zinc chloride solutions, a study 

of the chemistry of these solutions was conducted with the 

purpose of identifying the dominant ionic species present as 

a function of ZnC,n concentration, alkali chloride 

2 

concentration, and pH. Then, the electrochemistry of these 
solutions was studied by polarography and other techniques 
in various eoncentration-pH regimes which had been 
identified earlier. The optical apparatus to permit fluid 
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flow visualization by laser Schlieren and shadowgraph 
techniques was assembled and tested on simple electroplating 
cells. The present optical arrangement is capable of 
recording the laser Images in any of these formats: 35 mm 

photographic film, 8 mm movie film, or 1/2 inch 
videocassette. 

The results of the electrorefining of zinc in aqueous acid 
zinc chloride solutions can be summarized as follows. 

First, the laser Schlieren technique, while it does not 
provide quantitative information, is a very powerful tool in 
the study of fluid flow in electrochemical cells. The flow 
patterns observed under a variety of conditions show a 
significant amount of convection even at relatively low 
current densities. Such convection improves mass transport 
of zinc to the cathode well beyond the rate of diffusion, 
and thus, delays the onset of dendrite formation. This 
permits one to plate at higher deposition rates for a given 
surface quality. Furthermore, it is clear that the laser 
Schlieren technique can be used as a tool of analysis in 
electrochemical cell design. 

The next phase of the work is the adaptation of a cell to 
permit similar fluid flow studies of zinc electrodeposition 
in molten chloride solutions for comparison with the results 
of the aqueous system. 
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PROJECT B4 J HEAT P L OW CON T RO L A ND SEGREGATION IN 
I) I RECT I ONAL SOLIDIFICATION 


Principal Investigator : Prof. A. P. Witt 

Associate Principal Investigators s Pro! . W. M. 

Prof. P. K. 


Rohsenow 

Houpt 


Personnel: Mr. C. Wang 

Mr. T. Jan inski 
Mr. M. Wargo 


Mr. E. Bourret 


Mr. C. Herman 


RESEARCH ABSTRACT 


The thrust of this NASA-sponsored research program is 
directed toward the optimization of the vertical Bridgman 
technique for growth of electronic materials in single 
crystal form with the intent of advancing electronic 
materials processing at large and establishing a single 
crystal growth configuration through which the potential of 
reduced gravity environment for semiconductor crystal growth 
can be properly assessed. 

In the first phase of this research the limitations of the 
configuration were experimentally determined and heat 
transfer related deficiencies identified. Taking into 
account the necessity for complementary heat transfer calcu- 
lations in order to arrive at an optimized Bridgman-type 
grov,'th configuration, it was decided to base the design of 
an alternate system on the use of heat pipes separated by a 
gradient region since well-defined boundary conditions can 
thus be experimentally realized. Heat transfer analyses 
based on one- and two-dimensional models indicated the 
necessity of a flexible gradient zone configuration, a 
feature which was incorporated into the system presently 
under construction. 
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In a parallel study, directional melting of binary systems, 
as encountered during seeding in melt growth was analyzed 
for concurrent compositional changes at the crystal-melt 
interface. It was shown that steady state conditions cannot 
normally be reached during seeding and that the growth 
interface temperature (which determines the location of the 
crystal-melt interface) at the initial stages of seeded 
growth is a function of back-melt conditions. The theoreti- 
cal treatment was numerically applied to HgCdTe and Ga-doped 
germanium. 

Considering the increasing use of direct current flow across 
the crystal-melt interface for both analytical purposes and 
growth and segregation control, a theoretical and experimen- 
tal study of the thermal effects associated with current 
flow was conducted. It was found that experimental measure- 
ments of DC induced growth during crystal pulling can be 
used for the precise determination of the Peltier coef- 
ficient. The validity of the developed theoretical frame- 
work was verified by numerically solving the time dependent 
energy transport equations using an explicit finite dif- 
ference technique. The results were found to be in good 
agreement with experimentally determined current induced 
changes in the microscopic rate of growth. 


1. VERTICAL BRIDGMAN TYPE CRYSTAL GROWTH; A HEAT TRANSFER 
ANALYSIS 


This work, part of a comprehensive approach to the optimiza- 
tion of semiconductor crystal growth by the vertical 
Bridgman technique, is a heat transfer analysis based on 
one- and two-dimensional models. The present treatment 
neglects radial temperature gradients within the charge to 
yield relatively simple solutions of the axial temperature 
distribution. The analysis focuses on the axial temperature 
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gradient in the liquid at the growth interface and on the 
growth interface position within the furnace. Also con- 
sidered are the effects of charge motion, thermal coupling 
to the furnace, charge length, generation of latent heat of 
solidification, conductivity change at the growth interface, 
geometry and material of the crucible. 


Axial Temperature Distribution According to One-Dimensional 
Models 


The one-dimensional model used for the analysis of axial 
heat transfer is shown in Pig. 1. Making the assumptions 
described in Table I, the equation for the axial temperature 
distribution of the charge is the same as for a moving fin. 




4 BiG = -4Bi 0 

f 


[ 1 ] 


where 


C = z/D (non-dimensional axial coordinate) 

6 = (T-T ref )/(AT ref ) (non-dimensional temperature) 
0g = non-dimensional furnace temperature 
Pe = Peclet number = VD/a 
Bi = Biot number = hD/k 
a = thermal diffusivity 

h = heat transfer coefficient (charge-furnace) 
k = thermal conductivity of charge 
D = charge diameter 

V = mechanical charge displacement rate 


Equation [1] is solved for each zone using boundary con- 
ditions of equality of temperature and continuity of flux at 
the zone boundaries. 
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"Infinite 11 Charge Length 

The contribution to the solution of eq. [1] from the posi- 
tive characteristic root is normally small and is zero for 
an infinite charge length. If the charge is not "infinite" 
in length, the interface location will change as the charge 
moves within the furnace. When most of the charge io in the 
hot zone, the interface will be displaced toward the cold 
zone and vice-versa. The charge appears 


TABLE I 

ASSUMPTIONS FOR THE ONE-DIMENSIONAL MODEL 
(Equation [I] results from these assumptions.) 


• Temperature not a function of radius, 

• Hot and cold zone furnace temperatures uniform (heat pipe 
action) . 

• Gradient zone furnace temperature varies linearly. 

• Heat exchange between the furnace and charge is governed by 
a heat transfer coefficient, constant within each zone. 

® Transients are negligible, 
o No thermal convection in the liquid. 


"infinite" in length when the contribution from both the 
positive and negative roots become small. If the following 
criteria is chosen, 

exp r(mJGj<0.01 [2] 

(where C = length of charge in hot or cold zone for charge 
00 

to appear infinitely long) then, 
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If Pe is small, the characteristic roots become ; 

m ± « ±2 * / 5T [4) 

Substituting eq. (4) into eq. (3] results in a useful 

expression for determing £ . 

00 

5 

io„ > 2 Bi 

Motion of the interface may result in it being a location of 
unfavorable radial gradients. A method of stabilizing the 
interface location is to increase the length of the charge 
so that solidification in the center portion occurs with 
essentially infinite ends in the hot and cold zone. 

B -i — a. — . i *n -C r? — — a. 

AUU l\Umut;r &r.£C,’CG 


As Bi increases, the temperature follows more closely the 
furnace temperature distribution and the axial temperature 
gradient in the gradient zone increases. The relationship 
between Bi and axial temperature gradient can be 
demonstrated through the solution of eq. [1] which gives the 
axial temperature gradient zone as, 

|1 = - [X+(2Bi H 1 / 2 )“ 1 +(2Bi c 1 / 2 )” 1 ]" 1 [6] 

where : 

Bi H = Bi in hot zone 
Bi G = Bi in cold zone 

It is concluded that the establishment of large axial ther- 
mal gradients, frequently necessary to avoid constitutional 
supercooling, necessitate a large Bi. 
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The latent heat of nolidif ication enters the solution only 
in the boundary conditions at the interface. Assuming equal 
solid and liquid thermal conductivities, continuity of flux 
gives , 

,d6 t dO 

( d 5 #, liq e { d 5 } sol + Pe r h 

where i 

r H " ^ H / c p AT ref 

AH « latent heat of solidification 
The effect of the latent heat is small if, 

Pe R h << [X + (2Bi H 1 /2 ) -l +(2B i c V2)“l]“l. [8] 

Otherwise, the latent heat will increase the temperature 
everywhere in the charge and the interface will move toward 
the cold zone. 


Thermal Conductivity Change at the Crystal-Melt Interface 

The charge phase with higher conductivity will have lower 
axial gradients due to lower thermal resistance to heat 
transfer in the axial direction. 

The assumption of a large Bi forces the temperature drop 
from hot to cold zone to occur entirely in the gradient 
zone. The axial temperature gradient in the liquid becomes: 

d0 = 2 1 

( d? } liq \ * 1+R 

K 

where : 

r K = K liquid //R solid 
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In materials where ^ is greater than 2 or 3, a large reduc- 
tion in axial gradient in the liquid at the interface 
results. It is found that, for example in germanium, the 
difference in thermal conductivity between the solid and the 
liquid phase results in a reduction of 43% in the axial tem- 
perature gradient of the liquid in the gradient zone. 


Crucible Effect and Radial Resistance in the Charge 


Containment of the charge always decreases axial gradients. 

A crucible of low thermal conductivity adds thermal 
resistance between the charge and furnace effectively 
decreasing Bi; one of high conductivity "short-circuits" 
heat from the furnace within the crucible itself, again 
decreasing axial temperature gradients. The heat flow in 
the crucible are approximated by solving the heat conduction 
equation in the crucible: 


, 2 
1 JL 3Q_ 3 6 

r 3r Dr * + ,2* 

dt 


s s 


0 


[ 10 ] 


with the assumption that the 3 2 0/3s 2 term is equal to that 
of the charge in eq [1] and Pe approaches zero. 
Reformulating eq [1] results in a modified Bi number. 


Bi 


eff 


D 2 
1. C 


Bi 


Bi C lU q2 ' 


An 


k d 2 

•} "*!) + 


D -1 

k 1 C 

”("2 Bi c An D ~) ) 
c 

where : 

D c = crucible diameter 

k c = crucible thermal conductivity 

Bi C = h c D c /k 


Ul] 
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Bi e g f a Biot to bo used in eq. [1] 

h G b heat transfer coefficient between furnace and 
outer crucible surface 

k a conductivity of the charge (either solid or liquid) 


According to this relationship, the value of the conduce 
tivity ratio k c /k, which maximized the effective Bi is: 


Bi c * n ( D > 1/2 

rV— — > 

2((-~) 2 -l] 


( 12 ) 


For typical values of Pi G and r> G /P the optimum conductivity 
ratio is between 0.1 and 1.0. Increasing Bi^ and decreasing 
D g /D serve to increase the optimum conductivity ratio. 


An analysis similar to the one above is made for radial 
thermal resistance in the charge. The result of this analy- 
sis indicates that the Biot number should be modified as 
follows ; 


Bi 


eff 


Bi 



(13) 


where : 

Bi = Biot not accounting for radial temperature 
gradients in the charge 

Bi Q ££ » Biot accounting for radial temperature gradients 
in the charge 

It is found that if Bi/8 is small compared to unity, radial 
gradients are small enough so that eq. (1) reasonably pre- 
dicts the axial temperature profile. If Bi/8 is not small 
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compared to unity, the above relationship can be used to 
approximate bho effect of the radial gradients on the axial 
temperature distribution. 


2 . VERTICAL BRIDGMAN TYPE CRYSTAL GROWTH: AN EXPERIMENTAL 

STUDY OP GROWTH AND SEGREGATION 


The presently reported experiments were designed to 
establish the limitations of conventional, vertical Bridgman 
growth and to identify heat transfer controlled deficiencies 
of this growth configuration in conventional form. For this 
purpose, gallium-doped germanium crystals, approximately 7 
cm long, were grown in a vertical single -zone Marshall tube 
furnace. Two types of cooling were imposed at the bottom 
end of the crystal: direct cooling of the seed, for which 

the temperature at the bottom of the seed is expected to 
remain constant throughout the growth experiment, and 
indirect cooling, for which the temperature at the bottom of 
the seed is determined by the thermal environment, which 
changes continuously. 


Growth Behavior 


Depending upon the type of cooling, the microscopic growth 
rates are found to differ significantly from the lowering 
rate of the crucible and to never reach a steady state 
value, Consequently, during growth the position of the 
solid-liquid interface moved within the furnace. It is 
found that with direct seed cooling, the interface displace- 
ment (out of the furnace) for growth of 7 cm exceeds 8 mm. 
With indirect seed cooling, the growth interface displace- 
ment was in the opposite direction, i.e. into the furnace. 
The results for a lowering rate of 4 pm/s are summarized in 
Fig. 1, Heat transfer considerations support that an 
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Fig. 1: Top: Comparison of Growth Rate as a Function of Distance Grown for 

Different Seed Cooling Conditions. Bottom : Differential Change in 

Growth Interface Position in the Furnace as a Function of Distance 
Grown 
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increase in the rate of lowering increases the Peclet number 
(convective effect due to motion), thus increasing the tem- 
perature of the charge everywhere, and causing greater 
deviations from the initial interface position in the fur- 
nace. This effect was confirmed during a growth experiment 
with a lowering rate of 8 pm/s? the interface displacement 
was found to be 18 mm. The (experimentally determined) con- 
tinuously changing growth rates are interpreted primarily as 
transients due to end effects. 


Segregation Effects 

Anomalous segregation of gallium in germanium was observed 
in all instances where the microscopic rate of growth 
assumed very low values. During the initial stages of all 
seeded growth experiments, the concentration of Ga is found 
to decrease under a steadily increasing growth rate 
(Fig. 2). After approximately 50 pm of growth, the segre- 
gation behavior is observed to exhibit the expected rate 
dependence. Since seeding involves backmelting of the seed, 
this anomalous behavior was suspected to be related to 
varying seed concentrations. Consequently, seed con- 
centrations ranging from 0.1 to 10 times the equilibrium 
concentration for the solid in contact with the doped melt 
were studied. In all instances the anomalous gallium incor- 
poration was observed. 

The anomalous segregation behavior of gallium in germanium 
was observed also in all instances where growth was tem- 
porarily arrested and resumed after varying equilibration 
periods. For these experiments it was found through inter- 
face demarcation and spreading resistance measurements that, 
following the arrest, the microscopic growth decreased (as 
expected) and that initially the associated segregation 
increased (Ga concentration decreased). However, a gradual 
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Fig 2: Growth Rate and Segregation During Seeding. Note Period of 

Uncontrolled Growth and Anomalous Segregation at Oneset of Growth 







138 


ORIGINAL PAGE IS 
OF POOR QUALITY 


increase in Ga concentration was observed as the transient 
growth rate during the arrest period decreased to less than 
1 ym/s. Similarly, a decrease in Ga concentration was 
observed when growth was resumed and the growth rate 
increased from 0 to about 1 ym/s. As the growth rate 
exceeded about 1,5 um/s, the segregation behavior was found 
to be normal (Fig. 3). Experiments were made for varying 
lengths of arrest times, ranging from 5 minutes to 12 hours; 
in all instances, the same anomalous segregation was 
observed. The most conspicuous evidence of anomalous segre- 
gation is provided by Ga concentration maxima centered about 
growth rate convergences to zero. It should, moreover, be 
noted that the dopant concentration maxima assume values 
significantly in excess of values encountered during steady 
state growth at rates in excess of 6 ym/s . 

The presently reported "end effects" as well as the 
earlier uncontrollable, changing non-planar crystal-melt 
interface morphology, both of which have adverse effects on 
crystal growth and segregation, are heat transfer related 
deficiencies which appear to be inherent to growth by the 
Bridgman technique in conventional vertical configuration 
The elimination of these identified deficiencies have 
evolved as primary criteria in the conceptual design of a 
growth system currently under construction (Fig. 4). 

3. THEORETICAL ANALYSIS OF SEGREGATION DURING DIRECTIONAL 
MELTING OF SEMICONDUCTOR CRYSTALS WITH LARGE 
LIQUIDUS-SOLIDUS SEPARATION AND ITS IMPLICATIONS ON 
SEEDED GROWTH 


All seeded crystal growth from the melt involves as a first 
step a varying degree of backmalting of a crystal. This 
approach, generally successful in systems with small 
1 iquidus-solidus separation such as doped semiconductors, 
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must, however, be expected to be problematic if binary or 
multicomponent systems with large liquidus-solidus separa- 
tion (Hgj L _ x Cd x Te , Pbj__ x Sn x Te, Go-Si and the like) are 
involved. Any changes in composition that occur during 
melting at the crystal-melt interface are of fundamental 
importance since they will control the growth interface tem- 
perature and its stability. Both parameters affect the 
system's behavior during equilibration (melting arrest) and 
may control the degree of crystalline and compositional per- 
fection achievable during ensuing solidification. Seeded 
crystal growth experiments, for example by the Bridgman and 
Czochralski techniques, involve boundary conditions which 
were not considered in earlier treatments. It is the pur- 
pose of this work to analyze the evolution of the solute 
concentrations, C g (t) and C L (t), at a crystal-melt interface 
and the related interface temperature during melting, with 
emphasis on the initial transient behavior. The analysis 
applies to seeding during crystal growth of binary systems 
in a configuration where a charge of known composition C L is 
directionally melted in contact with a seed of composition 
Cg. The theoretical treatment is carried out in non- 
directional form; it is numerically applied to Hg^_ x Cd x Te 
and Ga-doped Ge to provide for comparative analysis of 
systems with small and large liquidus-solidus separation. 


Theoretical Analysis of Solute Redistribution During 
Directional Melting 

Making use of a series of basic assumptions, Fick's second 
law, formulated for a coordinate system which moves with the 
interface, assumes the form: 
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[2] 


Solute flux conservation through the interface yields: 
9C 


D 


s 3! 


ac 

RC* = D — 
S L Dz 


2=0 


R C 


2 = 0 


Using the change of scale: 

C = C - C 

o 

and the dimensionless variables: 


t > 0 


[3] 


t4] 


Rz 

Y m \ 


2 

R t 
D 


a 


D 


The system of eqs . [1] and [2] becomes: 


a 





3 2 C 3C DC 
L L L 

gy 2 + 3y = 3x 


[5] 


[ 6 ] 


This set of equations, with the appropriate boundary con- 
ditions and the initial conditions, can be solved using the 
Laplace transforms: 



and 


W 



D 
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yielding : 


\ r\f *>* 

d y dy 

s s 



5Y - 0 

C* 


>v »v 



with 

q ~ C 1) 
o 

Taking the general solutions: 





, -ay/a 
A o 


r*w 




+ 


q 


we obtain: 



142 

ORJO'NAL PAQEt 13 

or por^ oijAi.iTv 



and 



[ 10 ] 


These relationships depend mainly on the two parameters k 
and a = D S /D L . T!ie i nvei ^ se Laplace transforms are obtained 
in the general case by a numerical method. Analytical solu- 
tions were obtained for particular values of a. 

The inverse Laplace transform is made by approximating the 
integral : 

C+i«* 

— 1 / .AT t 

c ( Y r t ) = J e Y U)d4 

C-i 00 


by a trapezoidal rule using discrete values of the Laplace 
transform f(A) along the integration paths = c + iu>. For 
this purpose the integral [26] is approximated by Fourier 
cosine series which are evaluated using the Fast Fourier 
transform algorithm. The numerical inversion of [9] gives 
C s (y,T) for given values of a and k. Figure 1 is a plot of 
t 0 gg r the time necessary to reach 99% of steady state as a 
function of k, for four different values of a. 

[The real time (t) is related to the dimensionless time t 
through t = D l t/R^.] 



Functional Dependence of t q on k for Various Values 
of a. (Tq in non-dimensional form is the time re- 


quired for C* to reach 99% of the steady state value.) 
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At the hypothetical upper limit (a^l), the liquid and the 
solid behave symmetrically and a corresponding solution is 
given by; 


C, 


C 


/o=» r®» 

__ o y/r o y /t 

V y ' T) K T a ~ k)QVtc{ 2 rt ” T ) ex P ( ~y ) - 7” or£c( 77 - 7 


which reduces at the interface to; 


C * (0,T) b C (k-l)erf "7 
S o 2 


This relationship confirms the plotted data, 99% of steady 
state (erf /r/2 » 0.99) is reached for all values of k * 1, 
with t = 13.4. The corresponding curve on Fig. 1 is used as 
a reference line to determine the variation of t 0>99 with 
for any given value of k. 

Considering the special case of a+0, the inversion of [9] 
yields ; 


f L* 1 ft 

C S (Y,X) = ~ 2 ~ Hl-k)erfc(“^z “ ^ )exp(-y) + ■£ erfc("^“ + ~ 


2/7 


2 -k. , 2 -k 1 , , 1 -kx . _ . y 2 -k , 

— exp(— - “)yexp(“ 2 — ) erfc Hr + — /t } 
k 2 k 2 k 2/7 2 k 


[ 11 ] 


and 


C * ( 0 , t ) = “ {( l-k)erfc (“ ) + ~ erfc ( “] 
S 2 2 k 2 


( £_JS) exp T ) erf c(~~ t 

k Yr 


k 


[ 12 ] 
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The relationships [10] and [11] can be applied to solidifi- 
cation using the transform S defined by k 2*-l/k, S 4 -l , L*^S. 

It yields a solution which is identical with that given by 
Smith et a!, where is considered negligibly small. 

Figure 1 shows that t curves drawn for different values of a 
converge to aquation [10] for « £10“^. Generally accepted 
values for D s and D rj in semiconductors lie in the range 
lO"*^ - IQ" 8 cm 2 /soc for D s and SxLO” 4 - 10" 8 cm 2 /soc for D^» 
Therefore the Smith et al. solution (or its equivalent for 
melting [10] ) must be considered a good approximation. It 
should be pointed out. however, that there are as yet only 
limited data (of reasonable precision) for diffusion coef- 
ficients at melting point temperatures. While the uncer- 
tainty concerning the magnitude of is not expected to 

have a noticeable effect on the value of a, D g ca. *e con- 
siderably lower than assumed, as a consequence of the tem- 
perature controlled increase in vacancy concentration, for 
example, and may thus have a significant effect on the value 
of a. The data show that already for a - 10“ 3 the Smith et 
al. approximation no longer applies, and that for any given 
interface displacement rate, the compositional transients in 
the solid (t 0>99 ) (and in the liquid) increase with 
increasing a for k<2 and decrease for k>2. 


jation Behavior on Seeded 


'Stal 


Growth 


It is of interest to consider the present analysis in terms 
of real time (t) behavior, i.e. to consider, for example, 
the specific effect of the rate of directional melting (R) 
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on the time required to reach steady state. For this pur- 
pose, it is convenient to analyze the behavior of 

specific systems. An arbitrary choice is 
Hg i „ x Cd x Te of the composition x - 0.2 (C Q = 0.2), Pig. 2 
presents the concentration of CdTe in the solid (HgCdTe) at 
the crystal-melt interface (C s ) as a function of x for 
k s and k L . In real time (t), the CdTe concentration (C s ) is 
dependent on the rate of melting (R) which,, in order to be 
meaningful, must satisfy the stability criteria. At the 
beginning of melting, the system is at equilibrium [C g (z,0) 

= C o ,C L {z,0) - C 0 /k] , and therefore under conditions of 
absolute stability. As melting proceeds, the interface sta- 
bility is affected by out-diffusion and the subsequent for- 
mation of the boundary layer in the solid. Thus, while for 

C (0,z) = C /k melting can be initiated at any rate, out- 
L o 

diffusion and interface breakdown must be anticipated if the 
rate exceeds a critical value. If the stability criteria is 
strictly applied to a steady state melting configuration, 
the G s /R ratio required for stability is very high (=10° 
C/sec for Hg q gCd 0>2 Te )* Such a high value is achieved by 
either a very high thermal gradient or a very slow rate of 
interface displacement. Making use of Fig. 2 and applying 
numerical values for HgCdTe, the present analysis shows that 
under neither condition will steady state be reached during 
seeding in Bridgman or Czochralski growth where the seed is 
back-melted for several millimeters only; the crystal-melt 
interface temperature during melting must therefore be 
expected to lie somewhere between the solidus and the 
liquidus temperatures. Considering the back-melting of an 
ingot of Hg 0 gCd 0 2 Te for half an hour at the rate of 3.7 x 
10" 2 ym/sec (the rate required for absolute stability at 
steady state if Gg = 1000°C/cm), the value at the end of 
back-melting (i.e. after half an hour) is 4 x 10.“^ and, 
according to Fig. 2, the concentration in the solid (C g ) 
will be essentially unchanged at C c and the interface tem- 
perature will remain at about T g . 



T 


10 


20 


30 


T L 


. * 

Functional Dependence of Cg , the Interface Composition in 
Solid Hg q 8 cd 0 2 Te on T for ^ and (•) * 
assumed to be 4.7 x 10 ^ ; to assess the uncertainty associ 
ated with the assumed D value, the limiting case for 
a = 1 is also presented.) 
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Tt is noteworthy that for sy' ' ems with a liquidus-solidus 
separation smaller in h 9q . 8 ^ 0 . , ^ but larger than in very 

*k 

dilute systems, the temperature of the interface (T ) at the 
end of the backmelting is expected to be somewhere in bet- 
ween T s and T l . Its value, at a given time, can be calcu- 
lated from y curves of particular systems. The interface 
temperature is related to the composition at the interface 
through the relation T = T m + m fc C s , where T m is the melting 
temperature of the major constituent (for example, HgTe in 
Hg^ . 0 Cd 0 . 2 Te ^ and m t hhe sl °P e °f the straight line 
joining T m and the point on the solidus curve corresponding 
to the concentration Cg on the phase diagram. Therefore, 
the deviation of the interface temperature from the solidus 

•fa 

temperature during transient melting (AT = T - Tg) in 
dimensional notation is given by: 



m C 
t S 


m 


t=0 


C 

o 


fa 

As Cg approaches the value of kC Q , AT approaches a maximun 
value for which T* = T r . 


Segregation and Growth Interface Temperature 


As a consequence of segregation during transient melting or 

’ft 

solidif icaiton to growth interface, temperature (T ) is sub- 
ject to change and the interface displacement rate (R) 
varies with time and differs from the mechanical displace- 
ment rate of the system (Rgyg-j-) provided either through the 
translation of the charge relative to the furnace or through 
a shift in the imposed thermal gradient. A simple 
expression for R can be obtained in system with a constant 
axial thermal gradient (G= dT/dz) . R can be broken down 
into the two components: R = R S y St + R int' where 

R int ^eprecents the incremental rate due to the change in 
the interface temperature induced by the change in the 
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interface concentration [(3C/9t) z _g * 0] during the tran- 
sient. At any given time t, the interface temperature is 


T*(t) = T + mC*(t) [ A . 1 ] 

m S 

and its displacement rate of the interface due to its change 
in temperature is: 


dz 

- 12. 

dT 

R = 



int dt 

dT 

dt 


[A. 2] 


From [A. 1] : 


dT 

dt 



z=0 


and 


dz _ 1 
dT ~ G 


Therefore , 


R 

int 


m 

G 


dC 


dt 


z=0 


[A, 3 ] 


and , 


R 

int 


R 

syst 


dC 

m S 

G * dt 


z=0 


[A. 4] 


Equation [A. 4] shows that R = R s st if steady state is 
reached [(dCg/dt) =0], if G+°° or if m = 0 (no temperature 
change with composition). 
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4. THE EFFECT OF DIRECT CURRENT ON CRYSTAL GROWTH FROM THE 
MELT: InSb 


Experiments were conducted to investigate the effect of 
direct current on the growth of doped InSb. The results 
indicate that the growth rate is controlled on the 
microscale by the absorption of heat at the crystal/melt 
interface due to the Peltier effect. This effect is 
superimposed on longer time variations in growth rate caused 
by heat input imbalances due to Joule heating in the hot 
zone. A theoretical framework developed on the basis of 
conduction controlled heat flow predicts that the Peltier 
coefficient of InSb in contact with its melt can be calcu- 
lated from the current induced microscopic growth rate 
change. The results obtained using this technique are in 
agreement with the values of the Peltier coefficient calcu- 
lated (using the Kelvin relationship) from the published 
Seebeck coefficient. The theoretical framework was further 
verified by numerically solving the time dependent energy 
transport equations (with appropriate boundary conditions) 
using an explicit finite difference technique. The results 
are in good agreement with the experimentally determined 
current induced microscopic growth rate changes. 


The Effect of Long Duration (20 s) Current Pulses on Crystal 
Growth 


To determine the effect of current over long time periods, 
direct current pulses of 20 s duration Were passed across 
the liquid/solid interface at 120 s intervals. While appli- 
cation of direct current in the form of short pulses has 
been shown to enhance the microscopic growth rate (for 
Peltier cooling), with pulses of long duration the perturbed 
interface velocity must, with time, relax to the steady 
state pulling rate as the temperature distribution in the 
liquid and solid adjust to accomodate the new equilibrium 
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energy balance at the interface. Without flow of current, 
conservation of energy at the crystal/melt interface requires 

dT 

k — § 
k s ax 


- k 


dT 

I 

dX 


x=0 


*» pAHR 


x=0 


During growth with the flow of current: 


dT’ 
S 

dX 


x~0 


dT' 
L 

dX 


= ( pAHP ' + (IIJ) 


x=0 


where : 


ir = Peltier coefficient 
J - current density 


The primed values of the interfacial temperature gradients 
indicate that their magnitudes have changed to accomodate 
the new heat production (absorption) term irJ. 


Theoretical Framework for Calculating the Peltier 
Coefficient of InSb (in Contact with its Melt) from Current 
Induced Layer Thickness 

It has previously been shown by Lichtensteiger et al. that 
application of short duration (1 s) pulses of direct current 
across the liquid/solid interface of InSb during Czochralski 
growth caused an instantaneous increase in the microscopic 
growth rate which remained constant for the duration of the 
pulse. In addition, during current off growth conditions 
(3s), the microscopic growth uate also remained constant, but 
at a reduced magnitude. Their results were reproduced using 
experimental conditions that were the same as those employed 
to produce the previously discussed current induced layers: 
the graphite crucible was used with the temperature control 
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thermocouple located at the base of the crucible. The 
observation of constant growth rate in both current-on and 
current-off regions was verified for a variety of duty 
cycles (t on /t 0 ff^ and current densities. Interface demarca- 
tion (5 A, 10 ms duration, 100 ms interval) was superimposed 
on the current pulses. For values of current density below 
25.4 A/cm it was shown that the microscopic growth rate is 
effectively modulated between two values which remain 
constant for the duration of a pulse cycle. 


It could be shown experimentally and confirmed theoreti- 
cally, however that when a current pulse (Peltier cooling) 
is terminated, the liquid/ solid interface is repositioned 
to a point removed from its equilirium steady state (without 
current flow) . The result is a change in the temperature 
distributions in the liquid and solid which will influence 
the net heat balance at the interface and thus affect the 
microscopic growth rate. 


The change in growth rate, dR Q , caused by this change in 
temperature gradients (dG g , dG L ) is given by: 


dR = 


k dG 

s- 


k dG 
— L L 


pAH 


and is found to be -0.50 pm/s for one second pulses and a 
pulling rate of 12 pm/s. This effect is beyond the 
resolving capabilities of the interface demarcation tech- 
nique as applied here. For a single pulse cycle, the tem- 
perature gradients in the vicinity of the crystal/melt 
interface therefore remain virtually unchanged and it 
follows that the microscopic growth rate changes observed 
during current pulsing are the direct result of interface 
motion accomodating the absorption of heat due to the 
Peltier effect. The balance of heat at the interface in the 
absence of direct current is: 
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KG - k G « pAHR 
SB L L 1 

When current is applied, a term irJ, appears on the RHS, 
accounting for the Peltier effect: 

k G - k G - pAHR + )U 
S S L L 2 

(R^ and R 2 are the growth rates during current-off and 
current-on conditions respectively,) In the previous analy- 
sis it was shown that during a single pulse cycle the tem- 
perature gradients at the crystal/ melt interface remain 
virtually unchanged. Therefore, a direct relationship bet- 
ween the absorbed Peltier heat and the difference in 
microscopic growth rates (current-on vs current-off) is 
established : 

II J = pAH(R - R ) 

1 2 

Since the current density, J, is a predetermined parameter, 
it follows that the Peltier coefficient can be directly 
calculated from the experimentally determined, contiguous 
values of and R s : 

pAH( R = R ) 
n _ - 1 2 — 

J 

The absolute magnitude of the Peltier coefficient, as calcu- 
lated from this relationship, is directly dependent on the 
heat of fusion, H. The published values range from 67.7 J/g 
to 213.4 J/g. Consequently, the heat of fusion of InSb was 
measured using a Perkin-Elmer DTA 1700 differential thermal 
analyzer with a Perkin-Elmer System 7/4 thermal analysis 
controller. The heat of: fusion was measured to be 238.1 J/g 
(13,48 kcal/mole ) . 

The Peltier coefficient was calculated for contiguous pul$e 
sequences. The results are presented in Table I. 



154 


ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE I 

LISTING OF EXPERIMENTALLY DETERMINED 
VALUES OF THE PELTIER COEFFICIENT 


Sample 

t ON //t OFF 

B 

J (A/ cm2) 

y v) 

I-38a 

2/6 

18 

16.8 

-0.0822 

1-43 

2/6 

20 

13,0 

-0.0844 

1-64-10 

1/3 

10 

9.8 

-0.0968 

1-64-20 

M 

20 

19.7 

-0.0927 

1-64-30 

I! 

30 

29.6 

-0.0989 

1-64-40 

M 

40 

39.7 

-0.0871 

1-68-10 

20/100 

10 

9.5 

-0.1064 

1-68-20 

m 

20 

.•17.3 

-0.0960 

1-69-10 

20/100 

-10 

-12.2 

-0.1134” 

1-69-20 
[_ 

n 

-20 

-25.1 

Remelt 
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C. ADAPTIVE MATERIALS PROCESSING 


PROJECT Cl s ADAPTIVE CONTROL OF WELDING PROCESSES 


Principal Investigator: Prof. T. W, Eagar 



Personnel: Prof. J. Lang 

Dr. fi . Gershwin 
Mr. T. Lynch 
Mr. E , Raible 
Mr. K. Leong 
Mr. P. Kemp 

RESEARCH ABSTRACT 


Growth of automated welding is limited at present by 
applicable sensor technologies. This project involves study 
of new sensors for both arc welding and resistance spot 
welding processes. For arc welding, the arc itself is used 
as the sensor. Minor changes in both composition and 
geometry of the plasma result in fluctuations in the arc 
voltage. These fluctuations can be studied at frequencies 
from 50 to 20,000 Hz by recording the signal on audio 
equipment and later digitizing and analyzing the results 
with a computer. While the process as currently employed is 
an off-line system, development of an on-line process 
controller would be possible once the optimum analysis 
technique is developed. For resistance spot welding, it has 
been shown that the dynamic resistance of the weld joint can 
be correlated to growth of the fusion zone. Previous 
investigators have used similar techniques, but they have 
not analyzed the structure of the dynamic resistance curve 
in detail nor nave they correlated the results with actual 
weld geometry. Studies during the past year indicate that 
this technique is not as useful for projection welds as for 
spot welds. 
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1 , Arc Wo 1 fling 


This project was initiated during the post year. Ac a 
result much of the effort has been devoted to acquiring, 
setting up, and developing software for the new computer 
which is being used for signal analysis of the weld noise 
voltage. The system has both A/D and D/A converters capable 
of operation to 10 kHz. Paster sampling rates should be 
possible with further developments in the system software. 
Most of the sampling and graphing programs are now 
operational . 

Isolator and filtering circuits have been developed for 
interfacing the arc with the audio recorder and the 
computer. Work has begun on circuits which will allow 
direct recording of the arc to the computer, thus 
eliminating the need for the audio recorder. 

Much of the noise resulting from the welding arc is due to 
voltage fluctuations in the welding power supply. These 
fluctuations are uniform and reproducible, hence it should 
be possible to eliminate the power supply noise thereby 
improving the signal to noise ratio of the fluctuations 
resulting from the arc. A Thevenin equivalent circuit for a 
commercial welding power supply was developed and an attempt 
was made to incorporate this into a Kalman filter. The 
results were not encouraging so subtraction of signals using 
FFT techniques is now planned. 

During the past year it has been shown that welding over 
disturbances as small as 1.5 mm will produce measurable 
differences in the welding noise voltage. The optimum 
method for analysis of these fluctuations remains to be 
studied. Following this, the sensitivity, selectivity, and 
reproducibility of these disturbances will be investigated. 
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This work in being performed in conjunction with an office 
of Naval Research study of compositional changes in welding 
arc plasmas. That program has provided a transistorized 
welding power supply which will be used in future studies. 
This power supply will introduce much less extraneous noise 
into the arc voltage which should aid significantly in 
analysis of the data. 

2. Resistance Spot Molding 

During the past year the resistance spot welding machine has 
been instrumented with equipment for monitoring dynamic 
resistance, dynamic force, and electrode displacement. 
Acoustic emission equipment will be acquired in the near 
future. All of this information is being or will be used to 
study the mechanics of weld nugqet formation. Studies on 
projection wolds during the past year suggest that dynamic 
displacement may be the most useful technique for monitoring 
projection wolds. The equipment is being modified for 
studies in the coming year of spot welding of galvanized 
s tee Is • 

PUBLICATIONS 


1. E. Elias and T. W. Eagar, "Signal Analysis of Voltage 
Noise in Welding Arcs," to be published in Proc . of ASM 
Conference on Welding Consumables and Process Develop- 
ments , Peoria, IL, September 1981. 

2. J. G. Kaiser and T. V?. Eagar, "The Effect of Electrical 
Resistance on Nugget Formation During Spot Welding," 
submitted to the Weld. J. 
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XV, MATERIALS PROCESSING RESEARCH, MASSACHUSETTS INSTITUTE 
OF TECHNOLOGY 

This section summarizes ongoing Materials Processing re- 
search at Massachusetts Institute of Technology being 
conducted in a variety of Departments, Laboratories, and 
Centers. The research involves a total of 54 faculty 
members, 76 staff, 15 graduate students, and 57 undergrad- 
uates. The Table of Contents at the beginning of this 
report lists the projects and the page on which the research 
abstract will be found. 
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A. RAPID SOLIDIFICATION PROCESSING 



2V394 


Faculty: 


M 0 

M. 

No 

Ro 

R. 

Fo 

Ko 

Jo 

Jo 

G, 


Cohen Staff: L. Arnberg 


C. Flemings 

J. Grant 

Latanision 

Kaplow 

J* McGarry 

C. Russell 

Szekely 

Bo Vander Sande 
Yurek 


P. Domalavage 
Vr Franetovic 
A. Garrett-Reed 
Y. Gefen 
Mo Kurkela 
Hr, Ling 
R n Mandell 
J„ Megusar 
R. O' Handley 
G. Olson 
Y. Shiohara 
To Tsuru 
S Zhang 


SUMMARY 


Rapid solidification continues to comprise a major thrust 
area of Materials faculty and staff at MIT. Work is 
proceeding on crystalline and non-crystalline metals, and on 
polymers. Work is planned on ceramics. Research includes 
fundamental studies on solidification mechanism and 
resulting structure, innovative processing techniques, 
properties and applications. 

Professor Grant and co-workers continue their broad program 
focused on alloys which offer high strength at elevated 
temperatures, and on the processing methods necessary to 
achieve these goals. Professor Cohen and co-workers are 
concentrating on structure and structureproper ty relations 
in ferrous alloys. Professor Vander Sande, interacting with 
most other faculty working in this area, continues his work 
on structure. Professor Latanision is working on corrosion 
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related aspects and Professor Yurek is studying the 
high-temperature oxidation resistance of fine-grained 
alloys. Professor Flemings and co-workers are concentrating 
on relation of solidification theory to structures produced, 
on innovative processes for rapid solidification, and on 
solidification at high undercoolings. Professor McGarry has 
initiated studies on rapid < olidif ication of polymeric 
materials. Some current research activities are given 
below. 


1 . STRUCTURE AMD PROPERTY CONTROL BY MEANS OF RAPID 

SOLIDIFICATION TECHNOLOGY : GLASSY AND MICROCRYSTALLINE 

METASTABLE ALLOYS 

A rather broad range of programs constitutes the activities 
of this group. Emphasis remains focused on alloys which 
offer high strength at elevated temperatures and the 
processing methods necessary to achieve such goals. The 
application and study of rapid solidification (RS) to 
produce new classes of alleys and to achieve superior 
structures and properties are broadly applied to both 
microcrystalline and glassy alloys. 

Numerous microcrystalline alloy systems have been studied to 
date, including A1 , Cu, Fe , stainless steels, Ni-base super- 
alloys and Co-base superalloys. Glassy alloys are based on 
transition metal-transition metal combinations and on 
transition metal-nonne tal alloys. Phenomena of interest 
include strength, plasticity, toughness, corrosion resist- 
ance, superplasticity, fatigue improvements, irradiation 
damage resistance (fusion reactors), etc., etc. 

Brief paragraphs follow which attempt to attach a flavor to 
the specific programs underway. 
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Glass Formation, Alloying Effects on the Glass Transition 
and Crystallization Temperature of Ni-Nb Alloys 

Alloys in the 60-40 to 40-60 atom pet, range for Ni and Nb 
were carefully studied for embrittlement on heating and on 
crystallization, with measurement of attendant properties. 
Major alloying studies were undertaken to determine causes 
of increases or decreases in Tg, Tx and on toughness. Ease 
of glass formation through alloying and the effect of 
alloying on mechanical behavior were of particular interest. 

Stabilization and Strenqthen.inq of Pd Si Metallic Glasses 

J. A 8Q 2Q — 

The effect of composition on thermal and mechanical 
stability was studied by adding B and Zr to Pd Si singly 
and .in combination. Combined additions of B + Zr of 3 and 
1.5 at. pet., respectively, result in a stable glassy 
structure with an increased T g of 73 K. Additionally, 
high-temperature strength was significantly improved. 

Deformation Induced Dilatancies in Metallic Glasses at how 
Temperatures 


As a means of studying the deformation and fracture of glassy 
alloys, the compression mode was chosen since failure in the 
tension mode leads to infinitesimally small strain values. 
Using a PdSiCu alloy, it was possible to produce 2 mm 
diameter fully glassy rods which were ideal for miniature 
compression specimens „ Deformation values up to almost 75% 
were readily attained, without cracking, and with very large 
numbers of shear bands. The fine structures of these highly 
deformed specimens are undergoing further study. Large 
changes in free volume on crystallization have been observed 
as a function of the amount of prior deformation. 



iei 

a. ORIGINAL PAGSs J'-f 

0F POOR QUALITY 

Production and Study of Superfine Crystalline Alloys 
Produced by Conversion of Selected Glassy Alloys 

The attainment of grain sizes as fine as 100 to 500 A is 
demonstrated from prior glassy (homogeneous) alloys. 
Resultant grain sizes are very uniform. The aim is to 
produce ductile, tough, strong alloys which may also have 
other attractive characteristics such as corrosion 
resistance, for example. The alloys are of the 
Fe-Cr-Ni-Mo-B type with low boron content. 

The Structure and Properties of RS Lithium Modified 2024 
Aluminum Alloys 

Lithium contents up to 3 percent are readily possible in RS 
atomized alloys and result in fine-grained ( 2 to 3 microns), 
high specific modulus alloys. Studies include measurements 
of tensile properties, fatigue behavior, fracture toughness, 
crack propagation rates, etc. The ratio of Cu:L.i is a major 
alloying variable. 

The Structure and Properties of RS 2618 Aluminum Alloy 

This alloy, one of the best high-temperature ingot alloys, 
is being studied as RS atomized and hot consolidated 
material, and again after RS atomization plus attrition 
grinding to fine flakes to produce an oxide-dispersed 2618 
product, which is expected to have significantly improved 
strength properties from 150 to 25 0° C . 

The Structure and Properties of RS X2020 


The X2020 type of alloy ( Al-Cu-Mn-Li ) is recognized as 
having particularly attractive properties, such as a high 
specific modulus and high specific strength. This 
continuing study attempts to optimize such properties while 
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maintaining or improving such other features as notch 
toughness, notch fatigue behavior, etc. The Cu:Li ratio is 
a major variable since this ratio controls the amount and 
type of ternary intermetallic compounds formed. 

The Structure and Properties of RS Al-Mg-Li Alloys 

Whereas the Al-Mg-Li alloys have not shown comparable 
strength values to those of RS X2020 (Al-Cu-Li), the 
potential for major gains in specific strength and modulus 
favors the low specific gravity Al-Mg-Li alloys. Detailed 
studies of structure-property relationships are underway 
with special emphasis on notched strength values and 
fracture toughness behavior. 

The Structure and Property Relationships in X7091, 705 with 
Zr + N.i , and Other Alloys 

Maximum strength values are synonymous with the 7XXX alloys. 
Major additions of transition metals have been shown to be 
effective in RS 7XXX alloys by changing the chemistry of the 
phases formed and their size and distribution. In 
particular, additions of Fe , Ni , Ti , Zr , Mn , Co in various 
combinations and amounts have had beneficial effects 
especially on strength, with evidence that toughness and 
fatigue performance can be maintained at high levels, A 
number of these alloying combinations is being studied. 

High-Strength, High-Temperature, High-Conductivity Copper- 
Base Alloys 

Three distinct areas are being investigated: 

A. MZC Type Alloys ( Cu-Mg-Zr-Cr ) . These alloys with 

excellent strength at 20°C and excellent strength at 
400°C in creep rupture are being examined as conven- 



164 


ORIGINAL V, 

OF POOR QUALITY 


tional ingot product, RS atomized powders, and RS 
atomized and attrition milled flakes. The latter two 
conditions represent oxide-dispersed alloys at two oxide 
contents. Thermal conductivity values are near 80 
percent of those of pure copper. The presence of small 
amounts of ZrO^ 1:193 ^ een observed to prevent softening 
(at 20°C) after exposures as high as 1000°C for 1 hour. 

B. Cu-Ni-Ti Alloys . These alloys in ingot form are 
stronger than the MZC alloy, but have only about 60 
percent of the conductivity of pure Cu; they are also a 
bit weaker in creep at 400°C, As RS alloys, both in 
powder form and as fine attrited flakes, they have exhi- 
bited unusually attractive strength and ductility at 
20°C and outstanding creep-rupture properties at 400°C. 
Studies of thermomechanioally processed structures are 
continuing in a search for further improvements. These 
alloys with small amounts of TiO^ due to particulate 
surface oxidation do not soften at 20°C even after one 
hour at 1000°C. 

C . The Solution Treating, Ageing and Overageing Behavior of 
MZC and Cu-Ni-Ti Alloys Produced by Rapid Solidifi- 
cation . These copper alloys undergo complex ageing 
behavior patterns. The MZC alloy, with precipitates of 
both Cu^Zr and Cr , presents additional structural 
complication when ZrO^ and Cr^O^ are also formed during 
powder and flake processing. The Cu-Ni-Ti alloys, which 
harden unusually during quenching from the solution 
temperature, are further complicated when prepared from 
RS powders and attrited flake materials which produce 
TiO^. The presence of the fine oxide dispersions serves 
to pin dislocations and thus to retain the energy of 
cold work due to extrusion. The structures are thus 
fully stable up to 1000°C, do not show re-crystal- 
lization and maintain strength in spite of such high 
temperature exposure. 
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The Development: of Cu-lONl- and Cu~3QNl-Ba.se Alloys by 
Tertiary and Quarternary Alloying Additions 

These alloys,, produced by rapid solidification (ultrasonic 
gas atomization), are being studied to establish the role of 
alloying elements such as Fe, Cr , Ti , Mn and others on 
strength, toughness, salt water corrosion, etc, A number of 
the alloying elements enhance salt water corrosion as long 
as they do not form coarse excess phases. The RS approach 
is being used to enhance the control of structure. First 
results show significant increases in strength and ductility 
over current commercial ingot-based alloys. Thermomechan- 
ical treatments and corrosion studies are just getting under 
way. 

The Structure and Properties of Highly Alloyed Nickel-Base 
Superalloys Prepared from RS Atomized Powders 


The aim here is to determine if fully alloyed, superfine RS 
powders can be suitably attrited to near-micron' thick flakes 
to provide, for example, 1 tc 2 vol pet of .d^O in a 
complex alloy such as IN-100. Additionally, fully alloyed, 
RS IN-100 of very fine powder size is mechanically blended 
with fine yttria powders to attempt to produce a homogeneous 
oxide dispersion. The application, broadly, of O.D. type 
alloys is restricted by extreme costs and inability to 
readily reproduce the best reported properties of such 
alloys made by more conventional processing methods. 

The Superplastic Behavior of Duplex Ferrite-Austenite 
Stainless Steels 

This study has as its aims the determination of the optimum 
content of ferrite to austenite on superplast.ic deformation, 
on strength of the final product, on corrosion resistance 
and other properties. The role of composition is equally 
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Important since the various ratios of ferrite to austenite 
can be achieved by many alloying combinations. The limited 
solubility of the ferrite for carbon and nitrogen will per- 
mit precipitation of carbides and nitrides in the ferrite. 
The RS process guarantees a fine grain size (3 to 20 pm) in 
the extruded product. Thus a broad variety of structures 
and structural characteristics will be studied for their 
influence on superplasticity. 

The Initiation and Growth of Intercr ys talli ne Cracks in 6-6* 
Nickel-Base Superalloys, in < S-Carb,idc Cobalt-Base 
Superalloys, and in < S-Oxlde N.ickel-Base Oxide-Dispersed 
Alloys 


following the characterization of crack initiation and early 
growth in the creep mode, over a range of temperatures and 
stresses, samples will be prepared for fatigue testing as 
pre-cracked specimens to determine the further growth modes 
of cracks under a wide range of conditions of temperature, 
fatigue cycle frequency, hold times, etc. The aim is to 
determine which of these three basic alloy systems is better 
suited for further alloy development for service from 650 to 
1200°C. 

First-Wall Fusion Reactor Materials Prepared by RS 
Technology 

The extreme demands on materials to operate as first-wall 
materials with high dpa, helium generation and precipita- 
tion, radiation creep, etc., call for unusually stable, 
strong alloys. RS techniques can produce fine-grained, 
homogeneous alloys with selected second-phase precipitates 
or dispersions uniformly distributed in the matrix, and of 
controlled particle size and shape. Extensive alloying, far 
beyond that possible through ingot technology, has already 
been denonstrated . Selected experimental alloys suitable 
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for fission reactor exposure to simulate a fusion reactor .is 
another important issue. Using RS technology, a range of 
type 316 stainless steels with increasing Tie content have 
been prepared and are in test; type 316 S3 with a fine 
dispersion of AI 2 O 3 has also been prepared; ferritic Fe~BeO 
alloys, and selected copper-base alloys, have also been pre- 
pared and are undergoing neutron Irradiation and heavy ion 
irradiation. Further plans call for study of several RS 
ferritic alloys with and without oxide or carbide disper- 
sions . 

New Magnetic Metallic Glasses with Reduced Metalloid Con tent 


A series of new magnetic metallic qlassy alloys of the types 
COgQNbj^Bg, COg^Nb^gBg and Feg^NbgB-^ is being studied for 
nagnetic properties and stability against crystallisation. 
These alloys are observed to exhibit a larger correlation 
length (approx. 2.0A) which is shorter than is typically 
observed in "late transition 8 Q-motalloid 2 o" glasses 
containing more than one transition metal species (approx. 
2.3 to 2.5A). the magnetic behavior of these alloys will 
also be studied in the crystalline form, it is hoped as a 
ductile, formable shape. 

Magnetic Properties of Some Mew Co-Nb-B Metallic Glasses 

These new glasses show attractive soft ferromagnetic 
properties as do the more common transition metal/metalloid 
glasses, while retaining ease of fabrication and good 
thermal stability despite the significantly reduced 
metalloid content. One of the patterns observed is that the 
presence of boron in the nominally alloy increases 

the moment reduction due to Nb. Extensive studies of the 
magnetic moments and hybridization in Co-Nb alloys are 
continuing . 
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Sponsors : National Science Foundation (via the M.R.L. 

program at MIT) , Department of Energy 
(several divisions), Army Research Office, 
International Copper Research Association, 
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U.S. Navy (David Taylor Naval Ship Research & 
Development Center), Michelin Fellowship, Bowmet 
Fellowship, Cabot Fellowship, DARPA 

Faculty ; N. J* Grant 

Staff ; L. Arnberg, P, Domalavage, V. Franetovic, 

Y. Gefen, J. Megusar, R. O’ Handley 
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R. McCormick, C. Smith, E, Ting, W. Wang, 

W. Webster 

Publications ; 
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2. J. Megusar and N.J. Grant, "Stabilization and 

Strengthening of Pd Si Metallic Glass," Mats. Sci. 

BO 20 

and Eng. 49, No. 3 (1901) 275. 

3. J. Megusar, L. Arnberg, J.B. Vander Sande and N.J. 
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A Prime Candidate Alloy (PCA) by Rapid Solidification," 
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4. J.P. Clark, N.J. Grant and T.B. King, "The Market for 
Manganese Derived from Deepsea Nodules," Natural 
Resources Forum, 5, (1981) 249. 

5. J. Megusar, A.S. Argon and N.J. Grant, "Deformation 
Induced Dilatations .in Metallic Glasses at Low 
Temperatures," presented at 4th International Conference 
on Rapidly Quenched Metals, Sendai, Japan, August 1981. 
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6. N.J, Grant, J. Mogusar, and L, Arnborg, "Structure and 
Proportion of Rapidly Solidified Austenitic Stainless 
Steels for the Pus ion Reactor Environment ," presented at 
4th International Conference on Rapidly Quenched Metals, 
Sendai, Japan, August 1981. 

7. L> Arnberg, J. Mogusar, D. Xmeson, H.J. Frost, J.B. 
Vender Sande, Q.K. Marling, and N.J, Grant, "The 
Microstructure of Neutron Irradiated Rapidly Solidified 
Path A Prime Candidate Alloys," presented at Materials 
Research Society meeting on Rapid Solidification, 
November 1981. 


2 - RAPID SOLIDIFICATION PROCESSING OF IRON-BASIS ALLOYS 


The Office of Naval Research (ONR) has been supporting a 
research effort to determine the relationships between 
metallurgical structure and properties in rapid solidifica- 
tion processed iron-base alloys. This research program 
involves three faculty, one senior staff member, and four 
students . 

The research work being undertaken in this program can be 
divided into four parts: microstructurai characterization 

of rapidly solidified (RS) steels, grain growth retardation 
in RS steels, fracture toughness of RS steels, and decom- 
position of Fe-base amorphous alloys. Each of these areas 
is briefly described below. 

Microstructurai Characterization of Rapidly Solidified 
Steels 


In an effort to obtain a qualitative assessment of the 
fundamental solidification variables associated with solid- 
ification at rapid solidification rates, two austenitic 
steels have been observed in the form of powders and after 
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consolidation. The two steels, a high-sulfur, 303 stainless 
type and a high-phosphorous steel (Fe-Cr-Ni-P) have been 
prepared by the centrifugal atomization technique pioneered 
by Pratt and Whitney Aircraft, West Palm Beach, Plorida. 

Scanning transmission electron microscopy (STEM) has been 
used for microstructural and microchemical analysis of these 
steels. Both fee and bee are found to be primary solidifi- 
cation phases in the as-solidified powder, where the smaller 
powder particles (<70y dia.) tend to be bcc. In the high- 
sulfur steel, cellular solidification structures with 
sulfide precipitates (100 to 200 nm dia.) at the cell walls 
are observed in both fee and bcc particles. The bcc 
structure, however, has many small sulfide precipitates (10 
to 20 mm dia.) in the cell interior with few larger sulfide 
precipitates at the cell walls. The small precipitates, 
observed only in the bcc structures, form on cooling from a 
supersaturated solid solution that results from reduced 
solute partitioning during solidification. Partitioning of 
chromium and nickel is minimal in these cellular structures. 
A noncellular bcc structure is also observed with small 
sulfide precipitates throughout the entire structure. This 
non-cellular bcc structure results from smooth-front massive 
solidification. Analysis of the nucleation process for 
solidification indicates that a transition from fee 
nucleation to bcc nucleation occurs with increasing wetting 
angle in heterogeneous nucleation. Thus bcc should nucleate 
in the smaller droplets of a liquid dispersion where 
catalytic surfaces of low potency (large wetting angle) tend 
to be the only heterogeneous nucleants available. 

After hot consolidation, the high-sulfur steel contains a 
uniform dispersion of fine sulfides with an average sulfide 
size two to three orders of magnitude finer than conven- 
tionally-processed materiel. 
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In addition to the fee and bee primary solidification 
phases, rapid aolidif ieation of the high-phosphorus steel 
produces a Cine cellular phase at the cell walls. The 
amorphous phase, whleh is stab]e to ~5QQ°C, is enriched in 
phosphorus and chromium, but contains significantly less 
phosphorus than conventional glass-forming alloys. Hot 
consolidation of powders produces a chomically-uniform 
metastable austenite which can be effectively precipitation 
hardened by phospho-carbides. 

Grain Growth in Rapidly Solidified Steels 

Rapidly solidifed steols have been found to exhibit an 
unusual resistance to grain coarsening at high temperatures. 
Comparison of conventionally processed and rapidly 
solidified 9Ni-4Co and 2Mo steels revealed that the rapidly 
solidified material retains a grain size of ~20 pm at 1200°C, 
where conventionally processed material coarsens to several 
hundred microns. Thus far the only possible boundary 
pinning nechanism suggested by transmission and scanning 
transmission electron microscopy analysis is a fine 
dispersion of inclusions, principally MnS, Analysis of the 
grain coarsening behavior in terms of the theory of grain 
boundary pinning and particle coarsening indicates that the 
observed inclusion dispersions can account for the grain 
coarsening resistance. A thermodynamic survey of potential 
stable dispersed phases, based on both available 
experimental information and estimated subregular solution 
interaction parameters for austenitic iron, has identified 
TiN and rare-earth sulfides as particularly promising for 
alloy development via rapid solidif icaton. Application of 
TiN dispersions is being examined in rapidly solidified 
microalloyed steels. 



Fracture Toughness of Rapidly Solidified Steels 


The use of high austenitizing temperatures is known to 
increase the sharp crack toughness (K^ o ) of some steels, 
attributed primarily to dissolution of void initiating 
particles; simultaneously, t ne blunt notch toughness (C ) 
and ductility are found to decrease due to excessive grain 
coarsening. The high temperature grain coarsening 
resistance of the rapidly solidified steels described above 
may allow improvement of K without loss of energy and 
ductility. Preliminary toughness measurements on 
9Ni4C°-0.6C and 4M°-0.3C steels indicate that the KIC of 
rapidly solidified material is equal to or superior to that 
of conventionally processed material (after high austeni- 
tizing treatment). Furture efforts will focus on 4340, 

300M, and HP310 steels where the benefits of high austeni- 
tizing treatments are well established. Both I< and Charpy 
energy will be studied as a function of austenitizing 
treatment. The behavior of matrix tool steels (based on M-2 
and M-50 matrix compositions) is also being investigated. 

Decomposition of Fe-Base Amorphous Alloys 

The prospects for using the glassy phase of materials as the 
starting point for achieving (by thermomechanical treatment) 
technically advantageous phases and phase-morphologies which 
cannot be attained otherwise have provided the motivation 
for this work. Research has been focused on Fe-based 
alloys, especially Fe-B alloys, which are of intrinsic 
interest in the glassy state because of their unique 
magnetic properties. 

The as-quenched Fe-B glass is found to possess tetragonal 
Fe^B -type short range order by Mossbauer Spectroscopy. 
Stress relaxation in the glass has been studied by the 
changes of stress-induced magnetic anisotropy visible in the 



Mossbauer spectra. The anisotropy is affected by thermal 
relaxation and crystallization. 

The overall transformation kinetics have been studied by 
Differential Scanning Calorimetry (DSC). There are two 
crystallization steps involved in glasses with boron 
concentration less than 17%, while only one crystallization 
step is apparent in the glasses with boron concentration 
more than 17%. The primary crystallization of a- Fe in a 
8 4.5 b 15.5 9l ass was a -lso studied by Mossbauer Spectro- 
scopy and Transmission Electron Microscopy (TEM) . The 
growth of isolated Fe^B particles was studied in-si tu on the 
hot stage of a TFJM in the Feg^B^g glass. In low boron 
glasses, the crystallization of a - Fe is characterized by 
a diffusion-controlled mechanism (activation energy ~2 
ev/atm) with a limited particle size and decreasing 
nucleation rate. A transformation model with decreasing 
nucleation rate and limited particle size has been deve- 
loped. The polymorphous crystallization of FegB in low 
boron alloys is controlled by an interfacial mechanism. In 
high boron alloys the process is diffusion controlled with 
near simultaneous crystallization of Fe g B and a- Fe . 

From analysis of diffuse features in the electron diffrac- 
tion patterns of Fe^B crystallized from the glass, a defect 
structure is proposed, which is comprised of rod shaped 
domains of el structure in a matrix of e. While only the 
tetragonal forms of FegB are obtained by crystallization 
from the glassy state, the orthorhombic form crystallizes 
from the melt. 

a - Fe crystallized from the glass exhibits a smaller 
lattice parameter than pure a - Fe and, in the Mossbauer 
spectrum, evidence of a perturbed Fe- atom site. These 
effects are attributed to the metastable extension of 
solid-solubility of boron in a - Fe , reaching about 6% boron 
(in the Feg^ g ^15 S alloy). 
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3 . IRRADIATION EFFECTS IN RAPIDLY AND CONVENTIONALLY 
SOLI DIP Ili D ALLOYS 


Phase Stability In Rapidly Solidified N.i-Nb under Nl Ion 
Irradiation 


Two alloy compositions in the Ni-Nb system (Nig 0 Nb 40 and 
Nis^Nbig) were produced by rapidly quenching from the melt 
with the piston-anvil technique. The Ni 60 Nb 40 samples 
(originally amorphous) were transformed to a metastable, 
partially crystalline state by heat treatment in a differen- 
tial scanning calorimeter. The Nig 5 Nb 15 samples were fully 
crystalline, with the majority of the grains composed of 
collections of primary dendrite arms. Samples of both com- 

4 ”}* 

positions were irradiated with 4 MeV Ni ions in a Van de 
Graaff accelerator at temperatures from .48 to .62 T L for 
Nig 0 Nb 4 g (698 to 898 K) and at two dose levels: 2 and 20 

displacements per atom. 

The resulting irradiation induced microstructures were 
examined by transmission electron microscopy and compared 
with thermally aged samples. Under thermal aging con- 
ditions, the N.ig 0 Nb 4 g samples evolved towards the 
equilibrium twophase m.icrostructure (NiNb and NigNb) . Ni 
ion irradiation at 698 and 708 K produced a marked change in 
crystallite morphology, giving a greatly increased aspect 
ratio of the NigNb crystallites. The thermal evolution was 
arrested by ion irradiation in the temperature range 
studied, possibly by inhibiting the nucleation of the NiNb 
phase. No irradiation-induced voids were observed for this 
composition. 

Thermal aging of Nig 5 Nb 15 at up to 3 hour at 898 K produced 
no sifnificant changes in the microstructure. Aging for 
longer ties or at higher temperatures produced a complicated 
precipitate structure composed of large plate (rod) -like 
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Ni 3 Nb precipitates in a Ni matrix. The ion irradiated 
Nig 5 Nb 15 samples showed an increase in stacking fault and 
dislocation density, breakup of the grains into smaller 
collections of dendrite cells, dissolution and redistribu- 
tion of the ordered Ni 3 Nb precipitates and formation of 
oriented dislocation tangles in the dendrite cell walls. 

Ion irradiation under the above conditions was found to 
drive the microstructure along a different path than thermal 
evolution. 


Feasibility of Studies of Ion Irradiation-Induced Phase 
Stability Changes in Fe-Ni-C Alloys 

Neutron and ion irradiation experiments have shown that 
liquid metal fast breeder reactor and fusion reactor 
irradiations can be expected to produce changes in the phase 
stability of candidate alloys. The focus of this research 
is on the contribution of the chemical vacancy effect to 
phase stability. Excess vacancies are expected to contri- 
bute to the thermodynamic driving force for precipitation 
primarily when the microcluster rmatrix interface is par- 
tially or wholly disordered and able to furnish a vacancy 
sink and when the atomic volume of the precipitate is larger 
than the atomic volume of the matrix. The high vacancy 
supersaturations produced by irradiation may provide a suf- 
ficient thermodynamic driving force to induce precipitation 
at alloy compositions and temperatures at which a single 
phase is normally present. We are presently investigating 
this effect in Fe-32 wt % Ni-C alloys (0.03, 0.1, 0.3 wt % 

C) by irradiation of transmission electron microscopy (TEM) 
disks with high energy Ni ions. 

These Fe-Ni-C compositions were chosen for several reasons. 
First, for alloys containing 32 wt% Ni the martensite start 
temperature is approximately -100°C even at very low C con- 
tents and thus martensite formation (which might deform the 
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TCM disk surfaces, making the location of the narrow damage 
reg ion within the disk difficult), will bo suppressed. 
Second, austenitic steels are of practical interest for 
radiation environment applications and an Fe-32 wtt Ni-C 
alloy will remain austenitic at the temperatures of interest 
to this project. Third, in Fo-Ni-C alloys possible precipi- 
tates include graphite and (Fe,Ni) i C whose nucleation rates 
should be considerably enhanced by the excess vacancy eon- 
c q n tr a t i o n d u e to i t* r ad i a t i o n . 

Preliminary calculations have shown L*hat at 500°C a vacancy 
supersaturation of 1x10° may bo obtained with a displace- 
ment. rate of lxl o’"* 5 s’* 1 . Further nucleation rate calcula- 
tions indicate that the thermodynamic chemical vacancy 
drivinq force under heavy ion irradiation conditions may 
provide a sufficient driving force tor nucleation in the 
single phase austenitic region of the Fo^Ni—C phase diagram. 

Sponso rs: National Science Foundation through CMSE, 

National Science Foundation 
Facul ty: K. C. Russell 

G r a dju a t e S t u d e n b s : S, Boat, R. Cher nock, C. Parker 


4. RAPID RATH SOEIMFICATION PROCESSING 
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These programs are on (1) development and application of 
solidification theory to rapid solidification and solidifi- 
cation at high undercoolings, (2) experimental studies on 
structures obtained by rapid solidification, and relation to 
theory, and (3) development of innovative processes for pro- 
ducing rapidly solidified materials. 

The research is sponsored by National Aeronautics and Space 
Administration, National Science Foundation, and Army 
Materials and Mechanics Research Center. 
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Rapid Solidification of Binary Alloys 

Simple binary alloys are being rapidly solidified by splat 
cooling, ribbon quenching, and laser or electron beam 
surface melting. Cooling rates and interface growth 
velocities are calculated and structures produced related to 
solidification theory. As example, the amount of residual 
second phase is measured in rapidly solidified 
aluminum-copper alloys, using the SEM and quantitative 
metallography? also by using x-ray methods. Composition 
profiles in other binary alloys are measured using STEM, 
Results are compared to the developing theories of 
solidification of rapidly solidified alloys. 

Undercooling and Rapid Solidificatio n 

One Part of this study comprises dispersion of droplets of 
molten metal alloys in fluxes and slag to obtain large 
degress of undercooling following the methods of Turnbull 
and Perepezko. Methods are being studied of combining the 
large undercoolings obtained with rapid solidification. 
Structures are examined by electron microscopy and related 
to solidification theory. Work to date is on both low 
melting point alloys and on high temperature iron and nickel 
base alloys. A second part of the study involves levitation 
melted samples which are undercooled by varying and 
controlled amounts in inert gas, or in a surrounding fluid 
slag. The levitation work is being conducted on nickel 
base alloys. 

Sponsors : National Science Foundation, National Aeronautics 

and Space Administration, Army Materials and 
Mechanics Research Center 
Faculty : M. C. Flemings 

Staff : G. Abbaschian, Y. Shiohara 

Graduate Students : J.T. Burke, G. Chu, R. Ewasko, 

D. Maclsaac, L. Masur, Y. Wu 


179 


Publications : 

1. M.C. Flemings, "Segregation and Structure in Rapidly 
Solidified Cast Metals," Metallurgical Treatises , 
Co-editors, J.K. Tien and J.P. Elliott, Metallurgical 
Society of AIME , 1981. 

2. L. J . Masur and M.C. Flemings, "Solute Redistribution in 
Rapidly Solidified Al-Cu Alloys," Presented at 4th 
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Sendai, Japan, August 1981, 
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5. OXIDATION AND CHEMICAL STABILITY OF RAPIDLY SOLIDIFIED 
MATERIALS 


A broad range of programs are underway in the corrosion 
laboratory on corrosion resistance, hydrogen embrittlement, 
chemical stability, and oxidation behavior of rapidly 
solidified crystalline and non-crystalline materials. 

Corrosion Resistance of Microcrystalline Fe-Base Alloys 

In some instances microcrystalline alloys produced by rapid 
quenching are found to have superior corrosion resistance 
when compared to otherwise equivalent wrought or cast alloys 
of the same composition. During the past year we have 
extended our study of microcrystalline stainless steels 
having grain diameters ranging from 2 (microcrystalline) to 
100 ym (conventional wrought material). Alloys, produced 
with the aid of Professor N.J. Grant, were either fully 
austenitic or of an austenitic/ferritic duplex microstruc- 
ture. The m.icrocrystalline alloys are found to be con- 
sistently more resistant to localized corrosion in acid 
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chloride environments than their nought counterparts. 
Motallographic examination indie a on that pitting initiates 
more uniformly on mieroerystall ino surfaces, but pits do not 
propagate. VJo believe that the high degree of chemical 
homogeneity characteristic of the microcrystalline alloys, 
which resulted from /apid solidification, is the principal 
reason for their resistance to localized corrosion. 

Hydrogen Permeation and Embrittlement in Metallic Glasses 

Metallic glasses are considered to represent a good 
structural and compositional analog of grain boundaries in 
polycrystalline metals and alloys. Structurally, grain 
boundaries and glassy alloys may both be modeled atomically 
in terms of the packing of polyhedral units. Likewise, 
transition metal-metalloid type glasses are similar in 
composition to solute segregated grain boundaries. During 
the past year we have continued our examination, by means of 
electrolytic hydrogen permeation experiments, of the 
absorption and transport of hydrogen in Ni^P^ and 
NiSi 4>5 B 2 >8 glasses. These glasses approximate grain 
boundaries in polycrystalline nickel and nickel-base alloys 
(Incenel 600, Hastelloy C-276, for example) in which P 
segregation has been observed by Auger electron spectroscopy 
through the use of the Central Surface Analytical Facility 
in CMSE. We had observed that such glasses absorb con- 
siderably more hydrogen than is typical of bulk poly- 
crystalline material because, we believe, of the effect of 
metalloids such as P in stimulating hydrogen absorption. 

More recent work shows evidence of hydrogen trapping in 
these glasses when the charging overpotential exceeds a 
critical value. 

Chemical Stability of Metallic Glasses 


Glassy metals have remarkable mechanical, magnetic and 
electric properties, and, accordingly, great potential 



181 


commercial utility. Without chonH^M stability in the 
environments to which the no alloys may he exported while in 
service, their unusual properties may not, however, he 
appreciated. Home metallic q lasses are remarkably inert 
chemically. The objective of this program io to examine the 
corrosion resistance oL’ such alloys and to develop a 
lundamental basis for understand inn ibis behavior. In 
particular, studies are underway to evaluate the role of 
structure (crystalline vs. glassy), homogeneity, and alloy 
composition (the largo concentration of metalloids in some 
inert glasses, the role of film forming alloy additions, 
etc.) on the kinetics of anodic and cathodic reactions on a 
variety of glassy surfaces. During the past year wo have 
examined film growth kinetics on Fe-Cr-Ni-P-B and Fe-Ni-P~B 
glassy alloys by -cans of ellipsometry . Likewise, the 
composition and valence states of the elements present in 
thin passive films have been characterized by means of Auger 
electron spectroscopy and ESCA, We have observed that P is 
present in the passive films. This is in contradiction to 
earlier work in Japan which sgggests that P stimulates 
initial reactivity of the surface and subsequent formation 
of a protective film but is not itself incorporated .into the 
film. At any rate, matters appear rather complex and our 
current effort is oriented toward separating the effects of 
metalloid elements, such as P, from the influence of film 
formers, such as Cr , wit’ regard to passivation and film 
growth kinetics. Work is in progress to examine resistance 
to localized corrosion and stress corrosion cracking. 

A Study of the Relat i onship Between the Microstructure and 
Chemical Stability of Rapidly Quenche d Ir on Base Alloy s 

The principal objective of this research is to examine the 
possibility that through rapid quenching technology and 
subsequent microstructural manipulation, microcrystall ino 
iron-base alloys with improved corrosion resistance might be 
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developed* The alloyo of Interest at this stage include the 
Fq-'H and Po-Ti-P systems. Research to just underway. 

T he Oxidation Behavior of. Fine-Grained Rap idly So lidified 
Iron-Base Alloys 

This program, which began on 1 September 1981, is directed 
to determining the oxidation behavior of very fine-grained 
alloys produced by compaction and heat treatment of rapidly 
solidified alloys. The rapidly solidified alloys will be 
made by both melt spinning and gas atomization. The 
advantages and disadvantages of these two techniques will bo 
evaluated relative to the structure and the properties of 
the compacted alloy. Compaction will be accomplished by hot 
extrusion. The alloy systems to be investigated will be 
based on Fe-Cr-Al, with and without additions of Kare-ear th 
elements (e.g., Y) and with and without dispersions of very 
fine (<0.5 ym) nonmetallic phases (e.g., TiC, AIN, etc.) 
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8th InH. Congreas on Metallic Corrooion, p. 400. 
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4. R, M. Latanioion, M. Kurkelo and F. Lee, "The Role of 
Grain Boundary Chemistry and the Environment on 
Intergranular Fracture", Proc. Third International 
Confer once on Effect*? of Hydrogen on Behavior of 
Materials, in press. 
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Transactions , in press, 
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Structure of Grain Boundaries in Inconel 600", 
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6. OTHER RELATED PROJECTS 


1 . Rapid Solidification of Thermoplastics 
F. J. McGarry 

See — Effect of Processing on Polymer/Composito 
Structure an-i Properties, page 192. 

2 . The Modelling of Rapid Solidification Processes 
J. Szekely 

See — Mathematical and Physical Modelling of 
Materials Processing, page 236. 



185 


\ N82' 2 7 396 

B . E FFECT OF PROCESSING ON POLYMER/COMPOSITE STRUCTURE AND 
PROPERTIES 
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SUMMARY 


After many decades of regarding polymer processing as a 
series of unrelated, often empirical, disciplines dealing 
with the formation and shaping of materials, the industry 
has come to realize that further advances in the vitality 
and economic health of the field will require a more 
consistent and rational point of view in which processing is 
viewed as a participant in the underlying triad of relation- 
ships which comprise materials science and engineering. 

This triad includes processing as it influences material 
structure, and ultimately properties. Although several MIT 
faculty are presently involved in research dealing with 
various aspects of what might be termed polymer processing, 
a group of faculty has formed under the auspices of the 
Department of Materials Science and Engineering whose work 
is centered specifically around the theme of processing- 
structure-properties, and this group has developed enough 
coherence to justify its inclusion in this report under a 
common heading. Their work spans the full range of 
subdisciplines of polymer science and engineering: Professor 

Wnek is concerned primarily with polymer chemistry and 
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synthesis, Professor Sung treats structure modification and 
optimization through processing as the v nifying theme of her 
research; Professor Roylance has been developing methods of 
melt flow modeling as an important aspect of his work in pro- 
cessing structure-property relations of polymers; Professor 
McGarry’s longstanding interest in mechanical properties of 
composites has always included processing as an important 
aspect; and Professor Rose has extended his expertise in 
biomedical materials research to include polymer processing 
effects. Although any classification of the diverse research 
interests of this group is of course somewhat arbitrary, they 
are in fact unified by a common interest in processing and 
its role in shaping structure, and the following descriptions 
are intended to summarize their principal goals and approaches. 

An analysis of the design technology of advanced graphite/ 
epoxy composites being conducted by Professors Mar, Dugundji, 
Pian, and Crawley of the Department of Aeronautics and 
Astronautics is also reported herein. 


1. POLYMER SYNTHESIS 


Electrically Conductive Polymers 

Composites of polyacetylene and low density polyethylene 

(LDPE) have been prepared by impregnating the LDPE with a 

Ziegler-Natta catalyst followed by exposure to acetylene 

gas. The resulting materials soften upon heating provided 

that the polyacetylene content is <ca. 5 wt% and remain 

flexible and tough indefinitely. The composites can be doped 

-1 -1 

with iodine to afford conductivities as high as 10 ft cm 
Their conductivities decay less rapidly upon air exposure 
compared with ^-doped polyacetylene film. The construction 
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o £ batteries from the composites having open circuit volta- 
ges of <ca. 2.'2V has been demonstrated. The microwave 
absorption characteristics of the composites are being 
investigated . 

We are currently exploring the generality of our synthetic 
approach with regard to the types of polymeric materials 
which may act as media for the in-situ polymerizations of 
acetylene. The properties of composites prepared by simply 
dispersing polyacetylene powder in various matrices are also 
being studied. 

Polymeric materials are also being synthesized in which 
rigid, conjugated pendant groups are attached to an alipha- 
tic backbone. The goal is to produce materials which are 
processable (at least before doping) and which can be doped 
to high conductivities with both electron doners or accep- 
tors. The potential influence of ordered domains derived 
from the rigid pendant groups on electrical properties will 
be addressed. 

Spons or : MIT Center for Materials Science and Engineering 

Faculty ; G. E. Wnek 
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Hirsch, "Electrically Conducting Acetylene-Methylacety- 
lene Copolymers - Synthesis and Properties," 
Macromoleoules , 1_4, 479 (1981). 

3. G. E. Wnek, J. Capistran, J. C. W. Chien, L. Dickinson, 
R. Gable, R. Gooding, K. Gourley, F. E. Karasz, C. P. 



Li Ilya and K.-D. Yao, Adv. Chem, Se ries/ Las Vegas ACS 
Meeting, in press. 

4. J. C. W. Chien, G. E, Wnok, F. E . Karasz, J. M. 
W*:rakomski, L. C. Dickinson, A. J. Heeger and A. G. 
MacDiarmtd, "Electron Paramagnetic Resonance Saturation 
Characteristics of Pristine and Doped Polyacetylenes ," 
Macromo lecules , in press. 

5. M. E. Galvin and G. E. Wnek, "Electrically Conductive 
Polymer Composites. Polymerization of Acetylene in 
Polyethylene," P olym er, in press. 

2 • PROCESSING -STRUCTURE _ RELATIONS IN PO LYMER S 

Characterization of Polyme r Surface Struct ure 

In order to characterize polymer surface structure as a 
function of different processing conditions (injection 
molding, compression molding, etc.), last year we developed 
a modified attenuated total internal reflection FT-IR 
dichroism technique by using a symmetrical double-edged 
internal reflection crystal. We have now extended our work 
to measure three-dimensional orientation and other struc- 
tural parameters, by rotating the polarizer and further 
improving the optical arrangements on the ATR set-up for 
increased sensitivity and accuracy. By using uniaxially 
orientated polypropylene and polyethylene terephthalate 
films, we demonstrated that we can obtain not only three- 
dimensional orientation functions for each absorption band 
(crystalline as well as amorphous), but also other surface 
structural information such as degree of crystallinity of 
trans/gauche conformer content. 

When the surface of the polymer is rough or brittle, FT-IR 
ATR technique is difficult to apply. In collaboration with 



Dr. Krishnan at Digilab, we developed an FT-IR pho- 
toacoustic (PA) dichroism technique to determine molecular 
orientation from such surfaces. In order to demonstrate the 
feasibility of this new technique, we used uniaxially drawn 
PET and compared PA and ATR dichroism techniques. For three 
well-known parallel bands, the dichroic ratio is greater 
than unity in both ATR and FA, as expected from the drawing 
process . 

Sponsor : Office of Naval Research 

Faculty : C, S. P. Sung 

Graduate Student : J. P. Hobbs 

Publications : 

1. C. S. P. Sung, "A Modified Technique for Measurement of 
Orientation from Polymer Surfaces by Attenuated Total 
Reflection IR Dichroism," Macromolecules, 14, 591 
(1981) . 

2. J. P. Hobbs, C. S. P. Sung, K. Krishnan and S. Hill, 
"Three Dimensional Characterization of Polymer Surface 
Structure by Modified Attenuated Total Reflection IR 
Dichroism," Proceedings of the 40th Annual Technical 
Conference, Soc. Plastics Eng., May 1982 (San 
Francisco ) . 

3. K. Krishnan, S. Hill, J. P. Hobbs and C' S. P. Sung, 
"Orientation Measurements from Polymer Surfaces Using 
Fourier Transform Infrared Photoacoustic Spectroscopy," 
Applied Spectroscopy , May-June, 1982. 


Mechanistic Studies of Adhesion Promoters in A1 2 0 3 /Poly - 
ethylene Joint 

Several spectroscopic techniques (FT-IR, ESCA, SEM and 
SEM-EDX) were employed to characterize the molecular struc- 
ture of a deposited organo silane-type adhesion promoter and 
to probe the interactions at the polymer/silane interface as 
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well as at the sapphire/silano .interface. This structural 
information was correlated to properties of the dry joint by 
evaluating joint strength and failure mode and locus. 

Results are; first, direct confirmation of the heteroge- 
neous, multimolecular nature of deposited silane on 
sapphire; second, a strong possibility of interdiffusion at 
the polyethylene/silane interface; third, the possible pre- 
sence of a covalent bond at the sapphire/silane interface; 
and fourth, reasonable correlation with joint strength and 
failure mode. 

Spon s ors : Air Force Office of Scientific Research, Army 

Research Office 

Faculty: C. S. P. Sung, N. H. Sunq (Tufts University) 

Graduate Students: I. J. Chin, A. Kaul (Tufts University) 

Publications : 

1, N. H. Sung, A. Kaul, S. Ni, C. S, P. Sung and I. J. 

Chin, "Role of Organo Silanes and Organo Titanate in 
Promotion of Adehsion Strength of Al 2 03 /Polye thylene 
Joint", Proceedings of 36th Annual Technical Conference, 
Reinforced Plastics/ Composites Institute, SPI, Session 
2-B , 1981. 

2. A. Kaul, N. H. Sung, I. Chin and C. S. P. Sung, 
"Mechanism of Adhesion Promotion Through Organosilanes 
in a-Al 2 0 3 / Polyethylene Joints," Proceedings of 37th 
Annual Technical Conference, Reinforced 
Plastics/Composites Institute, SPI, Session 2-E, 1982. 

Molecu l ar U nders t and in g_ _ o f Physical Aging Phenomenon via 
Azochromorp horic Lab elling Stud ies 

As far as physical aging is concerned, the processing of 
plastics by extrusion, injection molding, etc. may be 
regarded as a quench. Therefore, the mechanical properties 
of the glassy plastics are expected to become more brittle 
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as a function of time. Wg use an azochromophoric label in 
the main chain of amorphous polymers in order to provide 
some insight into the physical aging phenomenon. Analysis 
of photochemical transacts isomerization kinetics of such 
azo labels reveals that kinetics can be characterized by two 
rate constants, The fast rate is similar to the rate 
observed in dilute solution, while the slower rate is a 
hundred times slower. As physical aging proceeds, the frac- 
tion (a) which is characterized by the fast rate decreases 
until it reaches equilibrium. We propose to interpret a to 
reflect the fraction of free volume above a certain size 
necessary for the azo label to isomer ize , and thus the 
decrease in a as a function of physical aging means the 
collapse of free volume of such a size. To the best of our 
knowledge, this is the first demonstration of the use of 
molecular labels by which we can estimate a fraction of free 
volume above a critical size at a given temperature and time 
of aging. This provides more information about the distri- 
bution of free volume and its change with aging. 

Sponsors: NSF Polymers Program, Whitaker Health Sciences 

Fund, Sun Kyong Fibers, Ltd. 

F aculty : C. S. P. Sung 

Graduate Students: L. Lamarre, K. II. Chung 

Public ation s : 

1. C. S. P. Sung, L. Lamarre and K. II. Chung, "Studies of 
Thermal Transitions and Physical Aging Monitored by 
Azochromophoric Labels Attached .in the Mainchain of 
Amorphous Polymeric Solids," ACS Polymer Preprints , 
22-2, 277 (1981). 

2. C. S. P. Sung, L. Lamarre and K. H. Chung, "Use of 
Azochromophor ic Labels as a Molecular Probe of Physical 
Aging in Amorphous Polymeric Solids," Macromole cules , 

14, Nov-Dec (1981). 
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3 • PROCESSING OF THERMOPLASTICS 
C old Form i ng of Po 1 y v inyl Chi geld e 

Polyvinyl chloride, one of the largest volume plastics, is 
especially sensitive to heat so molt forming operations 
require close controls to avoid thermal degradation of the 
polymer. The situation is exacerbated with formulations 

intended for higher temperature applications, therefore 
great benefit would derive from forming methods which could 
be carried out economically at or near room temperature. 

The purpose of this research is to identify and explore such 
methods suitable both for powder and sheet PVC. 

Considerable success has now been achieved in powder forming 
of homogeneous PVC samples with promising mechanical proper- 
ties . 

Sponsor : B. F. Goodrich Company 

Faculty : F. J. McGarry 

Staff : J . F. Mandell 

Graduate Stude nt: J. H. Chen 


Rapid So l I d i f lea tlon of Therm oplastics 

The free volume content in glassy thermoplastics depends 
upon the rate of cooling through the glass transition tem- 
perature. It strongly affects certain macroscopic mechani- 
cal properties via the process known as aging. Very rapid 
cooling and solidification of metals has produced amorphous 
structures with unusual properties. The purpose of the pre- 
sent research is to look for similar effects in glassy poly- 
mers containing abnormally large free volume. Some 
exploratory work with crystalline polymers may be done also. 
Careful measurements on quenched thin films of polyvinyl- 
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chloride and polystyrene do show a significant increase in 
free volume a a compared with control specimens, 

^onsors National Aeronautics and space Administration 

Faculty: F. J, McGarry 

Sta ff : J, F. Mandoll 

Graduate Students : A , Ag r a wa 1 , H , Le o 

Numerical Model ing of Polymer Melt Processing, 

Finite element techniques are being developed to model the 
velocity, stress, and temperature fields which exist during 
such polymer melt processing operations as extrusion and 
injection molding. Computer codes have been written and 
verified for steady and transient nonisothormal flows of 
Newtonian and shear-thinning fluids in which the effects of 
viscous heat generation and heat transfer by conduction and 
convection are incorporated. Research is presently ongoing 
aimed at extending the models to include chemical reaction 
and free-surface effects. The ability to model these flows 
in an economic and accurate manner will be of considerable 
value in designing processing operations so as to produce 
polymers of optimal properties and in developing equipment 
of improved efficiency. 

Sponsors: National Aeronatics and Space Administration, 

Army Materials and Mechanics Research Center, 
Draper Laboratories, Hysol Corp. 

Faculty : D, K. Roy lance 

Gradua te Students: E. Bronaugh, C. Douglas 

Publ icat ion : 

1. D, K. Roylance, "Finite Element Modeling of 

Nonisothermal Polymer Flows," to appear in ACS Symposium 
Series, "Computer Applications in Coatings and 
Plastics . " 
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Processing of Ultrahigh Molecular Weight Poly ethyl ene 

The performance of ultrahigh molecular weight polyethylene 
of total joint replacements depends crucially on how the 
reacted polymer is processed into the solid implant com- 
ponent, as processing determines the ultimate molecular 
weight and defect distribution. 

Evaluation of performance depends on accurate total joint 
simulations and assessments of wear and mechanical and 
chemical breakdown, and comparisons with clinical 
retrievals, This is done in several ways with the polymer 
processing as the independent variable, 

Sponsor: NIH 

£^ulty: R, M. Rose, A. M. Crugnola (University of Lowell) 

Staff: G. Arndt, E . Goldfarb, I. M. Puffer 

Grajluate Student: A, Litsky 

Undergraduate Students: A. Casavant, N. Goldberg 

Publ icat ions : 

1. R. M. Rose, A. M. Crugnola, W. R. Cimino, and M. D, 

Ries, "The In-Vivo Performance of Polyethylene 
Components of Total Joint Replacements," Proceedings of 
Implant Retrieval Conference at the National Bureau of 
Standards, May 1-2, 1980, 

2. R. M. Rose, II. J, Nusbaum, H. Schneider, M. Ries, I. 
Paul, A. Crugnola, S. R. Simon and E. L. Rad.in, "On the 
True Wear Rate of Ultrahigh Molecular Weight 
Polyethylene in the Total Hip Prosthesis," J. of Bone 
and Joi nt Surg., Vol. 6 2-A , No. 4, June, 1980, pp. 
537-549. 

3. R. M. Rose, A. Crugnola, M. Ries, W. Cimino, I. Paul, A. 
and E. L. Rad in, "On the Origins of High In-Vivo Wear 
Rates in Polyethylene Components of Total Joint 
Prosthesis", Clinical Orthopaedics and Related Research, 
1979, pp. 277-296. 
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4. R. M. Rose, M. J, Nuobuum, I, Paul, A, Crugnola, and 
E . h. Radin, "Wear Moehantnmo Coir Ultrahtgh Molecular 
Weight*, polyethylene in the Total Hip Prosthesis", J. of 
Appl. Polymer Bci., Vol. 2d, 1979, pp, 777-789, 

5. h, E. Lanyon, I, L. Paul, C. T. Rubin, fi. L, Thrasher, 

R, DeLaura , R. M, Rone and E* L, Radin, "In-Vivo Strain 
Measurements from Bone and Pr on then in Following Total 
Hip Replacement," J. = of Bone and Joint Burg., 63 -A , July, 
1981, pp. 989-1001. 

6. R. M, Rone, M, ivies, I. Paul, A. Crugnola and E. Bills, 
"The True Wear Rate of Ultrahigh Molecular Weight 
Polyethylene in the Total Knee Prosthesis," submitted to 
J. Bone and Joint Sure., 1981. 
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7. R. M. Rone, W. R. Cimino, B. Ellis and A. M, Crugnola, 

Exploratory Investigations on the Structure Dependence 
of Wear Resistance of Polyethylene," J.^of Wear, 1981, 
in press. 

8. R. M. Rose, "Wear Mechanisms and Wear Debris," Bulletin 
of the Hospital for Joint Diseases, 1980. 

9. E. L, Radin, I. L, Paul, R. M. Rose, "Osteoarthrosis as 
a Final Common Pathway," The Aotiqpathogenesis of 
Osteoarthrosis, G, Nuki, Pitman Publishing Co., Ltd., 
Kent, England, I960, pp. 84-89. 

10. R. M. Rose, A. M, Crugnola, W. R. Cimino and M. D. Ries, 

"The In Vivo Performance of polyethylene Components of 
Total Joint Replacements," Implant Retrieval: Material 

and Biological Analysis, NBS Special Publication #601, 
U.S. Department of Commerce, 1981, pp, 3-28. 


jExplqratqry Research in Nondestructive Testing 

This project is aimed at exploring the properties of defects 
in polymers and polymer-matrix composites as dielectric 
singularities and at identifying convenient methods of 
detecting defects in Ti weldments, based on similar 
approaches . 
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Sponoor: Office of Naval Research 

Fac ulty! U» M. Rose 
Sta||s I * M. Puffer 
Graduate Stude n fe s II. La nd i o 
Undergraduate .Student: J. Parse 

4, DESIGN TECHNOLOGY OF ADVANCED GRAPHITE/E POXY COMPOSITES 

Design, fabrication, and testing of graph! he/epoxy laminates 
and simple structural components. Evaluation of mechanical 
properties Including strength, fracture, fatigue, stiffness, 
vibrations, and buckling. Failure mechanisms and damage 
tolerance of components. 

Sponsor: Air Force Materials Laboratory, APWAL 

Facu lty: J, W, Mar, J, Dugundjl, T, H, H. Plan, E* Crawley 

Graduate students: 0. Bauchau, M. Graves, S. Hollowell, 

D. Jensen, P. Lagace, H. McManus, J. Rogers, M. Wong 
Undergraduate Students: About 15 part-time students 

Theses : 

1. J. C. Rogers, "An Investigation of the Damage Tolerance 
Characteristics c Graphite/Epoxy Pressure Vessels," 
M.I.T., Department of Aeronautics and Astronautics, 

S*M., September 1981. 

2. 0. A. Bauchau, "Desiqn, Manufacturing and Testing of 
High Speed Rotating Graphite/Epoxy Shafts," M.I.T., 
Department of Aeronautics and Astronautics, Sc.D., June 
1981. 

3. S. J. Hollowell, "Aeroelastic Flutter and Divergence of 
Graphite/Epoxy Cantilevered Plates with Ben Hng-Tor sion 
Stiffness Coupling," M.I.T., Department of Aeronautics 
and Astronautics, S.M., January 1981. 

4. D. W. Jensen, "NaturalVibration of Cantilever Graphite/ 
Epoxy Plates with Bonding-Torsion Coupling," M.I.T., 
Department of; Aeronautics and Astronautics, S.M., 
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August, 1981. 

5. 11. L. N. McManus, "Failure of Modes in a Family of 

Graphite/Epoxy Laminates", Department of 

Aeronautics and Astronautics, S.M,,June 1981. 

Publications : 

1. 0. A. Bauchau, "Experimental Measurement of Elastic 
Shear Modulus of Graph! te/Epoxy Tubes," J. of Composite 
Materials , 15, March 1981, pp» 151-156. 

2. J, W. Mar, M. J. Graves and D. P. Maass, "Effect of 
CompressionCompression Fatigue on Balanced 
Graphite/Epoxy Laminates with Holes," J. of Aircraft , 
18, 9, September 1981, pp. 744-747. 

3. P. Bar-Yoseph, T. H. H. Plan, "Calculation of 
Interlaminar Stress Concentration in Composite 
Laminates," J, Composite Materials, 15, 1981, pp. 
225-239. 


5. OTHER RELATED PROJECTS 


1 . 


Computational and Experimental Studies of 
Viscoelastic Flows 
R. A. Brown 

See — Mathematical and Physical Modelling of 
Materials Processing, page 238. 
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C. CERAMICS PROCESSING RESEARCH 


Faculty; H. K. Bowen 

Staff ; 

J. S. Haggerty 

R. L. Coble 


R. L. Pober 


R. M. Cannon 

A. Bleier 

B. J, Wuensch 

SUMMARY 

A new emphasis on processing of ceramic materials was begun 
at M.I.T. four years ago with the establishment of the 
Ceramics Processing Research Laboratory. The recognized 
purpose of this effort was to develop a processing science 
base for the ceramics industry. Because of unique 
institutional problems (at universities, national 
laboratories and corporations), there has not been a strong 
base in ceramics processing science as there has been in 
physical ceramics. Thus, part of the effort at MIT is to 
form a strong and interactive coalition between MIT, the 
government agencies, and the ceramics producers and users 
industries. Sixteen companies have joined the consortium, 
over 50 companies have made visits to the campus during the 
past year to discuss research results, and several companies 
have sent staff members for "internships" in the laboratory. 

The research on processing of ceramics is led by 5 faculty 
members and 1 senior research scientist and involves 4 post- 
doctoral scientists, 24 graduate students, and 8 under- 
graduates. Professor Coble continues his innovative research 
on experiments and models for sintering and microstructure 
evolution; Professor Cannon has developed models and data 
for understanding the kinetics of grain boundary motion and 
of pore coarsening? Professor Wuensch has developed new 
materials for fuel cells and batteries; Dr. Haggerty has 



developed now experimental techniques us i nq lasers for 
form inq fine, monos i zed powders and i or iorminq stable and 
metaatable a true t urea from liquids, and solar materials 
research in conjunction with the MIT Haora y Laboratory} 
Professor lUoier, who joined the faculty last year, has boqun 
unique research on colloidal science and ceramics proeossinq; 
and Professor Bowen has developed research programs on the 
presint erinq science necessary for controlled qroon micro- 
structures. Work is conducted in laboratories in Buildinq 
13 and in the Ceramics Process inq Research Laboratory located 
in Buildinq 12, The laboratory and tuoi l il ies manaqer ol the 
Ceramics Process inq Research Laboratory is Dr, Richard L. 
Pober . 


1 . LASP.R PROCRBB INU 

A tacUlty, consist Inq ot *1 hiqh power CO, lasers and 

peripheral equipment, has been established to eonduet 

directed enerqy materials process tnq research. The 

equipment is be inq used to conduct processinu research in 

four distinct areas with non-mot a 11 to materials# One custom 

deslqned, 1.300 watt laser Is inter, rated into a unique 

crystal qrowth IncilUy, it is be inq used to produce hiqh 

purity, hiqh molt inq point crystals ol oxide, carbides, and 

borides in luqhly controlled atmospheres. The same 

apparatus will be used to splat quench various classes and 

metastable materials which have been equilibrated in 

previously inaccessible atmospheres. Another system has 

been developed to synthesize small diameter, uniform 

spherical powders ol Hi, Hid, and These powders are 

used to supply the powder process inq research. The third 

apparatus is bo inq used to invest iqate a novel laser drive, 

photo- induced chemical, vapor deposition process, A tunable 

00, laser is beinq installed to conduct laser induced 
#* 

chemistry research . 
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dponnorn: 01 i' too oi Naval Ror.oaroh, Polemic Advanced 

Ror.oaroh Proieetn Aqonoy, Standard Oil ot 
Indiana, M IT Cabot Fund, Nat ional Science 
Pound a l ion , National Aeronaut ton and Space 
Adminintrat ton, IM Corporal ion 

Pacnltv : d, d. Uaqqorty 

d tat l ; d. 0. Pantorth, d. FUnt. 

Oraituato student n: R. Mavra, 11. dawhill, L. dohiolor 

Undorijraduat o Student ns K. Rloonr., M, C.abriot, C. Kerwin, 

M, dehmair, H. Sheldon f J . Hillman, 

M . Kn it! in , 0 . do o kh a r t , P . T i o 1 boh 1 , 

'Phonon : 

1, Premia Kinq, "Cryntal l i/.at ton Kt no lion in Amorphous? 
d i 1 i oon Th v n F i l mn , " Ma y , 1 l Hi l , 

2, a. Hillman, "Rlocl route dlnvrnct orient ton oi Lanor CVP 
Pi linn. " 

PuAUieat ions? : 

\. w, u . Oannon, d.C. Pantorth, d. d. Haqqerty, R, A. Marra, 

" dinterable Ceramic Powdorn t rom Lanor driven Reactions, 
Part I? Prooonn Peneript ion and Mode l inn." Submitted to 
a. Am, Ceram. doo. 

2. W. R. Cannon, d.C. Pan forth, a. d. Haqqerty, R. A, Marra 
"diniorablo Coramio Powdorn Irom Lanor Privon Reactions, 
Part It? Powder Charnel orinl. ion and Prooonn Variables.," 
Submit tod to H, Am. Co ram. doo. 

2, d. C. Pantorth, a.d. Hnqqorty, "Synfcheisn oi Ceramic 

Powdorn by Lanor Privon Reactions," To bo published in 
Co ram i o Pmj inoorinq and Science procoodinun. 

»l . d. Miauta , W. R. Cannon, A. Ploior, and a.d. Haqqerty , 
Piapornion and Cantinq oi Silicon Powder Without 
no t loocul ant n , " nuhmi tied to a, A m. Ce ra m, doc. . 

A. a. d. Hnqqorty, and VJ. R. Cannon, "dinterable Powdorn 
irom Lanor Privon React ionn." chapter in Laser Induced 
Chemical Pyoco sno n t Hd . , J.l, Sto infield , Plenum Pub. 

Co., New York, NY (1081). 

b. w. R. Cannon, d.C, Pantorth, a. 11. mint, a. d, Haqqerty, 
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anti R. A. Marta, "Synfchcsin oJ: Ceramic Powders 1’tom 
Laser Heated Gas Phase Reactants, " Proceedings Society 
Cor Photo-Optical and Instrumentation Engineers, Vol . 


198, Laser Applications in Ma ter lai n Process ing , Sa n 
Diego, CA 1979. 


2 • p RO C. E SSING OF C E RAM I C _ POWDE RS 

The most critical problem, limiting the extensive use of 
ceramic materials, is that high technology ceramics cannot 
be reliably and reproducibly manufactured. The focus of 
our research program is on the physics and chemistry of 
controlled particulate formation (size, size distribution, 
dispersion, phase and chemical composition, and shape) and 
of controlled packing of particles (no particle density 
gradients). The ceramic systems being studied are; AI 2 O 3 , 
BaTiOj , SrTiOj, DaPe 6 O ig , ZnO, Si, SiC, Si 3 N 4 , Si0 2 , T.iO^, 
TiB 2 and ZrQ 2> Powders are made by homogeneous gas phase 
nucleation (laser driven) or by controlled growth in aqueous 
and nonaqueous solutions. The colloidal dispersion 
parameters are being studied and models and techniques are 
being developed for controlled coagulation of the powders 
into shapes which are subsequently densified. tow firing 
temperatures and shorter sintering times are necessary in 
addition to better control of the microstructure and thus, 
the physical properties. 

The technological goals of this research are directed 
towards the use of ceramics as multilayer Si-chip carriers, 
capacitors, turbine and diesel engine components, permanent 
magnets, and varistors. 

Sponsors : AVX, Corning, Solid State Dielectrics, Union 

Carbide, Vitramon, Ford, GM-Delco, Hitachi, 
Tektronix, IBM, EXXON, Department of Energy, 

ILZRQ . 
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Faculty : A. Bleior, H, K. Bowen and J.S, Hagyerty 

Staff : R. L. Pober , J. Blendell, S. Mizuta, Y. Suyama , 

H. Itoh, N. Levoy, 

Graduate Students : E. Barringer, C. Huynh, M. Parish, 

M. Strauss, M. Green, R. Rir.an, N. Jubb, 
E. Tormey, R. Garcia, T, Kramer. 

Undergraduate Students ; N* Levoy 

Theses : 

1. N. Levoy, "Dispersion, Classification and Sintering of 
AI 2 O 3 ," B.S., June, 1981. 

Publications : 

1. E, S. Tormey, et al., "Absorption of Dispersants from 
Nongqueous Solutions," Surface and Interfaces in Ceramic 
and Ceramic-Metal Systems, Eds., J. Pask and A. Evans, 
Plenum Press, NY, 1981, pp. 121-136. 

2. E. Barringer, "Disorder-Order Transition in Monodisperse 
T.itania Sols," M.I.T. Ceramics Processing Research 
Laboratory, Report # 8 , December, 1980. 

3. MATERIALS FOR SOLAR ENERGY 


Materials processing research is being conducted in two 
solar energy areas. In the first, graded-.index of 
refraction anti-reflection (AR) coatings are being developed 
on silicate glass by preferentially removing one phase of a 
phase separated glass. When used on covers of flat plate 
collectors, these broad band AR coatings increase the 
extractable heat by 30-50%. The second area involves 
several materials processing programs for photovoltaic 
devices. Our analyses indicate that amorphous processes 
offer the best chance of achieving required cost objectives. 
Consequently, we have limited our research to this focus. 
Processes to improve carrier life times in chalcogenide 
glass semiconductors have been investigated to take 
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advantage o£ their low cost and adjustable band gap. 
Processes to control the nucleation and growth processes in 
amorphous silicon are being studied to make use of the low 
cost of amorphous deposition processes and the high 
efficiency of crystalline Si, The laser CVD deposition 
process is being developed because of low manufacturing cost 
and superior properties through better process control. 

Sponsors ; Department of Energy, SERI, MIT Cabot Fund, 

3M Corporation 

Faculty : J. S. Haggerty, D. Adler 
Staff : S. C. Danforth, T. Gattuso, D. Imeson 

Graduate Students : A. Iqbal, M . Meunier, F , Van G.ieson 

Undergraduate Students : J. Hillman, C. Kerwin, B. King, 

B. Sheldon, 0, Estell 

Theses : 

1. A. Iqbal, "Determination of Chemistry of Gr aded-Index , 
Antlreflective Films on Glass," S.M. , February 1981, 

2. B. King, "Crystallization Kinetics in Amorphous Thin 
Films," May, 1981. 

3. V. Tengzelius, "Development of Antiref lect.ive Films on 
N0 2 /0/Ca0/S.i02 Glass," January, 1982. 

4. F. VanGieson, "Kinetics of Crystallization for Selected 
Amorphous Silicon Films," January, 1982, 

Publications ; 

1. A. Iqbal, J. S, Haggerty, and S. C. Danforth, "Surface 
Chemistry of Anti-Reflective Films on Borosilicate 
Glasses," submitted to the J . Am. Ceram. Soc. 

2. J. S. Haggerty, B. Sheldon, and A. G. Emsl.ie, "Exact 
Computation of the Reflectance of a Surface Layer of 
Arbitrary Refractive Index Profile and an Approximate 
Solution of the Inverse Problem," submitted to the 
J, Opt. Soc, Am, 

3. J. E. Ritter, Jr., K. Jakus, K. Buckman, G. Young, nad 
J. S. Haggerty, "Strength and Fatigue Behavior of a 
Borosilicate Glass with an Antiref lective Surface," 



submitted to Glass Technology . 

4. S. C, Danforth, F. VanGieson, I, Kohatsu, and J. S. 
Haggerty, "Laser Induced Controlled Nucleation and 
Growth Process for Large Grained Polycrystalline 
Silicon," American Institute of Physics, Conference 
Proceedings, Laser and Electron-Beam Solid Interactions 
and Materials Processing , Ed, Sigmon, Gibbons, Hess, 
Boston MA. (1980) . 

5. W. R. Cannon, S. C. Danforth, J. Flint, J. S, Haggerty, 
and R. A. Marra, "Synthesis of Ceramic Powders from 
Laser Heated Gas Phase Reactants," Proceedings Society 
of Photo-Optical and Instrumentation Engineers, Vol. 
198, Laser Applications in Materials Processing , San 
Diego, CA (1979). 

4 . PROCESSING OF OXIDE POWDERS 


A program to examine the application of current 
colloid-chemical models for single oxides to the processing 
of technical ceramic powders has been initiated. Its 
purpose is to extend the models to heterogeneous systems 
containing multimetallic oxides. Objectives are to prepare 
model, single and multimetallic oxides, to characterize the 
materials prepared using crystallographic, chemical, 
physical, and surface-chemical procedures, and to test 
recently proposed models of the electrical double layer. 

The approach consists of evaluating surface reactions and 
properties, determining intrinsic surface ionization 
constants, comparing surface and bulk solution chemical 
equilibria, and comparing computed and experimentally 
determined specific conductivities for dispersions relevant 
to high performance ceramics. Results will relate surface 
ionization reactions and electrical double layer properties 
to fundamental interactions among solutes, solvents 
(principally water), and solid x>articles. 
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The overall technical (goal of Ulna research is to develop 
predictive capabilities lor the processing of real and 
"model" ceramic oxides . 


Sponsor : Department of Energy 

Faculty: A. Bleier 

Graduate Student : W. C. Ilasz 

Publ ica tions : 

1. S. M.ii'uta , W. R. Cannon, A. Bleier, and J. S. Haggerty, 
"Wetting and Dispersion of silicon Powder Without 
Deflocculants," submitted for publication, Am, Ceram. 
Soc . 

2. A. Bleier, "Role of van der Waals Forces in Deter- 
mining the Wetting and Dispersion Properties of Silicon 
Powder," submitted for publication. 

3. A. Bleier and R. M. Cannon, "Synthesis and 
Characterization of Uniform Zirconium Dioxide," to be 
submitted tor publication. 


5. SINTERING AND MICRQSTRUCTURE EVOLUTION 

Most of the important technological developments for 
processing ceramic materials have resulted from theoretical 
models of sintering and grain growth coupled with an enor- 
mous number of experimental measurements. The efforts at 
MIT over the past 25 years have been at the frefront of the 
development of these concepts, and most recently has suc- 
ceeded in modelling the very complex simultaneous processes 
which occur during the firing of ceramic greenware. 

Because of the extreme sensitivity of ceramic materials to 
.impurities and additives, and because of the difficulties in 
characterizing fine powders, these theoretical analyses pro- 
vide guidance to the selection of experimental parameters 
(T, P, composition, particle size and distr ibu-tion) for a 
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critical set of! measurements to tost the basic theories as 
well as providing guidance to the ceramic processors. These 
analyses are also demonstrating the importance of the local 
distribution of porosity, impurities or dopants, and grain 
or particle size in determin-ing the densif ieation process 
and microstructure evolution. 

Numerous model materials have been chosen which have widely 
varying difference in diffusion coefficients, vapor 
pressures, surface and grain boundary energies, etc. These 
include: Al 2 0 3 , Zr0 2 , Au , ZnO, Si, Sic, Si 3 Ni 4 , MgO, LiF, 

B9Ti0 3 , and Th02-Zr0 2 . In each case, fundamental studies 
are performed which relate the basic thermochemical forces 
and transport coefficients to the observed evolution of the 
sintered microstructure. 

Some of the particularly exciting recent accomplishments 
are : 

1. analysis of the grain boundary grooving technique as it 
applies to surface diffusion studies, in particular, the 
effects of faceting and simultaneous lattice and surface 
contributions ; 

2. modelling of sintering kinetics based on virtually mono- 
dispersed, 1 m gold powder; 

3. investigation of the complex sintering hehavior of 
ZnO-doped ceramics to determine the nature and extent of 
the sintering mechanism; 

4. partial demonstration of the notion that the Herring 
scaling laws also apply to nonspherical powder as well; 
and 

5. analysis of the transport processes which occur in the 
Ni-W system using isothermal creep tests. 

Sponsor : Department of Energy 

Faculty : R. L. Coble, R. M. Cannon 

Graduate Students : J. Dynys, W. Hong, E. Rothman, H. Song, 

B. Zelinsk.i 



Theses : 


1. VI. Coblenz, "Physics of sintering of silicon," Sc.D,, 
MIT, Department of Materials Science and Engineering, 
June 1981. 

2. A„ Roshko, "Grain Coarsening in Doped BaTiOg," S.B., 

MIT, Department of Materials Science and Engineering, 
June 1981. 

Publications : 

1, J, M, Dynys , R. L. Coble, W. S, Coblenz, and R. M. 
Cannon, "Mechanisms of Atom Transport during Initial 
Stage Sintering of Al 2 0g," Sintering Processes , Edited 
by G. C. Kuczynskl, Plenum Publishing Corporation, 1980. 

2. W. S. Coblenz, J. M. Dynys, R. M. Cannon, and R. L. 
Coble, "Initial Stage Solid State Sintering Models. A 
Critical Analysis and Assessment," Sintering Processes , 
Edited by C. G. Kuczynski, Plenum Publishing 
Corporation, 1980. 


6. SYNTHESIS AND PROPERTIES OP FAST- ION CONDUCTORS 


Fast-ion conductors, materials which display ionic 
electrical conductivities up to six or eight orders of 
magnitude larger than normal ionic compounds (i.e., up to 5 
reciprocal ohm-cm), find important application in fuel 
cells and battery systems. We are conducting a broad study 
of such materials, which includes synthesis and power 
processing of potential new conductors, fabrication of 
sintered specimens for measurement of electrical properties, 
growth of single crystals, and precise neutron and x-ray 
diffraction measurements to provide insight into the 
mechanisms of the fast-ion conduction process. The 
materials which are currently the subject of exploratory 
synthesis and processing are alkali metal conductors of 
relevance to battery systems. Potential new conductors 
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include several alkali metal silicates and titanates, and 
nitrogen-based ceramics. In addition, novel solid-solution 
series are being prepared with systems which have previously 
bean demonstrated to display Cast-ion conduction. Among 
these are NASICON {Na 3 Zr 2 SiP0 12 ) and LISICON 
(Li 3 t 5 ^ 0 . 25 ^ 04 ) related systems. Phases Cor which single 
crystals have been synthesized and Cor which diffraction 
studios have begun are primarily the prototype 
cation-disordered Ag and Cu conductors. These are ideal 
systems for critical examination of models for the 
distribution and thermal motion of the mobile cations as 
they constitute a large fraction of the scattering density 
in these materials. These phases include Cu 2 S, CuAgS, Ag 3 SI 
and Ag 2 Se. 

A variety of processing techniques are employed, depending 
upon the system of interest. Solid state reaction of 
component oxides or salts are frequently employed, as are 
hydrothermal reaction and preparation of gels. The growth 
of single-crystal materials employs (again depending on the 
system) hydrothermal techniques, flux growth or vapor 
transport. 

Sponsors ; Department of Energy, Lawrence Berkeley Laboratory 

Faculty : B.J. Wuensch 

Staff : I. Kohatsu, A.-K Ekholm 

Graduate Students : C.L. Skarda 

Publications ; 

1. R.J. Cava, F. Reidinger and B.J. Wuensch, "Conductivity 

Mechanisms in the Superionic Phases of Agl and Ag2S as 
Determined by Neutron Diffraction," in Fast Ion 
Transport in Solids: Electrodes and Electrolytes , P. 

Vashishta, J.N. Mundy and G.K. Shenoy, Eds. Elsevier 
North Holland, 1979, pp. 217-220. 

2. R.J. Cava, F. Reidinger and B.J, Wuensch, "Single- 
Crystal Neutron Diffraction Study of the Fast-Ion 
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Conductor -Aq 2 S Between 18G end 325°C, J. Solid State 
Chom,_, Vol . 31, 1980, pp, 69-80, 
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D. MATERIALS SYSTEMS ANALYSIS 


Faculty: M.B. Revet 

J.P. Clark 
A.M. Church 
J.F. Elliott 

M. C. Flemings 

N. J. Grant 
T.B. King 
J. Szekely 
M.D. Zimmerman 


SUMMARY 


Staff: G.B . Kenney 

S.M. Mathur 


Recognizing the need to better understand the societal, 
economic, and policy tradeoffs associated with materials 
processing and utilization, the Materials Systems Analysis 
Group was established within the Department of Materials 
Science and Engineering in 1975. The principal purpose of 
the materials systems effort is to provide the materials 
engineer with the systems analysis required to formulate 
sound materials processing, utilization, and resource 
development policies and strategies for private industry and 
government. 



The interdisciplinary research on materials systems analysis 
involves the collective efforts of 9 faculty members, 2 
research staff members, n graduate students, and 1 under- 
graduate from the Departments of Materials Science and 
Engineering, Sloan School of Management, and Centers for 
Materials Processing and Policy Alternatives. Prof. Clark 
continues to expand his innovative materials system 
simulation and modeling research program which includes 
assessments of materials substitution dynamics, public 
policy implications, and materials process economics. This 
effort includes several collaborative programs with 
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materials engineers, economists, and policy analysts. Prof. 
Bever coordinates an ongoing program concerned with the 
technical and socioeconomic issues of materials recycling, 
input-output analysis, and technological change and 
productivity. The major thrust areas in materials systems 
research are outlined below. 


1. MODELING OF MATERIALS SUPPLY, DEMAND, AND PRICES 


Imbalances in the supply and demand of basic materials has 
periodically led to short term shortages and gluts which 
have sometimes resulted in substantial fluctuations in 
prices. Although an efficient free market system can be 
expected to balance supply and demand in the long run, there 
is a need to anticipate the conditions that lead to these 
short term fluctuations to facilitate more efficient 
planning at the level of the nation, industry, and firm. We 
have been developing detailed engineering and economic 
models of materials industries for these purposes. The 
models are usually built from the supply side from an 
engineering process analysis basis and then transformed to 
cost and supply functions, which are estimations of the 
relationship between output and price. The demand models 
are usually estimated by statistical analyses as a first 
approximation and then modified with engineering judgment. 
Thus far we have or are in the process of applying 
variations of this approach to the following industries: 
iron and steel (carbon steel, alloy steel, stainless steel 
and foundry products), aluminum, copper, manganese nodules, 
magnesium, and cobalt. 

Sponsors; U.S. Office of Surface Mining, U.S, Bureau of 

Mines, National Science Foundation 
Faculty : M. B. Bever, J. P. Clark, J. F. Elliott, M. C. 

Flemings, N. J. Grant, T. B. King, J. Szekely 
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Staff : G, B. Kenney, S. M. Mathur 

Graduate Students ; P. Baverstam, M. Cummings, P. Field, 

P. Foley, S. Kometani, J, Newman, 

D. Richards, D. Wadekar, R. Upaa 

Publications : 

1. F. E. Katrak, T. B. King and J. P. Clark, "Analysis of 
Supply of and Demand for Stainless Steel in the United 
States," Materials and Society , 4_, No. 4, p. 427, 19 80 . 

2. F, E. Katrak and J. P. Clark, "An Analysis of the 

Determinants of Stainless Steel Pricing and Imports in 
the United States: Part I - Pricing," Materials and 

Society , 4_, No. 4, p. 437, 1980. 

3. J. P. Clark and F. E. Katrak, "An Analysis of the 

Determinants of Stainless Steel Pricing and Imports in 
the United States: Part II - Imports," Materials and 

Society , 4_, No. 4, p. 447, 1980. 

4. J. P. Clark, N. J. Grant, and T.B. King, "The Potential 
for Utilization of Manganese Derived from Deepsea 
Nodules," Natural Resources Forum , 5_, p. 249, 1981 . 

5. J . P N Clark and G. B. Kenney, "The Dynamics of 

Intermaterial Competition in the Automotive Industry, 
Part I: A Framework for Analysing the Dynamics of 

Intermaterial Competition," Materials and Society , 5_, 

No. 2, 1981. 

6. G. B. Kenney and J. P. Clark, "The Dynamics of 

Intermaterial Competition in the Automotive Industry, 
Part II: A Case Study of the Demand for Magnesium," 

Materials and Society , 5_, No. 2, 1981. 

7. P 0 Foley and J. Clark, "U.S. Copper Supply - An 
Economic/Engineering Analysis of Cost-Supply 
Relationships," Resources Policy , 7_No. 3, p. 171, 1981 

8. B. J. Reddy and J. P. Clark, "Effects of Deepsea Mining 
on International Markets for Copper, Nickel, Cobalt, and 
Manganese," in Deepsea Mining , J. T. Kildow, ed . , M.I.T. 
Press, 1980. 
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2 . PUBLIC POLICY 

Public policy as enunciated by federal, state, and local 
regulations shapes corporate policies which determine 
processing strategies and resource development, or the lack 
thereof. Consequently, public policy is an important 
endogenous variable and/or tool to the materials systems 
analyst. One specific study evaluates the effect of public 
regulations on the U.S. copper industry by developing a 
methodology for assessing the costs and benefits of 
regulations . 

Sponsor : National Science foundation 

Faculty ; J. P, Clark, M. B. Zimmerman, A. M. Church 
(University of New Mexico) 

Graduate Student* P. Foley 
Publication ; 

1. A. M. Church, P. T. Foley, II. B. Zimmerman, and 3 . P. 
Clark, "The Effects of Taxation of the U.S. Coal and 
Copper Industries," Materials and Society , 5_, No. 1, 
p. 81, 1981. 

3. RECYCLING 


Research on the technological and socioeconomic aspects of 
materials recycling was continued. The research was 
concerned with resource recovery from waste, the recovery of 
metal scrap, the recycling of automotive materials, and the 
recycling of chromium. 

In the area of resource recovery, an undergraduate student 
worked in the laboratories of the Raytheon Company, Research 
Division under the supervision of Dr. Ernst F. Schloemann as 
part of the Undergraduate Research Opportunities Program. 
This work contributed to the redesign and testing of an eddy 


current separator with the objective of Increasing its 
capacity. Work planned for the next phase will involve the 
magnetic separation of light and heavy scrap recovered from 
municipal solid waste and the separation of aluminum from 
shredded polyester bottles. 

The analysis of technical, economic, and institutional 
aspects of the recovery of resources from waste was 
continued. This analysis was carried out in support of a 
study by a panel on the "Recovery of Energy and Materials 
from Solid Waste" of the Building Research Advisory Board of 
the National Research Council, 

As a continuation of earlier research on the disposal of 
discarded automobiles, detailed planning was carried out for 
a study on the "Impact of Materials Substitution on Energy 
Conservation, Environmental Goals and Resources Recycling; 

A Case Study of Recyclability in the Automotive Industry." 
This study will be conducted as part of the research program 
of the American Society of Mechanical Engineers with funding 
from the National Science Foundation, 

As an extension of the work on recyling, the strategies of 
conservation of materials efficient utilization, substitu- 
tion, recycling and research and development were reviewed. 

An analysis made previously for the Congressional Research 
Service was updated for publication. 

Sponsors : Raytheon Company; National Science Foundation 

(through American Society of Mechanical Engineers) 
Facul ty ; M. B. Dover 
Undergraduate ; D. G. .Facinelli 
Pu 1 ications ; 

1. M. B. Bever, "The Potential for Reduction of Import 
Dependency; Efficient Materials Utilization, 
Substitution, Recyling, and Research and Development," 
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Congressional Seminar, in emerging Issues in Science and 
Technology, Congressional Research Service, U. S. 
Government Printing Office, Washington, 1980, pp. 

88 - 100 . 

2. D. G. Facinelli, "The Development of Non-Ferrous Metal 
Separators Using Permanent Magnets and Their Application 
to Municipal Solid Waste," S.B. Thesis. 

3. E. Schloemann and D. Facinelli, "High-Capacity 
Nonferrous Metal Separator Using Permanent Magnets," 
Proc. 2nd International Symposium on Materials and 
Energy from Refuse, Antwerp, Belgium, October 1981 (to 
be published). 

4. APPLICATIONS OF INPUT-OUTPUT ANALYSIS 

The production and consumption of twenty-six basic minerals 
has been analyzed within the framework of input-output 
models at the institute of Economic Analysis, New York 
University, under the direction of Professor Wassily 
Leontief. In research carried out in collaboration with Ms. 
Sylvia Nasar, the recycling of six metals was analyzed for 
three scenarios. This analysis will be continued with 
revised data and expanded to include environmental and 
energy aspects. 

A study of the recycling of chromium by input-output 
analysis will be completed. It will be complemented by an 
investigation of structural aspects of the recovery of 
secondary chromium and its use by industry. 

As a continuation of the work mentioned in the foregoing, 
planning was completed for a project on "Evaluating the 
Impact of Prospective Changes in Materials Use: An 

Input-Output Approach." This project will be conducted at 
New York University and will be sponsored by the Bureau of 
Mines . 
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Sponsor ; Bureau of Mines (at New York University) 

Faculty ; M. B. Bever 

Publication ; 

1. S, Nasar-0 ' Br ien and M. B. Bever, "Applications 

of Input-Output Analysis to Investigations of Metals 
Recycling," Proceedings, Council of Economics, AIME , 
1981, pp. 27-40. 

5 . MATERIALS, TECHNOLOGICAL CHANGE AND PRODUCTIVITY 


The relations of materials and productivity are of critical 
importance for economic performance and long-term resource 
availability. These relations are greatly affected by 
technological change. An investigation of materials, 
technological change, and productivity was undertaken in 
collaboration with Ms. Sylvia Nasar. Results were 
presented at the 3rd Conference on Productivity at Rutgers 
University, April 23, 1981. The research is continuing. 

Faculty ; M. B. Bever 
Publication ; 

1. M. B. Bever and S. Nasar, "Materials, Technological 
Change and Productivity," (accepted for publication). 
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E. WELDING RESEARCH 

Faculty: T. W. Eagar 

D. E, Ilardt 
J ■ H • Lang 
K. Masubuchi 
H. M. Paynter 

J. Szekely 

K, Terai (Visiting) 
W. C. Unkel 
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Staff: c. Allemand 

A. Block-Bolten 

A. Xrnakita 
(Visiting) 

G. Oreper 

B. Russell 
A. Zona 


SUMMARY 


Although continuously present Cor many decades, the welding 
effort at MIT has expanded rapidly in the past five years. 
Presently, this constitutes the largest academic program in 
welding in the United States. It is led primarily by Pro- 
fessors Masubuchi and Eagar who study Welding Fabrication 
and Welding Processes, respectively. A number of other 
faculty are actively involved in many ways, creating an 
interdisciplinary effort centered in the Departments of 
Ocean Engineering, Mechanical Engineering, Materials Science 
and Engineering, Electrical Engineering, and Materials 
Processing Center. 

Personnel include eight faculty, four research scientists 
and engineers, two technical staff members, twenty-one 
graduate students, and three undergraduate students. A 
selection of the ongoing research is presented below. 


1. WELDING FABRICATION 


The research effort during the period from October 1, 1980 
through September 30, 1981 centered around the following 
projects : 
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- Residual Stresses and Distortion in Structural 
Weldments in High-Strength Steels 

- Improvement of Reliability of Welding by In-Process 
Sensing and control (Development of Smart Welding 
Machines Cor Girth Welding oC Pipes) 

Development oC Fully Automated and Integrated 

( "Instnmatlc" ) Welding Systems for Marine 
Appl ications 

- Further Advancement of Welding Technology 

- Research on Metal Working by High- Power Laser, 

Further details of the research activities are described in 
the following pages. The Kawasaki Heavy Industries Research 
Fund also provided financial assistance in conducting 
research on woldinq fabrication. 

The final report of a proposal entitled "Development of 
Joining and Cutting Techniques for Deep-Sea Applications" 
was issued from the M.I.T. Sea Grant Office in June 1981 as 
MITSG 81-2. The objective of the research program from July 
1976 through Juno 1980 was to (1) generate basic information 
pertinent to joining and cutting techniques for deepsea 
applications and (2) to develop some prototype tools 
suitable for underwater joining and cutting tor deep-sea 
appl ications . 


Residual Stresses and Distortion in Structural Weldments in 
High Strength Steels 

The objective of the original three-year research program, 
which started on December 1, 1977, was to analytically and 
experimentally study residual stresses and distortion in 
structural weldments in high-strength steels, especially 
HY-130 steels. The program included the following tasks: 
Task 1: Research on thick butt welds 
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Task 2: Research on girth v/e Ids of cylindrical 

shells 

The program has been extended Cor another year to cover the 
following additional tasks; 

Task 3; Development of information on weldments in 
bronze 

Task 4; Research on thermal stress relieving of 
weldments 

Task 5; Development of improved computer programs, 

Most of the efforts under Tasks 1 and 2 were completed 
before November 1980 and the results wore included in the 
first and second technical progress reports issued in 
November 1979 and November 1980, respectively? The 
remainder of Tasks 1 and 2, primarily the completion of two 
Ph.D. theses, were completed by May 1981. Task 3 was also 
completed by June 1981. These results and results of the 
work under Tasks 4 and 5, which will be completed by 
November 1981, will be included in the third technical 
progress report which v/e plan to issue in November 1981. 
However, we may not be able to complete all the intended 
efforts by November 30, 1981. Those efforts which are 
likely to be unfinished include (1) experiments under Task 4 
and (2) manuals of some two-dimensional finite element 
programs for analyzing transient thermal stresses and 
residual stresses in weldments. We are asking the sponsor 
for an extension of six months; therefore, the entire 
research program will be completed by May 30, 1982. We plan 
to issue a final report at the conclusion of the program. 

Sponsor : Office of Naval Research 

Faculty : K. Masubuchi 

Staff : A. Imak.ita (Visiting), A. J. Zona 

Graduate Students: J. Agapakis, E. Goncalves, R. S. McCord, 


V. J. Papazoglou 



1, V. J, Papazoglou and K. Masubuchi, Second Technical 
Progress Report on "Study of Residual Stresses and 
Distortion in Structural Weldments in High-Strength 
Steels," Office of Naval Research, Contract No. 
N00015-75-0469, November 1980. 

Theses : 

1. E. Conceives, "Fracture Analysis of Welded Structures". 

2. R. S. McCord, "An Investigation of Strain, Distortion, 
and Heat Distribution During Welding of Nickel-Aluminum 
Bronze" , 

3. V. J. Papazoglou, "Analytical Techniques for Determining 
Temperatures, Thermal Strains, and Residual Stresses 
During Welding". 

Publ icatlons : 

1. V. J. Papazoglou and K. Masubuchi, "Analytical Methods 
for Determining Temperatures, Thermal Strains, and 
Residual Stresses Due to Welding," a paper presented at 
the 1981 Annual Meeting of the American Welding Society, 
Cleveland, Ohio, April 6-9, 1981. 

2. K. Masubuchi, "Welding Stresses," presented at the 28th 
Sagamore Army Materials Research Conference, Residual 
Stress and Stress Relaxation, Lake Placid, New York, 

July 13-17, 1981, 


Improvement of Reliability of Welding by In-Process Sensing 
and Control (Development of Smart Welding Machines for 
Girth Welding of Pipes) 

The overall objective of this three-year research program 
which started on June 15, 1979, is to improve the reliabi- 
lity of welding by developing "smart" welding machines. A 
smart welding machine is equipped with sensors, artificial 
intelligence, and actuators with the goal of reducing welding 
errors by one or two orders of magnitude. Although the 
concepts and techniques which will be developed in this 
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program will have more general applications to welding, thin 
research program is focused specifically on welding of 
pipes. 


Sp onsor : Department of Energy 

Fac ulty: D. E, Hardt, K, Masubuchi, H, M. Paynter, and 

W. C. Unkel 

Graduate Students : J. Convert!, Y, Dr or, S. Liang, J. Moore, 

M. Zacksenhouse 

Research Report : 

1. Third Progress Report on "Improvement of Reliability of 
Welding by In-Process Sensing and Control (Development 
of Smart Welding Machines for Girth Welding of Pipes),", 
to the Department of Energy under Contract No, 
DE-ACQ2-79ER10474 , A000 , June 1981. 

Theses : 

1, Y, Dror, "Calculations and Measurement of Weld Puddle 
Geometry in Gas Tungsten Arc Welding", 

2. J, Convert!, "Plasma- Jets in Arc Welding". 


Development of Fully Automated and Integrated ( "Instamatlc" ) 
Welding Systems for Marine Applications 

The objective of this two year research program which 
started in July 1980 is to develop fully automated and 
integrated welding systems, which may be called "instamatic" 
welding systems for various marine applications. These 
systems which will be developed through careful engineering 
will be able to perform certain prescribed welding jobs by a 
person with no welding skill. In many cases welding will be 
performed in a completely enclosed system so that no spark 
and very little fume will be generated from the welding 
systems , 

Sponsors : National Sea Grant Office, National Atmospheric 

and Oceanic Administration, Department of Commerce, 
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Faculty ; K. Masubuchi 

Staff; A. J. Zona 

Graduate Students ; H. L. Gustin and D. W. Schloerb 

Publications ; 

1. K. Masubuchi, "Underwater Factors Affecting Welding 
Metallurgy," a paper presented at the Seminar on 
Underwater Welding of Offshore Platforms and Pipelines 
by the American Welding Society , New Orleans, LA, 
November 5-6, 1980. 

2. K. Masubuchi, "Review of Underwater Welding Technology," 
Oceans '81 - Proceedings of a Conference held in Boston, 
MA, September 16-18, 1981, pp. 649-651. 


Further Advancement of Welding Technology 

This three-year research program, which started in October 
1980, covers the following three subjects; 

1 . Automation of welding processes 
2 , Computer applications to welding 

3. Improvement of reliability of welded structures. 

Sponsors ; A group of industrial sponsors including Hitachi, 
Ltd., Hitachi Shipbuilding and Engineering Co., 
Ishikawaj ima-Har .ima Heavy Industries, Kawasaki 
Heavy Industries, Kubota Ltd., Kobe Steel Works, 
Mitsubishi Heavy Industries, Mitsui Engineering 
and Shipbuilding Co., Nippon Kokan Kaisha, Nippon 
Steel Corp., Osaka Transformer Co., Sumitomo 
Heavy Industries, and Tokyo Sh.ibaura Electric Co. 
Faculty ; K. Masubuchi and K. Tera.i (Visiting) 

Staff ; A. Imakita (Visiting) and A. J. Zona 

Graduate Students ; J. Agapak.is, E. Goncalves, and V. J. 

Papazoglou 

Publications: 


1. K. Masubuchi, "New Directions in Welding Reserach and 
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Development," a special lecture given at the 
International Con lor once on Welding Research In the 
,1980 1 s , Osaka, Japan, October 27, 1980. 

2. A. Imakita, V. J. Papaaoglou , and K, Masubuahi, 
"Annotated Bibliography on Numerical Analysis of 
Stresses, Strains, and Other Effects Due to Welding," 
Document x-99t-81 presented at the 1981 Annual Assembly 
of Commission X of International Institute of Welding, 
Oporto, Portugal, September 8-11, 1981. 


Research on Metal Working by High-Power baser 


This research program includes the following tasks: 

Task 1: Study on application? of high-power laser for 

motai working. 

Task 2: Experiments on characteristics of high-power 

lasers , 


Sponsor : Japan Welding Engineering Society 

Faculty : K» Masubuchi and K. Torai (Visiting) 

Staff: A. IrnakUa (Visiting) 


2 . 


WELDING PROCESSES (T.W. EAGER) 


The major thrusts of the new program include: 

Fluxes 

Sensors for automation (see Materials Processing 
Research Base Abstract) 

Heat and fluid flow 

Mechanical properties of weldments. 

The majority of the research is concerned with arc welding 
although studies of resistance spot welding and laser 
welding are underway. These topics are described in greater- 
detail below. 
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Flux Development 

Flux shielded welding processes account for the largest 
quantity of welded products; yet our understanding of the 
flux chemistry and slag-metal reactions remains poor. 

Reactions involving manganese, silicon, chromium, and oxygen 
are being studied in fused and bonded submerged arc fluxes. 

A thermodynamic analysis of the slag-metal equilibrium has 
been developed and tested with remarkable success. 

Investigation of active fluxes and flux-cored electrodes is 
planned. Studies of halide fluxes for use with titanium 
have shown a number of advantages. Submerged arc welding, 
gas metal arc, and electroslag welding of titanium are all 
being studied. 

Sponsors : National Science Foundation, Office of Naval Research 

Faculty ; T. W. Eagar 

Staff : A. Block-Bol ten , B. Russell 

Graduate Students : S. K. Fan, U. Mitra, D. Ries 

Undergraduate Students : T. Lynch, R. Schoder 

Publ ications : 

1. C. S. Chai and T. W, Eagar, "Slag-Metal Equilibrium 
During Submerged Arc Welding," Met. Trans. , 12B, 1981, 
p. 539. 

2. C. S. Chat and T. W. Eagar, "Slag Metal Reactions in 
Binary CaF 2 -Metal Oxide Welding Fluxes," accepted by the 
Welding J. 

3. C. S. Chai and T. W. Eagar, "Prediction of Weld Metal 
Composition During Flux Shielded Welding," to be 
published in the Proceedings of the ASM Conference on 
Welding Consumables and Process Developments , Peoria, 


IL, September 1981. 
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Heat and Fluid Flow 


A number of investigators have provided heat flow models for 
welding involving a moving point source of heat. While this 
approximation works well away from the source, it produces 
erroneous predictions near the weld pool. A comprehensive 
model allowing for a more realistic distributed heat source 
is being developed. Both theory and experiment show that 
the size of the weld pool and the HAZ can vary by a factor 
of two or more for equivalent heat inputs. 

Studies of fluid flow in the welding arc plasma and the 
molten weld pool are underway. The bulk of this work is 
theoretical. Somer experimental studies of the effects of 
convection on the size and shape of the weld pool are 
underway . 

A comprehensive model of heat and fluid flow in electroslag 
welding has been completed and tested experimentally. This 
study indicates that there is little opportunity to vary the 
size of the heat affected zone by altering the welding 
process parameters. 

The mechanism of heat transfer during laser welding of alu- 
minum is being studied by measurements of reflectivity and 
plasma formation as functions of surface preparation and 
alloy content. Initial results indicate differences of :i 
factor of 20 in the initial heat transfer rate depending on 
the particular alloy and surface preparation. 

Sponsors : Department of Energy, Office of Naval Research 

Faculty ; T. W. Eagar, J. Szekely 

Staff : C. Allemand, A. Blcck-Bol ten , G. Oreper 

Graduate Students ; M. Lin, N. S. Tsai 

Undergraduate Student ; G. Dunn 

Theses ; 

1 . 


N. A. Dudziak, "Spectrograph.ic Measurement of Metal 
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Vapor Concentration in a Welding Arc,’ 1 S.B., June 1981. 
2. C. A. Huntington, "The Reflectivity of a C0 2 Laser Beam 
by Aluminum," S,B , June 1981. 


Mechanical Properties of Weldments 

Factors influencing the fracture toughness of HY80 weld 
metal are being studied through investigation of high heat 
input variations with commercial submerged arc fluxes, and 
the effect of nitrogen on strain aging of the weld metal. 

In the former, the effects of flux composition and welding 
process parameters on inclusion composition, size, and 
distribution as well as phase transformation behavior is of 
interest. In the latter project, the effect of weld metal 
composition on the state of nitrogen is being studied with 
concern for the effects of residual stresses and thermal 
cycles in multi-pass welds. Internal friction is being used 
to follow the combination of nitrogen in submerged arc and 
gas metal arc welds. 

The effects of weld repair on the fatigue properties of cast 
titanium have been studied. It is concluded that weld 
repair is not detrimental but that the final surface 
finishing technique should be chosen with care. In 
particular, chemical milling of the final part is inferior 
to either mechanical polishing or peening of the surface. 

Sponsors ; Office of Naval Research, General Dynamics 

Corporation, U.S. Army Materials and Mechanics 
Research Center 
Faculty ; T. W. Eagar 

Graduate Students ; 0. Boyd as, B. Wilson, G. Hunter 
Thesis ; 

1. G. B. Hunter, "Processing Effects on the High Cycle 
Fatigue Life of Weld Repaired Cast Ti-6A1-4V Parts," 

S ,M. , June 1981 . 
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3. EXPLORATORY RESEARCH ON NONDESTRUCTIVE TESTING 


This project has two goals: 

(a) the exploration of the dielectric singularities 
associated with defects in polymer-matrix composites 
in order to ascertain the degree of detectability 
and differentiation of the various defect types. 

(b) exploration of the electrical properties of titanium 
weldments to ascertain the feasibility of electro- 
magnetic or electrical NDE as a method of monitoring 
interstitial gas content. 

Sponsor : Office of Naval Research 

Faculty : R. M. Rose 

S taff : I. M. Puffer 

Undergraduate Student : J. Parse 


4. OTHER RELATED PROJECTS 


1 . 


Electromagnetically Driven Flows in Materials 

Processing 

J. Szekely 

See — Mathematical and Physical Modelling of 
Materials Processing, page 234 
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F. SOLIDIFICATION PROCESSING 


Faculty: 


T. W. Eagar 

Staff: P. Charreyron 

M. C. Flemings 

Z. Hu 

J. F. Elliott 

G . Kenney 

D. K. Roy lance 

M. Simpson 

D. R. Sadoway 

Y . Shiohara 

J. Szekely 

M. Yerebakan 


SUMMARY 


The Solidification Processing activities underway include 
the studies of Professor Flemings on behavior of semi-solid 
metals, purification and strengthening by fractional melting, 
continuous casting, control of ingot surface quality, and 
metal-matrix composites. Professor Elliott continues his 
work on formation of deoxidation products in steels. These 
studies are summarized below. 


In addition, there are ongoing studies by many faculty members 
on rapid solidification. There are studies on crystal growth 
being conducted by Professors Witt and Gatos. In addition, 
there are studies on experimental and mathematical modeling 
of solidification processes. These studies are summarized in 
separate parts of Section IV of this report. 


1. DEFORMATION BEHAVIOR OF SEMI-SOLID METALS 


This ongoing program is concerned with the fundamental rheo- 
logical behavior of semi-solid metals. Past work has been on 
non-dendritic "Rheocast" metals, current work is primarily 
on dendritic solids and on metal powders which are liquid 
phase sintered under pressure. Engineering applications are 
to better understand liquid phase sintering phenomena, and 
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to develop improved forming materials and purification meth- 
ods for materials which are partly liquid. 

Sponsor ; Army Research Office, Department of Energy 

Faculty ; M, C. Flemings 

Staff : P . Charreyron 

Graduate Student ; D. Pinsky 

Publications ; 

1. V. Laxmanan and M. C. Flemings, "Deformation of Semi- 
Solid Sn-15 Pet Pb Alloy," Met. Trans., Vol . 11A, No. 12, 
December 1980, pp. 1927 - 1937. 

2. T. Matsumiya and M. C. Flemings, "Modeling of Continuous 
Strip Production by Rheocasting , " Met. Trans., Vol. 12B, 
No. 1, March 1981, pp. 17-31, 

3. M. Suery and M.C. Flemings, "Deformation Behavior of 
Semi-Solid Alloys," accepted for publication, Met. Trans. 


2 ' STRENGTHENING OF METALS BY FRACTIONAL MELTING 

Ultra-high strength alloys are produced by partially melting 
alloy ingots and "squeezing" out interdendritic residual 
liquid, thereby drastically lowering impurity levels and 
residual second phases. Properties in wrought material pro- 
duced by these ingots are equivalent to, or better than, the 
best properties obtained in ingots of equivalent alloys pro- 
duced by other means, including consolidation of rapidly 
solidified powders. Work currently is on 7000 series alumi- 
num alloys. Mechanical properties obtained are in excess 
of 90,000 psi yield strength, 100,000 psi tensile strength, 
and 10% elongation. The process has been scaled up to 1.5 lb. 
ingots of aluminum, and plans are underway to extend the 
research to superalloys. 
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Sponsor : U.S. Army Armament Research and Development 

Command. 

Faculty : M. C. Flemings 

Staff : Z . Hu 

Graduate Students : M. Gungor, T. Piness 

Publication : 

1. F. E. Goodwin, P. Davami, M, C. Flemings, "Strengthening 
of Wrought Aluminum Alloys by Fractional Melting," 

Met, Trans. A., Vol. 11A, Nov. 1980, pp. 1777-87. 


3 . SURFACE QUALITY OF STEEL INGOTS 


This is an experimental modeling study of factors affecting 
ingot surface quality in steel ingots and continuous castings. 
Work has been on low melting point alloys in molds which have 
one transparent face and on steel solidifed around "dipped" 

bars. High speed photography has been used to record wave 

.¥ 

motion and meniscus behavior during filling, which is sub- 
sequently compared with surface and subsurface structures 
obtained. Wave motion and metal-mold heat transfer coeffic- 
ients have an important effect. The dampening effect of a 
steady D.C. magnetic field on convection alters the surface 
wave motion significantly, greatly improves surface quality, 
and enhances columnar grain growth. 

Sponsor : American Iron and Steel Institute 

Faculty : M. C. Flemings 

Graduate Students : M. Neff, D. Stemple 

Publication ; 

1. D. K. Stemple, E. N. Zulueta, and M. C. Flemings, 

"Effect of Wave Motion on Chill-Cast Surfaces," sub- 
mitted Met. Trans, B. 
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4, SECONDARY FORMATION OF DEOXIDATION PRODUCTS IN STEELS 


A study is underway of the influence of the principal para- 
meters of solidification, temperature gradient and rate of 
movement of tho solidification front on the composition and 
morphology of oxidos and oxysulfides that may form in 
steels after ladle deoxidation. The Fe-Si-O, Fe-Si-Mn-Q, 
and Fo-Si-Mn-O-S systems are being investigated. 

Sponsor : American Iron and Steel Institute 

Faculty : J. F. Elliott 

Graduate Student : D . -C . Hu 


5. METAL matrix composites 


Metal matrix composites, consisting of a light metal matrix 
and high strength ceramic fibers, show promise for engineer- 
ing applications where high modulus, high temperature 
strength, or high fatigue strength are required. Emphasis of 
this program is on fundamental aspects of setting and inter- 
face reactions in infiltration processes. 

Sponsor : U.S. Army Armament Command 

Faculty ; M. C. Flemings 
Graduate Student : D. P insky 


b . OTHER RELATED PROJECTS 

Solidification, in the generic sense, permeates much of 
the ongoing materials processing research. The following 
sections involve a variety of aspects of solidification 


processes : 


232 


h. Rapid Solidification Processing, page 159* 

B. Effect of Processing on Polymer/Composite 
Structure and Properties, page 185. 

C, Ceramics Processing Research, page 198. 

E. Welding Research, page 217. 

G. Mathematical and Physical Modelling of Materials 
Processing, page 233. 

H. Electroprocossing Research, page 241. 

I. Semiconductor Processing, page 245. 

J. Electronic Materials Research, page 260. 

K. Metals Processing, page 268. 
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G . MATHEMATICAL AND PHYSICAL MODELLING OF MATERIALS 
PROCESSING 

Staff: P* Brewor 

G, Chang 

M. Choudhary 

N. El-Kaddah 

I. Gabbalah 

J. McKelliget 
G. Oreper 

SUMMARY 

Professor Szekely's major area of research is mathematical 
and physical modelling of metals processing with major 
efforts underway on (1) turbulence. (2) elnctromagnetically 
driven flow, and (3) gas-solid reactions. Professor 
Roy lance continues his work on numerical modelling of poly- 
mer melt processes. Professor Sadoway is modelling electro- 
processing and Professor Flemings has modelling studies 
underway on electroslag casting and continuous casting. 


1. MATHEMATICAL AND PHYSICAL MODELLING OF METALS PROCESSING 
OPERATIO NS 

Turbulence Phenomena in Metals Processing 

There are many practical metals processing systems, where 
molten metal or slag phases undergo turbulent recirculating 
flow and where the characteristics of this flow play a major 
role in determining the overall process kinetics. The purpose 
of the research is to obtain an improved fundamental under- 
standing of the nature of these flows and then to apply this 
understanding to the solution of practical problems. 


Faculty: R. A. Brown 

T, W. Eagar 
M. C. Flemings 
R, M. Latanision 
D. K. Roylance 
D. R. Sadoway 
J. Szekely 
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This work involves mathematical modelling, viz the solution 
of the turbulent Navier- Stokes equations, physical model 
experiments, viz the use of laser anemometry to characterize 
the turbulence parameters and the flow fields in water model 
systems, and plate scale experiments to verify the models. 
The actual problems currently tackled include the erosion of 
blast furnace hearths, mixing and desulfurization kinetics 
in ladle metallurgical operations, AOD steelmaking, hot 
metal desulfurization, and dust removal from gases. 

Sponsors : National Science Foundation, primary sponsor; 

fellowship support from the Brazilian Government, 
the People's Republic of China, and Nippon Steel 
Corporation 
Faculty : J. Szekely 

Staff : G. Chang, N. El-Kaddah 

Graduate Students : L. Britto, R. Figueira, T. Kang, 

K. Shirable 


Electromagnetically Driven Flows in Materials Processing 

Electromagnetically driven flows play an important role in 
many metals processing operations, such as Electro-slag 
Refining, Electro-slag and Arc Welding, Electric Arc Steel- 
making, Induction Stirring, and the like. The purpose of 
this program is to develop a quantitative representation 
of the electromagnetically driven fluid flow phenomena in 
these systems. 

The current work includes the modelling of ESR systems 
concerning pool profiles and temperature profiles. The 
theoretical predictions have been found to be in excellent 
agreement with measurements obtained both in the laboratory 
and on industrial scale installations. 
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Work is also in progress on the modelling of Electro-slag 
and Arc Welding operations; the emphasis is to relate the 
role played by the principal operating parameters in determ- 
ining the structure and the properties of the welds pro- 
duced. Here again, the theoretical predictions were found 
to be in good agreement with the measurements. The research 
concerning the modelling of electric arc furnaces is aimed at 
representing the plasma region formed between the electrodes, 
with the ultimate objective of defining the optimum arc length 
for a given set of operating conditions. 

Electromagnetically driven flows also play an important role 
in the work aimed at modelling grain refining in space proces- 
sing applications. Here the specimens are positioned by 
electro-magnetic forces, which may also induce motion. The 
purpose of the work is to define the fluid motion resulting 
from the positioning forces, both during the flight and in the 
land based experiments. This understanding of the fluid flow 
field should be helpful in the development of grain refining 
modules . 

Sponsor : National Science Foundation 

Faculty : J. Szekely 

Staff : M. Choudhary, J. McKelliget, G. Oreper 

Graduate Student : A. Murthy 


Gas-Solid Reactions 


In terms of tonnages handled, the reduction of iron oxides and 
the gasification of coke in the iron blast furnace are perhaps 
the two most important metallurgical applications of gas- 
solid reaction systems. 


The work which is being carried out is largely experimental. 





aimed at developing an improved understanding, how the solid 
structure and the mode of preparation affects the rate at 
which coke particles react with C0 2 /C0 under conditions 

which are representative of the iron blast furnace. 

Parallel with this investigation the rate of iron oxide re- 
duction is also being studied, so that a composite picture 
may be developed of the coupled gas - solid reaction system, 
namely iron oxide reduction - coke gasification that controls 
the performance of the iron blast furnace. 

Sponsor : NATO - through Fellowship support 

Faculty : J. Szekely 

Staff : I. Gabbalah 


The Modelling of Rapid Solidification Processes 

This program of research involves the development of mathe- 
matical models of rapid solidification processes, with the 
objective of relating the processing conditions to the cool- 
ing rates that may be realized in given applications. 

Recent work includes the modelling of splat cooling, using the 
piston and anvil technique, and the operating of the twin roll 
quenching technique. 

In this latter work, allowance has been made for two dimension- 
al fluid flow, heat flow, and plastic deformation and the gov- 
erning equations were solved numerically . The theoretical 
predictions were compared with measurements and it was found 
that for this system there is only a very limited range of 
the process variables over which stable operation may be main- 
tained . 



Sponsor ; 

Faculty : 

Staff ; 

Student: 


Mexican Government (through Fellowship support) 
J. Szekely 
N. El-Kaddah 
M. Guttierez 


Electroslag Casting 

This program is a mathematical and experimental model of the 
electroslag casting process, including study of thermal be- 
havior, macrosegregation, and internal stresses. 

Sponsor ; National Science Foundation 
Faculty ; M. C, Flemings 
Stuff ; M. Simpson, M. Yerebakan 
Graduate Student: A. Jain 


2 • THE ROLE OF COPPER IONS AND OTHER CATHODIC DEPOLARIZERS , 

IN THE CORROSION OF ALUMINUM IN SEAWATER 

Four materials — titanium, plastics, stainless steel, and 
aluminum alloys — are considered as candidates for the con- 
struction of heat exchangers in ocean thermal energy conver- 
sion (OTEC) systems. Of these, aluminum alloys appear attrac- 
tive as the material of construction of heat exchanger tubes. 
Certain components (intake screens, etc.), however, must be 
made from materials which resist fouling such as alloys of 
copper. Copper ions, on the other hand, are likely to accele- 
rate the localized corrosion of aluminum alloys. In this pro- 
gram, the mechanism by which copper ions affect pit initiation 
on aluminum alloys in seawater is being pursued with a view 
toward corrosion control under service conditions. Likewise, 
the effect of bacterial fouling on this process will also be 
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studied. A scanning potential microprobe has been constructed 
for use in this work. This research to date has shown that 
the oxide films formed on aluminum in seawater are different 
from those formed in laboratory saline solutions, the most 
significant difference being the presence of substantial Mg 
concentration on the surface exposed to seawater. We have, 
during the past year, computer modeled the speciation of com- 
plex solutions and the adsorption of ions on aluminum surfaces 
in seawater. This program is part of the MIT/WHOI Joint Pro- 
gram in Materials Science and Ocean Engineering. 

Sponsor : MIT Sea Grant 

Faculty : R, M. Latanision 

Staff : P, Brewer (WHOX) 

Graduate Student: S. W. Smith 


3. COMPUTATIONAL AND EXPERIMENTAL STUDIES OF VISCOELASTIC 
FLOWS 


The ability to quantitatively predict the flow of viscoelas- 
tic fluids for describing the flow of polymeric solutions 
used in coating operations and in the design of dies for 
casting thermoplastics and thermosets is the long term goal 
of this project. Our research is aimed ultimately at devel- 
oping numerical methods for calculating the flow of visco- 
elastic fluids in such complicated geometries and at careful 
experimental measurements of velocities in these flows. 

These combined approaches will give us the ability to test 
constitutive equations for these flows. 

All present numerical methods fail when elastic effects are 
comparable to or greater than viscous effects. We have 
tracked this failure to inadequacies in the numerical schemes 
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for approximating the elastic portion of the stress, which 
may have steep gradients near corners and contractions in 
the flow. These steep gradients have been demonstrated for 
finite element based calculations with contravariant convect- 
ed Maxwell and second-order fluid models. Improvements to the 
finite element method that will allow accurate approximation 
of these stress fields are being developed. 

Sponsor : MIT Sloan Fund, union Carbide 

Faculty : r, A. Brown and R. C. Armstrong 

Graduate Students : A. Beris, M. Kim-E, P.-W. Yeh, M. A, 

Mendelson 

Thesis : 

1. M.A. Mendelson, "On the Numerical Solution of Viscoelastic 
Flow," M.S. thesis. 

Publications : 

1. M.A. Mendelson, P.-W. Yeh, R.A, Brown, and R.C. Armstrong, 
"Approximation Error in Finite Element Calculation of 
Viscoelastic Fluid Flows," J. Non-Newtonian Fluid Mech. , 
in press (1981) . 

2. H.M. Ettouney and R.A. Brown, "Effect of Heat Transfer on 

Melt/Solid Interface Shape and Solute Segregation in Edge- 
Defined Film-Fed Growth: Finite Element Analysis," J . 

Crystal Growth , submitted. 

3. C.J. Chang and R.A. Brown, "Effect of Steady Buoyancy- 
Driven Convection on Melt-Solid Interface Shape and Radial 
Solute Segregation in Vertical Bridgman Growth," J. Crystal 
Growth , submitted. 

4. G. Harriott and R.A. Brown, "The Fluid Mechanics of a 
Differentially Rotated Captive Drop," J. Fluid Mech. , 
submitted . 

5. Y . . Yamaguchi, C.J. Chang and R.A. "Jrown, "Axisymmetric 
Buoyancy-Driven Convection in a Vertical Cylinder Heated 
from Below," J. Fluid Mech. to be submitted. 


OTHER RELATED PROJECTS 


1. Numerical Modelling of Polymer Melt Processing 
D, K. Roylance 

See — Effect of Processing on Polymer/Composite 
Structures and Properties, page 193. 

2 . Modelling of Materials Supply, Demand, and Price 
J. P. Clark 

See — Materials Systems Analysis, page 211. 

3 . Heat and Fluid Flow (in Welding Processes) 

T. W. Eagar 

See — Welding Research, page 225. 

4 . Surface Quality of Steel Ingots 
M. C. Flemings 

See — Solidification Processing, page 230. 

5. Studies of Zinc Electrorefining 
D. R. Sadoway 

See — Electroprocessing Research, page 242. 

6 . Transport Phenomena in Improved Electrochemical 
Cell Designs for the Production of Magnesium 

D. R. Sadoway 

See — Electroprocessing Research, page 243. 

7 . Fluid Flow in Crystal Growth 
R. A. Brown 

See — Semiconductor Processing, page 258. 
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H. ELECTROPROCESSING RESEARCH 

Faculty: D. R. Sadoway Staff: J. Flint 


SUMMARY 


It is generally true that elec trodeposi ted metal is of 
better quality than the thermal reduction product. However, 
for many metals, notably the refractory metals, 
hydrometallurgical routes are closed because the 
decomposition potentials exceed that for hydrogen evolution 
and the metallic ions undergo redox reactions with the 
solvent. However, when these metals are deposited from high 
temperature non-aqueous electrolytes, the solid 
electrodeposites are powdery and/or dendritic. The 
challenge is to make coherent solid metal. 

Research in electroprocessing at MIT under the direction of 
Prof. D. R. Sadoway is concerned with the study of the 
fluid dynamics of the electroreduction process to determine 
how it may be modified to improve the quality of the 
deposit. A variety of experimental techniques, some 
uncommon to electroprocessing, is being used in this 
research. These include laser Schlieren photography, laser 
Doppler velocimetry, and frequency spectrum analysis. The 
work is sponsored both by government and industry. The 
specific projects involve fluid flow studies of zinc plating 
in aqueous and molten salt electrolytes (NASA) , cell design 
studies for magnesium chlorides electrolysis (DOE), digital 
signal analysis of manganese electrodeposition in molten 
chlorides (NSF) , and electroplating of molybdenum from low 
melting salts (Dow Chemical). This last project is 
particularly significant as it represents endorsement by 
industry of the electroprocessing effort and stands alone at 
MIT as a project involved in the primary processing of 
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••efractory metal. Looking to the future, one sees growth in 
the use of refractory metals as materials of construction in 
engineering. Their electrodeposition from molten salt 
electrolytes could weigh heavily in the extraction 
metallurgy of refractory metals. 

There is another emerging area of interest which will have a 
bearing on electroprocessing research in the Materials 
Processing Center. The problem of producing coherent metal 
deposits of high surface quality is one faced by electronic 
device manufacturers. According to industry sources, very 
little work is being done in the universities on this 
problem. For example, copper plating of multilayer circuit 
board interconnects from sulfate baths would fall into this 
category. It is expected that research activity will grow 
in this area almost exclusively through industrial support. 


1. STUDIES OF ZINC ELECTROREFINING 


To understand better the electrodeposition process, fluid 
flow patterns in the electrolyte are being followed by laser 
Schlieren photography in order to determine how mass 
transport affects the morphology of the metal deposit. Zinc 
is being studied as it is electroref ined both in aqueous 
electrolytes and in transparent molten salt electrolytes. 
Process variables such as current density and composition of 
the electrolyte are adjusted to change the morphology of the 
electrodeposit and, thus, to permit the study of the nature 
of electrolyte flow in relation to the quality of the 
electrodeposit . 

Sponsor ; National Aeronautics and Space Administration 
Faculty ; D. R. Sadoway 
Graduate Student; A. Abdelmassih 


Undergraduate Student; P. T. Rogers 
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Theses : 

1. A. Abdelmassih, "Fluid Flow Sutdieo of Zinc Electrolysis 
in Acid Chloride Solutions," due January, 1902. 

2. p. T. Rogers, "A Laser Schlieren Study of Electrolyte 
Flow Pattersn During the Eloctroref ining of Zinc from 
ZnC -LiClKCl Melts," S.B., June 1902. 

•I* M 

2 . TRANSPORT PHENOMENA IN IMPROVED ELECTROCHEMICAL CELL 
DESIGNS FOR THE PRODUCTION OF MAGNESIUM 

Increased use of magnesium as a structural material in 
transportation vehicles will significantly reduce fuel 
consumption. At present the restricted use of magnesium by 
the automotive industry is as much the result of inadequate 
capacity as price. To learn how transport phenomena in 
molten salt electrolysis cells affect magnesium production 
rates the effects of (a) forced convection of the 
electrolyte, (b) electrode configuration: vertical vs. 

horizontal; monopolar vs. bipolar, and (c) applied 
electromagnetic field will be studied. With the use of 
transparent laboratory magnesium electrolysis cells, not 
only will cell operating characteristics be measured, but 
the electrolysis process itself will be observed visually. 
The result of this work will serve to guide future 
developments in high yield cell designs. 

Sponsor : Department of Energy 

Faculty : D. R. Sadoway 

Graduate Student: S. Dokras 


3. HIGH PURITY MANGANESE BY FUSED CHLORIDE ELECTROLYSIS 


Electrochemical cells have among the lowest metal production 
rates of any metallurgical reactor. Molten salt 



electrolysis cells are particularly poor in this regard 
owing to their primitive designs Cor the most part. To 
determine whether digital signal analysis by fast Fourier 
transform of the cell voltage during operation can serve as 
the basis of diagnosis instrumentation for optimization and 
automation of electrolytic metal produetin processes the 
electrodeposition of manganese from molten alkali chlorides 
is being studied. 

Sponsor ; National Science Foundation 
Faculty : D. R. Sadoway 

Staff: J. Flint 


4. ELECTRODEFQSITION_ of molybdenum from low melting salts 

High purity molybdenum is a very desirable material due to 
its strength to weight ratio and corrosion resistance 
(particularly with respecc to molten metals). 
Electrodeposition of molybdenum is attractive because with 
this process come the prospects of plating a refractory 
metal coating onto less cost.ly base metal. 


Sponsor : Dow Chemical 

Faculty : D. R. Sadoway 
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I. SEMICONDUCTOR PROCESSING 



Faculty: R. 

C. 

Armstrong 

Staff: E. 

D, Bourret 

R. 

A. 

Brown 

T. 

Carlberg 

P. 

K. 

Iloupt 

J. 

R. carruthers 

T. 

B. 

King 

K. 

Crouch 

W, 

M. 

Rohsonow 

J. 

Fan 

A. 

F. 

Witt 

C. 

J . Herman 




Z. 

-J. Xing 




Q. 

-M. Zhou 


SUMMARY 

The primary thrust of the semiconductor processing group is 
directed (1) at advancing the theoretical basis for bulk 
growth of elemental and compound semiconductors in 
single crystal form and (2) at developing new experimental 
approaches by which semiconductor matrices with signifi- 
cantly improved crystalline and chemical perfection can be 
obtained. This effort, which evolved primarily with NASA 
and Air Force support, has recently been strengthened 
through DARPA which made it possible to acquire all instru- 
mentation required to engage in a meaningful and compre- 
hensive study of most advanaed approaches to silicon crystal 
growth. The research base, moreover, could be enlarged 
through an industrial grant (Mobil) , generous donations of 
instrumentation for both crystal growth and characterization 
by industry (T.I., BTL, Western Electric, and Semicon) and 
through industrial collaboration in sponsored research 
(H.P.), Implementation of the projected research expansion, 
directed toward the capability of growth of 4" diameter 
silicon crystals, has now been made possible through the 
allocation of additional laboratory space which is currently 
being modified to meet requirements. To broaden the theor- 
etical and experimental research potential of the group, 
both intra- and interdepartmental programs were established 
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which involve the areas of process metallurqy Professor 
King), heat transfer (professor Rohsenow) , mass transfer 
Professor Brown) and systems control (Professor Houpt) . 
Through NASA sponsorship, a cooperative research effort with 
the French Atomic Energy Commission (C.E.N. Grenoble) on 
solutal convection in melt growth systems could be expanded. 


1. ANALYSIS OF CRYSTAL GROWTH AND SEGREGATION IN VERTICAL 
BRIDGMAN CONFIGURATION 

Crystal growth and segregation in vertical seeded Bridgman 
configuration was experimentally studied using Ge-Ga as a 
test system. Using coded interface demarcation (CID) and 
spreading resistance measurements, it was possible to 
determine in post-growth analyses the shape and location of 
the solidification isotherm within the hot-zone as a 
function of time and to establish quantitatively the axial 
and radial segregation behavior. 

It is found that, depending on the axial thermal gradient in 
the growing crystal, the microscopic rate of growth does not 
reach a steady state value and deviates from the mechanical 
lowering rate by up to and in excess of ±50%. Upon arrest of 
ampoule lowering, thermal (and chemical) equilibration for 
periods ranging from 5 minutes to 120 minutes resulted in 
all instances in a relocation of the solidification isotherm 
(crystal-melt interface). 

Quantitative composition analyses based on spreading 
resistance measurements indicate that radial segregation 
increases with continuing growth, but is significantly less 
than predicted by the Corr iell-Sekerka theory. 

Anomalous, quantitatively reproducible segregation behavior 
is observed for gallium in germanium. Contrary to normal 
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segregation theories, it is found that for growth rates of 
less than 1 ym/s, the concentration of Ga in Ge increases 
with decreasing rates of growth. There is evidence that an 
inverse relationship for the interface distribution 
coefficient dependence on the growth rate exists. The 
presently reported results are used as a basis for the 
design of a Bridgman system in which two-dimensional heat 
flow is established about the crystal-melt interface. 

Sponsor : National Aeronautics & Space Administration 

Faculty : A. F . Witt 

Staff : J.R. Carruthers 

Graduate Student : C. Wang 

Publication : 

1. C. Wang, J.R. Carruthers and A.F. Witt, "Analysis of 
Crystal Growth and Segregation in Vertical Bridgman 
Configuration," Proc. of the ICVGE-5/ACCG-5 , Coronado, 
CA, July 19-24, 1981. 

2 . THE EFFECT OF DIRECT CURRENT ON CRYSTAL GROWTH FROM THE 

MELT: InSb 


Passage of direct current across a crystal-melt interface 
has been shown to affect the growth rate and segregation on 
a microscale. In this study, pulses of direct current with 
predetermined amplitude/time characteristics were passed 
across the liquid/solid interface during Czocbralsk.i pulling 
of doped InSb crystals for the purpose of producing spa- 
tially uniform alternating regions of high and low impurity 
concentration, [Such a "superlattice" is expected to func- 
tion aa a matrix for a tunable spin flip Raman laser (SERL) 
to operate in the far infrared region.] The influence of 
temperature control and crucible design on current-induced 
layer characteristics was quantitatively investigated 
through the use of simultaneous interface demarcation during 
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current pulsing. The experimental results were compared to 
a theoretical analysis utilizing a finite-difference numer- 
ical solution to the time-dependent energy transport 
equations and found to be in quantitative agreement. A 
theoretical transient segregation analysis was performed 
using experimental microscopic growth rate data in conjunc- 
tion with the numerical solution of the time-dependent mass 
transport equation calculated by J. Pa vie/.’,. The results 
are found to be in qualitative agreement Cor long duration 
(20s) pulses. It is shown that the application of direct 
current pulses during crystal growth from the melt provides 
a means for the in situ measurement of the Peltier coeffi- 
cient of a solid semiconductor in contact with its own melt. 

Sponsor ; National Aeronautics & Space Administration 
Faculty : A. F. Witt 

Graduate Student: M. Wargo 


3. THEORETICAL APPROACH TO THE DESIGN OF A VERTICAL 
BRIDGMAN GROWTH CONFIGURATION 

Heat transfer in a vertical Bridgman crystal growth system 
is analyzed with the intent of achieving optimized axial and 
radial temperature gradients in the vicinity of the growth 
interface through controllable experimental variables. The 
Bridgman system considered is composed of a hot and cold 
heat pipe separated by a region called the gradient zone in 
which the growth interface is located. The gradient zone 
comprises elements which prevent the establishment of 
satisfactory thermal gradients near the interface. A 
one-dimensional heat transfer model that neglects radial 
temperature gradients in the charge has been developed for 
the determination of the axial temperature distribution. 
Effects of charge motion, thermal coupling to the heat 
pipes, charge length, gradient zone length, latent heat, 
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conductivity change at the interface, crucible conductivity, 
and thickness could thus be demonstrated. Thermal coupling 
to the furnace (expressed by the non-dimensional Biot 
parameter, Bi), gradient zone length and conductivity change 
were shown to be parameters which strongly affect the axial 
gradient in the liquid at the interface. The origins of 
radial temperature gradients in the gradient zone has been 
considered in an extension of the one-dimensional model. It 
could be shown that Bi, as well as the temperature 
difference between the furnace and charge, in the gradient 
zone have a pronounced effect on radial temperature 
gradients. Strategies aimed at reducing critical radial 
gradients were developed and provide the basis for further 
analysis and systems design. 

Sponsor ; National Aeronautics & Space Administration 
Faculty : W. M. Rohsenow, A. F, Witt 

Graduate Student : T. Ja sin ski 

Publications : 

1. T. Jasins.i, W.M. Rohsenow, A. F . Witt, "Vertical Bridgman 

Type Crystal Growth: A Heat Transfer Analysis. Part I: 

Axial Temperature Distribution According to 
One-Dimensional Models," submitted to J, of Crystal 
Growth „ 

2. T. Jasinski, W.M. Rohsenow, A.F. Witt, "Vertical 

Bridgman Type Crystal Growth: A Heat Transfer Analysis, 

Part II: Radial Temperature Distribution According to 

Two-Dimensional Models," submitted to J. of Crystal 
Growth . 

4. THEORETICAL ANALYSIS OF DIRECTIONAL MELT-BACK 


A theoretical analysis of directional back-melting, which 
constitutes the initial step in all seeded crystal growth 
experiments, has been conducted for binary systems. It is 
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found that the achievenent of steady state conditions is 
strongly influenced by the magnitude of the liquidus-solidus 
separation; while steady state is readily reached in "doped" 
semiconductor systems, it cannot be reached in systems with 

large liquidus-solidus separation such as Hg A A Cd A „Te, for 
example. The analysis reveals, as a major complication in 
seeded vertical Bridgman growth, that upon initiation of 

growth the crystal-melt interface is not, as generally 
assumed, at the liquidus temperature, but at a temperature 
which is significantly lower than predicted from equilibrium 
phase diagrams. The equilibration process, following back- 
melting, has been investigated for the case where the seed- 
melt system prior to back-melting is not at equilibrium 
composition (CS = C L /k 0 ) . In this case, the behavior of a 
given system after arrest of back-melting is strongly 
affected by the respective values of the initial solid and 
liquid composition; the result will be spontaneous solidifi- 
cation, melting, or equilibration by solid and liquid state 
diffusion prior to controlled regrowth. The primary theore- 
tical conclusions of this analysis are currently subject of 
experimental verification. 

Sponsor ; National Aeronautics & Space Administration 
Faculty : A. F. Witt 

Staff : E. D. Bourret 

Publication : 

1. E. D. Bourret and A. F. Witt, "Theoretical Analysis of 
Directional Melting," Proc. of the ICVGE-5/ACCG-5 , 
Coronado, CA, July 19-24, 1981. 

5. DYNAMIC OXYGEN EQUILIBRIUM IN SILICON MELTS DURING 
CRYSTAL GROWTH BY THE CZOCHRALSKI TECHNIQUE 


A model for the dynamic oxygen level in silicon melts during 
Czochralski growth was established on the basis of balanced 
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oxygen fluxes through the melt. The theory predicts first 

order effects in good agreement with experimental results 

and accounts well for both absolute oxygen concentration and 

its axial variation in grown crystals. Basic approximations 

in the model limit its usefulness for a detailed analysis of 

second order effects. On the basis of this model it is 

proposed that axially uniform oxygen incorporation into 

growing crystals can be achieved by compensating the 

continuously increasing surface ratio (A s /A c ) through 

controlled immersion of auxiliary quartz surfaces or 

alternately through a controlled decrease in 5 , the 

c 

diffusion boundary layer thickness, at the crucible melt 
interface . 

Sponsor : Defense Advanced Research Projects Agency 

Faculty : T. B. King, A. F. Witt 

Staff : Te Carlberg 

Publication : 

1. T. Carlberg, T. B. King and A. F. Witt, "Dynamic Oxygen 
Equilibrium in Silicon Melts During Crystal Growth by 
the Czochralski Technique," J. of the Elec,' Soc,, 

129 ( 1 ) 189 ( 1982 ). 

6. DYNAMICS OF OXYGEN INCORPORATION DURING CZOCHRALSKI 
SILICON GROWTH 


The incorporation and distribution of oxygen during 
conventional Czochralski silicon growth has been studied in 
commercial crystals and was compared with a theoretical 
model recently developed. The comparison allows conclusions 
concerning the transport of oxygen from the crucible into 
the melt. It was found that the flux of dissolved oxygen 
into the melt is diffusion controlled with both melt flow 
and melt temperature affecting the dissolution rate. 
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The behavior of oxygen in silicon was examined by several 
experimental techniques? infrared spectroscopy was used to 
measure interstitial oxygen, spreading resistance to measure 
electrically active oxygen, and chemical etching to reveal 
crystallographic defects related to oxygen. The findings of 
all three techniques indicate non-uniform dissolved oxygen 
and precipitate distribution attributable to the crystal 
growth process and to automatic growth control functions. 

Sponsor ; Defense /advanced Research Projects Agency 
Faculty : A. F. Witt 

Graduate Student : D. Bliss 

Publication : 

1. D. F, Bliss, "Dynamics of Oxygen Incorporation During 
Czochralski Silicon Growth," S.M. Thesis, Department of 
Materials Science & Engineering, June, 1981. 


7 . AUTOMATIC DETERMINATION OF PIP PARAMETERS FOR TEMPERA- 
TURE CONTROL IN CRYSTAL GROWTH SYSTEMS 


The most common means of temperature control in crystal 
growing systems uses a thermocouple for temperature feedback 
and a standard proportional plus integral plus derivative 
(PID) series compensator. Owing to the length of the time 
constants in these systems, manual online tuning procedures 
tend to be tedious, time consuming, unrepeatable and 
unreliable. An effort was made to develop an off-line 
computer-aided procedure for choosing acceptable PID 
controller parameters. 

The process of choosing the controller parameters can be 
divided into the following six steps. (1) Quantify system 
performance: In this case, acceptable performance is 

defined in terms of the thermocouple response to step 
changes in commanded temperature. The requirements are zero 
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steady state error, no overshoot and settling time between 

one-tenth and one-half that of the uncontrolled plant, (2) 

Formulate a plant model: A simple first-order lag was used 

to relate thermocouple output to power amplifier input. 

Thus, (dx/dt) = x/T + (Ku)/T, where x = thermocouple output 

(volts), u = power amplifier input (volts), K = plant DC gain 

(volts/volt) and T = plant time constant (secs). (3) Design 

the controller: Since we are restricted to the PID control 

structure, the controller output (i.e. power amplifier 

t 

input) is given by u = -G^(x-x r ) - G 2 / (x-x r )dt - 
G 3 (d/dt) (x-x r ) where x r is the reference signal to be 
tracked and Gj_-G 3 are constant gains to be determined. The 
simplest approach is to use a PI compensator: choose the 

compensator zero to cancel the plant pole; choose the 
compensator gain such that the DC loop gain is equal to the 
desired closed loop pole. Therefore, G^ = 10/K, G 2 = 10/KT 
and G 3 = 0. (4) Formulate a discrete time plant model: The 

following difference equation can be written: x(nP+P) = 

Ax(nP) + Bu(nP), where A = exp(-P/T) , B = K(l-A), n = number 
of the sampling interval and P = duration of sampling 
interval (secs). (5) Evaluate the plant parameters: A 

computer with appropriate A/D conversion hardware is vsed to 
monitor input (u) and output (x) as arbitrary excitations 
are manually applied to the system. A least squares 
estimate for A and B is then found by simple linear 
regression. (6) Evaluate PID controller gains: Solving the 

d.i fference equation and substituting in G-^ and G 2 yields G-^ 

= ( 10/B ) ( 1-A ) and G 2 = (-G-j^ In A)/P. A digital simulation 
of the proposed closed loop system was implemented. A 
computer code was written to perform the linear regression 
parameter identification from stored data. 

Sponso r : National Aeronautics & Space Administration 

Facul ty : P. K. Houpt, A. F. Witt 

Graduate Student: G. Goodman 
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8 • APPLICATION of soft mold total liquid encapsulation to 

GROWTH OF CdTe BY THE VERTICAL BRIDGMAN TECHNIQUE 

An analysis of reported work on growth of CdTe suggests that 
the failure to achieve adequate crystalline perfection is 
largely related to deficient heat flow control (unfavorable 
and changing growth interface morphology) , heterogeneous 
nucleation at crucible walls, stresses induced by the 
confinement, and to inadequate control over stoichiometry. 

A growth system is currently being developed in which the 
crystal-melt interface morphology can be changed in a 
controlled manner from concave through planar to convex by 
making use of two coaxially aligned heat pipes separated by 
a "graded" heat guide zone. Preliminary studies indicate 
that complications arising from solid confinement materials 
can be avoided through the use of soft-mold, total liquid 
encapsulation with B 2 0 3 . This approach is based on the use 
of three coaxial cylinders consisting of BN (ampoule) , 

B 2 0 3 and porous high purity graphite surrounding the cast 
CdTe charge. The ampoule system is configured so as to 
provide for isolation of the liquid encapsulated crystal 
from the confining graphite and BN cylinder prior to 
solidification of the encapsulant. 

Sponsor : Defense Advanced Research Projects Agency 

Facul ty: A. F. Witt 


9. LEC GROWTH AND CHARACTERIZATION OF InP 


With the cooperation of the Electronic Materials Group at 
the Lincoln Laboratories, a comprehensive study of synthesis 
and LEC growth of InP was initiated at AFCRL . In initial 
growth experiments, four out of four undoped charges yielded 
multiple-twinned to polycrystalline crystals. After 
establishing rigorous cleaning procedures, revising the 
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out-baking process, replacing a conventional susceptor with 
a susceptor of UHP Poco graphite and using B 2 0 3 of higher 
purity, nine out of eleven charges yielded untwinned single 
crystals. The defect density achieved with the modified 
growth procedure in which the thermal configuration remained 
unaltered was in the low 10 3 range. Attempts to establish a 
meaningful relationship between the charge carrier 
mobilities of InP ingots and those of crystals grown by the 
LEC method were unsuccessful until quite recently when, with 
the modified growth procedure, it was found that a 
correlation does exist. Details of this correlation are 
currently under investigation. During this study it was 
found that Hall mobilities determined on the same samples in 
two laboratories yielded noticeably different values. A 
study of the causes for the differences has been initiated. 

The primary thrust of this investigation is aimed at the 
effects of solid solution hardening and the influence of 
liquid encapsulation on growth and segregation. A 
complementary study of the specific effects of liquid B 2 0 3 
encapsulation is carried out on campus. 

Sponsor ; Air Force - Rome Air Development Center 

Faculty ; A. F. Witt 

Staff ; J. Fan 

Graduate Student; B. Ahern 


10. EFFECTS OF LIQUID ENCAPSULATION ON CRYSTAL GROWTH AND 
SEGREGATION DURING CZOCHRALSKI PULLING 

The purpose of the present study was to determine the 
effects of liquid encapsulants on the segregation and growth 
behavior in Czochralski configurations. For this purpose, 
Ge-Ga was selected as a test system since it permits growth 
with and without liquid encapsulation and thus an 
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unambiguous analysis o f the effects oE encapsulants on the 
thermal melt characteristics and on the growth and 
segregation behavior. The melt characterization was 
conducted with an assembly o£ stationary and rotating 
thermocouples; the growth and segregation behavior was 
studied making use of coded growth interface demarcation and 
spreading measurements. 

Preliminary results obtained indicate the establishment of 
enhanced turbulent melt convection as a result of 
B 2 0 3 encapsulation. It is also observed that under constant 
thermal boundary conditions microsegregation in the central 
facet growth region changes from growth rate controlled to 
boundary layer controlled with use of B 2 0 3 as 1 "t u ^3 
encapsulant. A detailed study of the observed pronounced 
differences in growth and segregation behavior with and 
without liquid B 2 Q 3 is in progress. 

Sponsor : Air Force - Rome Air Development Center 

Faculty ; A. F. Witt 
Staff : Q.-M Zhou 

Publication ; 

1. Q.-M. Zhou and A, F. Witt, "Effects of Liquid Encapsula- 

tion on Crystal Growth and Segregation During 
Czochralski Pulling," Proc. of the ICVGE-5/ACCG-5 , 
Coronado, CA, July 1924, 1981. 


11. EFFECTS OF CONICAL INFRARED REFLECTORS ON CRYSTAL GROWTH 
AND SEGREGATION DURING CZOCHRALSKI PULLING 


In an effort to establish optimized crystal growth con- 
ditions for InP it was found that the heat transfer con- 
ditions across the crystal-melt interface and in the lower 
part of the crystal largely control the growth behavior, 
segregation, and crystal perfection. Controllability of 
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heat transfer in the critical regions is studied through the 
use of conical IR reflectors installed coaxially about 
growing crystals. To avoid complications associated with 
high pressure growth conditions, Ga-doped Ge is used as a 
model system. 

It is found that with the installation of the reflector, 
thermal asymmetry effects on growth become immeasurable and 
that the microscopic rate of growth becomes constant at all 
parts of the crystal-melt interface. Segregation studies 
indicate that the effective segregation coefficient during 
growth with a reflector increases significantly, which 
suggests a reduction in convective melt flow and a 
corresponding increase in the diffusion boundary layer 
thickness. Associated with the increase in the diffusion 
boundary layer thickness is an increased sensitivity of 
segregation to variations in bulk convective flow. Related 
to it, pronounced variations in dopant incorporation 
(segregation) are observed under conditions of constant 
microscopic rate of growth. These composition variations 
remain in all instances confined to the central facet growth 
region? peripheral off-facet regions exhibit uniform 
composition during all phases of growth. 

It is found that the use of IR reflectors decreases the 
sensitivity of the system to growth interface breakdown 
associated with constitutional supercooling. All crystals 
grown with IR reflectors were dislocation-free (EPD < 100/ 
cm ) while for growth without reflectors the dislocation 
density ranged from 300 to 2000/cirr. Detailed studies of 
the cause and effect relationships for growth with 
reflectors are in progress. 

Sponsor : Air Force - Rome Air Development Center 

Faculty : A. F. Witt 

Staff : Z.-J. Xing 
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1 2 • PREPARATION OF ORIENTED GaAs bicrystal layers by 
VAPOR-PHASE EPITAXY USING LATERAL OVERGROWTH 

A novel technique that utilizes vapor-phase epitaxy to grow 
bicrystal semiconductor layers with predetermined rotation 
axis, misorientation angle, and grain boundary plane was 
developed. The geometrical structure of the grain boundary 
in each layer is therefore completely specified. The 
technique has been demonstrated by using the 
AsC^ 3 ~GaAs-H 2 method to grow a series of GaAs bicrystals, 
each containing a [110] tilt boundary formed by a grain with 
a ( 111 ) B boundary plane and a grain rotated from (lll)B by a 
selected misorientation angle. The results of initial 
electrical measurements indicate that the height of the 
potential barrier associated with each grain boundary varies 
systematically with the misorientation angle. 

Sponsor : Air Force - Rome Air Development Center; Solar 

Energy Research Institute 
Faculty ; A, F, Witt 
Staff : J. Fan 

Graduate Student : J, Salerno 

Publication : 

1. J, P. Salerno, R. W. McClelland, P. Vohl, J.C.C. Fan, W. 
Macropoulos, C. 0. Bozler and A, F. Witt, "Preparation 
of Oriented GaAs Bicrystal Layers by Vapor-Phase Epitaxy 
Using Lateral Overgrowth," Proceedings of Annual 
Materials Research Society Meeting, Boston, MA, 

Nov. 16-19, 1981. 


13. FLUID FLOW IN CRYSTAL GROWTH 


This research program is a comprehensive theoretical and 
computational study ultimately directed towards a fundamen- 
tal understanding of che interactions of heat, mass, and 



259 


momentum transport in crystal growth from the melt. 

Emphasis has been placed on the floating zone, vertical 
Bridgman, Edge-Defined Film- Pod Growth processes, and on the 
microscale instabilities in directional solidification. 

In each research project, a combination of analytical 
modelling and computer-aided calculation is being used to 
develop qualitative understanding of the relevant physics 
and quantitative models of the selected crystal growth pro- 
cesses. To do this, state-of-the-art numerical methods have 
been developed for handling solidification problems. Five 
research projects are currently underway; each is described 
in detail in Chapter 3, Project B1 of this report. 

Sponsors : National Aeronautics and Space Administration, 

Mitsubishi Chemical Industries Ltd., Mobil 
Foundation, National Science Foundation 
Faculty ; R, A, Brown 

Graduate Students ; C. J, Chang, H. M. Ettouney, G. M. 

Harriott, L. H, Ungar, and Y, Yamaguchi 
Publications ; 

1. C. J. Chang and R. A. Brown, "Finite Element Calculation 
of Buoyancy-Driven Convection Near a Melt/Solid Phase 
Boundary," to appear in Proceedings of the Second 
National Symposium on Numerical Methods in Heat 
Transfer, Hemisphere Press (1981). 

2. L. H. Ungar and R. A. Brown, "The Dependence of the 
Shape and Stability of Captive Rotating Drops on 
Multiple Parameters," Philos. Trans. R. Soc. Lond. , in 
press ( 1981 ) . 

3. H. M. Ettouney and R. A. Brown, "Finite Element Methods 
for Steady Solidification Problems," J. Comp. Physics , 
submitted . 
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J. ELECTRONIC MATERIALS RESEARCH 


Faculty: H. C. Gatos 

J. Lagowski (Senior 
Research Associate) 


Staff: P. Becla 

S. Isozum.i 
M. Kaminska 
Y. Nanishi 
K. Wada 
X.-F. Yang 


SUMMARY 

The research program of the group is focused on semi- 
conductor materials and, aimed at the establishment of 
quantitative relationships underlying Crystal Growth 
Parameters - Materials Properties - Electronic Character - 
istics - Device Applications. The overall program evolves 
about the following main thrust areas: (1) Crystal Growth 

novel approaches to engineering of semiconductor materials; 
(2) Investigation of materials properties and electronic 
characteristics on a macro- and microscale; (3) Surface 
properties and surface interactions with the bulk and 
ambients; (4) Electronic properties controlling device 
applications and device performance. 


1. BRIDGMAN-TYPE APPARATUS FOR THE STUDY OF GROWTH-PROPERTY 
RELATIONSHIPS : GaAs 

In a high precision Bridgman-type apparatus the investi- 
gation of relationships between crystal growth parameters 
and the properties of GaAs crystals was continued. Key 
features of the system are a heat pipe for arsenic vapor 
pressure control and seeding without the presence of a 
viewing window. Pertinent growth parameters, such as 
arsenic source temperature, thermal gradients in the growing 

O V 
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crystal and .in the melt, and the macroscopic growth velocity 
can be independently controlled. During operation, thermal 
stability better than 0.02°C is realized; thermal gradients 
can be varied up to 30°C/cm in the crystal region and up to 
20°C/cm in the melt region? the macroscopic growth velocity 
can be varied from 50 pm/hr to 6.0 cm/hr. The density of 
dislocations was found to depend critically on As partial 
pressure and essen-tially dislocation-free crystals were 
repeatedly grown under As pressure precisely controlled at 
617°G. The free carrier concentration varied with As 
pressure variations. This variation in free carrier 
concentration was found to be associated with variations in 
the compensation ratio rather than with standard segregation 
phenomena. 

Sponsor : National Aeronautics and Space Administration 

Faculty : H. C. Gatos, J. Lagowski 

Staff: Y. Nanishi 

Graduate Student : J, Parsey 

Publication : 

1. J.M. Parsey, Y. Nanishi, J. Lagowsk.i and H.C. Gatos, 

"Br .idgman-Type Apparatus for the Study of Growth- 
Property Relationships: Arsenic Vapor Pressure-GaAs 

Property Relationship," accepted by J. Electrochem. Soc. 

2. OXYGEN- INDUCED LEVELS IN GaAs 

Utilizing our high precision Bridgman-type apparatus, the 
effects of oxygen (added in the form of Ga^O^ •'• n the welt) 
on the characteristics of bulk GaAs single crystals were 
investigated. It was found that oxygen induced or enhanced 
electron traps, particularly the level at 0.82 eV below the 
conduction band; it also decreased the carrier 
concentration, the minority carrier diffusion length and 
suppressed the 1.2 eV luminescence band. All these effects 
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were shown to ho ind ireetly induced by oxygon through the 
reduction of the Si concentration. SIMS analysis was 
extensively used in this study. 

Sponsors: National Aeronautics and Space Administration, 

National Science Foundation 
Faculty: ll. C. Gatos, O', Lagowski 

S tall ; M. Kaminska, K. Wad a 
Gr ndua tc S t ud e n t ; 0 . Pa r s ay 

P ubli cation : 

1. M. Kaminska, 0. Lagowski, J.M. Parsey, k. Wada, and 
U.C. Gatos, "Oxygen-Induced Levels in GaAs," presented 
at 1901 Symposium on GaAs and Related Compounds, Tokyo, 
Japan, 1901. 


3 , ORIGIN OF THE 0,82 eV ELECTRON TRAP IN GaAs AND ITS 
ANNIHILATION DY SHALLOW DONORS 


The concentration of the major electron trap (0,02 eV below 
the conduction band) in GaAs (Bridgman-grown) was found to 
increase with increasing As pressure during growth. It was 
further found that (for a given As pressure) the con- 
centration of this trap decreased with increasing con- 
centration of shallow donor dopants (Si, Se and To). Donor 

XV ~3 

concentrations above a threshold of about 10 cm led to 
the rapid elimination of the trap. On the basis of these 
findings the 0.82 eV trap was attributed to the antisite 
defect AsGa formed during the postgrowth cooling of the 

National Aeronautics and Space Administration, 
National Science Foundation 
H. C. Gatos, J. Lagowski 
M. Kaminska, K. Wada 
Graduate S tudent; J, Parsey 


crystals. 
Sponsor s : 


Staff: 
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J . I.aqownki, II. 0* Union, d.M, Par no y, R. Wad a, M. 
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In U.’.oohral nk l- type apparat un bulk UnAn qrowlh (ninule 
ovyrd aln with a thioknonn ot tho order v> l rail l imo torn) wan 
nueoenniully earrled out olootroepv taxi ally al 000°0 
( i no thermal liquid phano epitaxy initiated and muvtuined by 
pannlnu current thvouqh tho nubntrat o~uolut ion interlace) at 
qrowth eaten more than an or dor al maqnitude m out or than 
l* ho no obtained in t hoi mu l liquid phano epitaxy* Pool iminnry 
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proportion (oryntal lino par toot ion, carrier mobility, and 
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tho molt inn point ot UaAn . 

Sponsors National Aeronaut ion and Space Admlnintyat ion 
Faculty: It, C, da ton , d. Laqownki 

Uraduato Studentn: U, Boucher, 0, Wonq 


SELECTIVE EPITAXIAL UllOWTH OP UaAn BY LIQUID PHASE 
ELECmiKPITAXY 


Selective opitax ial qrowth ot UaAn wlndown (patterns) on 
SiO % ~mankod ( 100) oubot raton wan aohiovod by liquid phano 
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electroepitaxy (LPEE) . Prior to growth, wells (o£ the 
desired depth) were formed in the window areas by selective 
current-controlled dissolution. Overgrowth on the SiO^ 
layer was prevented by arresting epitaxial growth at the 
SiO^ layer height. Since both selective dissolution and 
growth rates are proportional to the current density, the 
depth of the wells and the thickness of the epitaxial layers 
can be precisely controlled by controlling the current 
density and process time. Another distinct advantage of the 
present approach is that, following the formation of the 
wells, electroepitaxial growth is initiated by simply 
reversing the polarity of the current passing through the 
substrate solution interface. Best results regarding the 
geometry of the wells and uniformity of the epitaxial layers 
were achieved^at relatively low current densities (less than 
about 15 A/cm ) . 

Sponsor : National Aeronautics and Space Administration 

Faculty ; H. C. Gatos, J. Lagowski 

Staff ; X.-F. Yang 

Graduate Stude nt; C. Wong 

Publications ; 

1. X.-F. Yang, L. Hwang and H.C. Gatos, ‘'Selective 
Epitaxial Growth of GaAs by Liquid Phase Epitaxy," 
accepted by J. Electrochem. Soc. 

2. J. Lagowski, H.C. Gatos, J.M. Parsey, K. Wad a, M. 
Kaminska and W. Walukiewicz, "Origin of the 0.82 eV 
Electron Trap in GaAs and Its Annihilation by Shallow 
Donors," accepted by Appl. Phys . Lett. 

6. "IN SITU" MONITORING THE LPE GROWTH 


The in situ continuous monitoring of the LPE growth rate and 
growth thickness of GaAs was achieved in an electroepitaxial 
growth configuration where the substrate back contact was 
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eliminated by positioning the substrate between two segments 
of GaAs solutions. Monitoring is based on the measurement 
of the total resistance of the substrate and the growing 
layer. Quantitative analysis indicated that this approach 
enables the determination of an LPE layer thickness as low 
as 0.1 pm. 

The in situ measurements of growth kinetics is uniquely 
suited for the study of growth-property relationships. 

Sponsors t National Aeronautics and Space Administration,. 

National Science Foundation 
Faculty : H. C. Gatos, J. Lagowski 

Staff : S, Isozumi 

Graduate Stud ent: A. Okamoto 


7. ELECTROEPITAXIAL GROWTH OF InP 

High quality InP layers were grown electroepitaxially from 
In-P solutions. The growth rate of these layers was found 
to be proportional to the current density. Growth rates 
greater by more than one order of magnitude than in thermal 
LPE were achieved. These results indicate that the growth 
of "bulk" InP crystals of epitaxial quality is possible. 

Sponsor : National Aeronautics and Space Administration 

Faculty : H. C. Gatos, J. Lagowski 

Staff : X.-F. Yang 


8. ENHANCEMENT OF INTERFACE STABILITY IN LIQUID PHASE 
ELECTROEPITAXY 


The thermodynamic as well as the dynamic constitutional 
supercooling criteria were derived for liquid phase 
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electroepitaxy (LPEB) taking into account the Peltier effect 
and electromigation, which control the growth interface 
temperature and the mass transport in the solution, 
respectively. It was shown that solute transport by 
electromigration enhances significantly the interface 
stability during growth. Thus, this treatment explains the 
experimentally attained stable electroepitaxial growth with 
velocities as high as 25 ym/min. Furthermore, it defines 
the geometry of the LPEE configuration and the growth 
parameters (current density and polarity, teirperature and 
substrate characteristics) which lead to the optimization of 
surface morphology. 

Sponsors ; DARPA, National Aeronautics and Space 
Administration 

Faculty ; H. C. Gatos, J. Lagowski 
Graduate Student ; A. Okamoto 
Publication ; 

1. A. Okamoto, J. Lagowski and H.C. Gatos, "Enhancement of 
Interface Stability in Liquid Phase Electroepitaxy," 
accepted by J. Appl . Phys. 


9. ISOTHERMAL GROWTH OF HgCdTe 


Hg Cd Te, n-type, layers with low carrier concentration 
* 1— x x 

and high carrier mobilities (approaching theoretical values) 
were obtained, for the first time, in the "as-grown state," 
i.e., without the standard post-growth annealing in Hg-rich 
atmosphere. This achievement was the result of the 
discovery that in isothermal growth, Hg Cd Te layers of a 

1"*X X 

given composition and optimum electronic properties can be 
grown only from a source of a well-defined composition and 
at well-defined temperature. Such layers also exhibit 
mirror-like surface morphology, and their composition is 
highly uniform radially and axially (to a depth up to about 
15 ym) . 
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Sponsor : DARPA 

Faculty: II. C. Gatos# J. Laqowski 

Staff: P, Do cl a 

t-sac kbsssc:.*':. q» 

Graduate Student: II. Ruda 

i ‘ If ~ ■ i‘ - :> r .» 

Publication : 

1. P. Becla, d, Laqowski, U.C. Gatos and II, Ruda, "A 
Modified Approach to Isothermal Growth of Ultra-High 
Quality HgCdTe for Infrared Applications," d. Electro- 
chem. Soc. . 128, 1171 (1981), 
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The metals processing effort is directed towards improvement 
of performance and usefulness of materials through modifi- 
cation and control of shape and internal structure. This 
includes the efforts of Professor King to analyze the 
interaction of reactive gases, introduced into a plasma arc, 
with iron alloys. Professors Flemings, Clark, Sadoway, 
Szekely, and Dr. Kenney have completed an assessment of the 
technology of magnesium production. The study of fast 
fluidized bed reactors is being continued by Professor 
Elliott. Aspects of separation processes are being investi- 
gated by Professor Ring. Also included in this section is 
Professor Rose’s program to develop high-field supercon- 
ducting composites for use in large magnetic fusion devices. 
Professors MacVicar and Bostock are concerned with the pro- 
cessing, structure, and property relationships of supercon- 
ducting materials and with the substitution of precious 
metals in standard electrical contact and connector applica- 
tions. Professor Pelloux is investigating the influence of 

* 
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processing procedures on the structure of zircaloy and 
nickel-base alloys with a goal of improving mechanical pro- 
perties and performance, Professor Kaplow's investigation 
of the potential of Motal-Insulator-Somiconductor junctions 
as photovoltaic devices is also reported, All of those 
research activities are presented in more detail below. 

Additional materials processing activities reported in other 
sections are cross-referenced at the end of this section. 


1 • KINETICS OF INTERACTION OF ARC PLASMAS WITH LIQUID 
METALS 


When a transferred plasma arc is used to melt metals , the 
gas used to support the arc is normally argon. Gases such 
as nitrogen are dissociated in the arc and the atomic spe- 
cies is, generally, much more reactive than the molecular 
species, it has been shown in this work that, when nitrogen 
is added to the plasma gas, it dissolves extremely rapidly 
in liquid rion and the amount of dissolved nitrogen reaches 
a steady-state solubility well above the value for 
equilibrium with molecular nitrogen. The steady state solu- 
bility increases as the percentage of the nitrogen in the 

gas increases, up to about 30% N , when a maximum value is 

2 

reached. The value is also increased markedly by increasing 
the oxygen content of the iron. A quantitative theory of 
the solution process has been established. Studies have 
also been carried out on iron-chromium, iron-nickel, and 
iron-manganese alloys. When hydrogen is also added to the 
plasma gas, the rate of nitrogen solution and the value of 
the steady-state solubility are markedly reduced. Hydrogen 
has also been shown to accelerate the rate of desorption of 
nitrogen. 

Work is continuing using hydrogen to deoxidize and decar- 
burize. 
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National Science Foundation 


Facul ty; T. B. King 
Graduate Student ; J. Katz 


2. PRIMARY P R ODUCTION OF MAGNESIUM 

m » « ■!■ | mitmmm mt wpw m — | — | , T i m il l i n 

Magnesium, with a density less than two thirds that of alu- 
minum, is an important domestically available light weight 
material which offers substantial potential energy savings 
to the transportation sector. The magnesium program was 
initiated in 1976 to consider the energy efficiency and eco- 
nomic viability of expanded magnesium production and use in 
the automotive sector, An international conference was 
also convened at MIT in 1977 to consider the production and 
use of magnesium, the most energy efficient automotive 
material with respect to net life cycle energy use. Work has 
recently been completed on a technical and economic 
assessment of current and proposed magnesium primary produc- 
tion technologies as a basis for future programming in pro- 
cess research and development. 

Sponsor: Department of Energy 

Facult y: J. P. Clark, M. C, Flemings, D. R, Sadoway, J. Szekely 

Staf f: G. B, Kenney 

Publication s : 

1. Report from the International Conference on Energy 
Conservation in the Production and Utilization of 
Magnesium, edited by Merton C« Flemings, efc al,, 

Conference held at MIT, Cambridge, MA. , May 1977. 

2. G. B. Kenney and J, P. Clark, "Magnesium: Energy 

Panacea?", American Metals Market - Minor Metals 
Section, August 18, 1977, 

3. G. B. Kenney and J. P. Clark, "An Analysis of Pricing in 
the Domestic Magnesium Industry," Proceedings of the 
Council of Economics of the JME, March 1978. 
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4. G. B. Kenney, An Analysis of the Energy Efficiency and 
Economic viability of Expanded Magnesium Utilizatio n, 
Garland Publishing Series of "Outstanding Dissertations 
Bearing on Energy," 1979. 

5. M. C. Flemings and G, B. Kenney, "Materials Research for 
the Fuel Efficient Automobile," Report to the 
Transportation Systems Center, U.S. Department of 
Transportation, October 1979. 

6. M. C. Flemings and G. 13. Kenney, "Materials Substitution 
and Development for the Light Weight, Energy Efficient 
Automobile," Report to the Office of Technology 
Assessment, Congress of the United States, February 
1980. 

7) Jl. B. Kenney, D. R. Sadoway and M. C. Flemings, "An 

Assessment of the Potential for Magnesium penetration of 
the U.S. Automotive Industry," Proceedings of the 37th 
International Magnesium Association. 

8. J. P. Clark and G. B. Kenney, "The Dynamics of 

International Competition in the Automotive Industry. 
Part I: A Framework for Analyzing the Dynamics of 

Intermaterial Competition," Materials and Society , 

Vol. 5, No. 4, pp. 383-309, (1981). 

9. G. B. Kenney and J. P. Clark, "The Dynamics of 

Intermaterial Competition in the Automotive Industry. 
Part II: A Case Study of the Demand for Magnesium" 

Materials and Society , Vol. 5, No. 4, pp. 391-406, 

(1981) . 

10. M. C. Flemings, et al , "An Assessment of Magnesium 
Primary Production Technology," Final Report to U.S. 
Department of Energy, Contract Number EX-76-A-01-2295 , 
February 1, 1981. 


3 • FLOW OF GASES AND SOLIDS IN A FAST FLUIDIZED BED REACTOR 

The fast fluidized bed reactor shows promise of being an 
important type of reactor in the non-fuels minerals 
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1. Kirkendall porosity, 

2. Residual stresses, 

3. Dewetting and Rayleigh coalescence of molten 
external tin layers, 

4. Diffusion kinetics in multif ilamentary 
materials. 

Sponsor : Department of Energy 

Faculty : R, M. Rose 

Staff : S. F. Cogan, I. M. Puffer 

Graduate Students : J. D. Klein, S-J. Kwon 

Undergraduate Students : L. Granick, P. Goldwhite 

Publications ; 

1. S.F. Cogan, D.S. Holmes, J.D. KLlein, and R.M. Rose, 
"Multifilamentary Nb^Sn by an Improved External 
Diffusion Method," Proceedings of ICMC Conference, 
Brookhaven national Laboratories, NY (May 1980). 

2. J. D. Klein, S.F. Cogan, G. Warshaw, N. Dudziak and R.M. 
Rose, "Properties of Microfilamentary Superconducting 
Composites Produced by a Modified External Diffusion 
Method," IEEE Trans, on Mag. , Vol. MAG-17, No. 1, pp. 
370-379 (January 1981). 

3. J. D. Klein, G. Warshaw, N. Dudziak and R.M. Rose, "On 
the Suppression of Kirkendall Porosity in 
Multifilamentary Superconducting Composites," IEEE 
Trans, on Mag, , Vol. MAG-17, No. 1, pp. 380-382 (January 
1981) . 

4. J. D. Klein, S.F. Cogan, S. Kwon and R.M. Rose, 

"Manufacture of Multifilamentary Nb Sn Superconductor by 
, 3 

the External Diffusion Method: Tin Coalescence and 

Diffusion," to be presented at the 9th Symposium on 
Engineering Problems of Fusion Research, Chicago 
Illinois, October 26-29, 1981. 
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6. SUPERCONDUCTING AND ELECTRICAL MATERIALS 

Rel ations hip of Materials Processing to Superconducting and 
Mechanical Properties 

»TC«-.m* a* 'rats; ’j»; rc.C :f- t-,— - -. i»«w «k i «i hh i ■» in n ja-w.-iri tnw «» 

In most of the high T^rH^ intermetallic compounds the 
micromechanisms by which the alloys become hard are such as 
to severely restrict dislocation substructure motion thereby 
limiting the possible plastic flow of the materials and 
cause these alloys to be very brittle. However, the speci- 
fic effects of internal defects and microstructure on the 
mode of fracture and degradation mechanisms of supercon- 
ducing properties are unknown. A microstructural charac- 
terization of the physical properties of CIS and A15 
superconductors is underway. Both Laves phase and com- 
parison A15 alloys are now under test. 

An analytical data base consisting of the mechanical proper- 
ties of as-grown single crystals with artificially induced 
defects, arc cast materials, and of various composite wire 
configurations is being developed specifically for ZrV^ 

( C15 ) and for Nb Sn (A15). Evaluations of surface fracture 

3 

and hardness numbers in pure and arc cast ZrV have been 

2 

already obtained. A four-point bending jig for determining 
the stress-strain relations in single crystal and selected 
polycrystalline specimens over the range of room temperature 
to 4.2K has been designed and constructed which permits 
elastic modulus and breaking stress determinations. Methods 
for producing very small specimens with geometrical 
integrity are being developed. 


Role of Structural Perfection and Texture in Refractory 
Superconductors 

Among the binary cubic Laves phase (C15) compounds, Zr V 2 r 
Hf V 2 , and their pseudobinary alloys, Hf x Zr-L_ x V2 have the 
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highest superconducting transition temperatures, in addi- 
tion, many of these materials undergo low temperature lat- 
tice transformations to lower lattice symmetry. 

Using a traveling heater solvent zone technique, we have 
grown, with groat difficulty, relatively large single 
crystals (2 mm - 7 mm) of SrV o , several extremely pure 
large-grained po.lycrystnlline rods, and a few 
polycrystalline rods containing approximately 1% by volume 
”r inclusions. Preliminary investigations of the resisti- 
vity, the low temperature structural transformation proper- 
ties (by neutron scattering and ultrasonics), and the 
critical field and superconduct 4 *g transition temperature of 
these samples showed that the properties of ZrV are extre- 
mely stress sensitive and that Cor all samples there appears 
to be an electronically driven transformation at about 1Q0K 
which is independent of the existence of the structural 
phase transformation. The lattice transformation of twinned 
single crystals to the rhombohedral structure at 96K is 
accompanied by the onset of an incommensurate charge density 
wave which becomes commensurate when the lattice transfor- 
mation is completed. Although there is considerable diffuse 
scattering in the "perfect" single crystal in the same tem- 
perature range, no charge density wave or structural trans- 
formation has been observed. Inductive (super-conducting) 
critical field measurements are underway on single crystal 
samples to determine if twinned and un twinned samples have 
different I!« (Q), If so, an attempt will be made to cor re- 
late the difference to Fermi characters. 

The absolutely hiqhost T values of any materials are 

G 

realized in the A15 class of alloys. The chief difficulties 
generally encountered in their fabrication all appear to 
relate to the development of a high degree of order and/or 
stoichiometry in the A15 phase. The resolution of questions 
concerning the exact nature of these materials is important 



both practically, in tounn ot undoi ntandirq the point lonntup 

v> t tnioronl ruotiuo to thoir both T nupoi conduct ivi ty » A 

d 

variety ot nppronohen in hoinu unod it thin laboratory to 

oharaotot n*o vnvioim AlMn ot int event tNh bn, Nb i«e, 

! t 

Nb AuW Thenc include ntiuUon ot (hr morphology ot urowth 
ot nputtered and ditiunlon i oaotod Al*> tilum and thou mu^ 
taro oxiden, invont iuat ionn ot epi t axial urowth median inmn 
and toxturlnvi In bulk lateriain, tahrioat ion ot ninule 
ovyntal and or ninule phuuc mnterialn, conn id o ra t tom* ot tho 
thermodynamic predict ionn ot order dinordoi ki notion and 
phane t told ntabi lily, and evatuationn ot oomponit ion, 


Wo havo found that tho otoatruwl inupeveonduot inq) r phyn p- 
oa X t and median i oa l pi oper t i on o t Nb bn d 1 1 1 mi i on l ayorn 
tabricatod on nimble ovyntal and polyovyntal lino Nb nub- 


nt rat on exhibit hchaviorn dor Ivina ivom nubntrato orieim 
tat ion and texturimt, reaped ively, Notably, when tho 

VlUu t i\n > .t i J . \ t.\f V iuk Al l. 1 i *i it .40 i 1 £ »vkvk V 1 »k 

[ivt \ i 1 v\t \u iviti uv isnl \ tu f ni ' iuyri m hiijiti 1 IvwHi |Y 

Alh-tUn, tho proport ion ot tho Nb bn tayor are undonirahly 
doaradod in all lenpeotn, Wo havo vooont ly tovmd that 
aitnonlimj treatment n ot niobium toil ntart inu dock drama t i > 
rally at toot tho phvuioal proport ten ot tho t Inal Nb bn 
layor and iln phane oontpon i t urn » A nyntomatio invont ivjat ion 
ot l ho ottootn ot Nb treatment*! amt nurtaoo condition** in 


underway. 


V'loot ion«boat\v Imauinvi Via 'I'hin Him v'ryoaonio tVvioen 


by inoorpoi at inn a ot yoqonio«dewar ntavic into a neanninq 
olootron mi ovouoopo , a two*Himomuonal imaqe ot tho non* 
equi l i hr him tvnponnon ot a mvpereonduetor to electron beam 
ir radial ion ran bo observed, irradiation ot a current** 
carry ! 11 q "weak 1 ink” t Urn with an olootron bourn cannon a 
ohanqe in tho renintive voltaqe ot tho link? thin ohamiw in 
voltaqe *a« a function of tho npat val pen it ion of tho boam 



278 


spot# can bo used to modulate the brightness of the 
microscope's CRT tube in much the same way as the emission 
of secondary electrons is used to produce conventional SEM 
pictures. However , the non-equilibrium response of super- 
conductors is not confined to resistive voltage changes but 
is a complex mixture of interchanges between Cooper pairs 
and quasi-particle populations (accompanied by emission or 
absorption of phonons). Past research using a variety of 
other detecting techniques has not sharply delineated the 
microscopic mechanisms underlying the non-equilibrium 
responses. This is especially true near the critical tem- 
perature for decay of the superconducting state. Using this 
new two-dimensional imaging technique to examine a variety 
of sample configurations (including "weak links" and super- 
conducting tunnel junctions), we expect to obtain new infor- 
mation concerning relaxation mechanisms in superconductors. 

A detailed investigation of the internal structures (or, the 
physical disposition) of metal films based on this mechanism 
also appears possible. Design and construction of a second- 
generation low temperature stage for insertion into an SEM 
is well underway. Electron irradiation effects on tunnel 
junction characteristics are the first experimental step. 
Planning and design of the experiments to be undertaken with 
the SEM stage are almost complete. 


Deformation Effects in Electronic Behaviors of Cryogenic 
Diodes 


In utilizing materials for practical devices considerable 
deformation results from ordinary processing techniques. 

The superconductivity of such materials also undergoes 
(sometimes dramatic) changes. Since Nb is the highest 
T^ element and the host matrix for the highest known 
T^ materials, it is particularly important to establish the 
connection between its microstructure and its supercon- 
ducting parameters. 
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Recently an acid-etch oxidizing technique has been developed 
Cor fabricating superconducting tunnel junction barriers on 
severely plastically-deformed (up to a true strain of 3.9) 
bulk Nb. The electron-coupled phonon spectra determined 
from the superconducting tunneling data shows two clear 
trends: those having a character similar to undeformed bulk 

Nb but predicting an increase of T^ over the bulk value; and 
those with distinctly different characters, predicting the 
decreased T^ of comparably deformed material. The effects 
of surface etching inherent in the barrier formation tech- 
nique with regards to those results is still to be deter- 
mined. The electron-phonon coupling parameter and effective 
Coulomb pseudopotential derived from the tunneling data 
indicate that a threshold of mechanical deformation exists: 
on one side of it, \ and a * do not agree with conventional 
theory, while on the other side of it, conventional theory 
is suitable for the results obtained. These data, and 
results on amorphous Nb films investigated by Stanford 
University and Bell Labs, suggest the possibility of a 
non-BCS mechanism. 

Since much of the tunneling we observe for the Nb-based A15 
substrates is very reminiscent of data from tunneling into 
highly deformed Nb surfaces and since very low effective 
tunneling T^'s indicate a damaged surface layer on the A15 
films, we have begun exploring specific configurations which 
can be reliably used to extract the energy gap information 
about superconductors. These include proximity effect junc- 
tions and point contact microbridge configurations. 

Precious Metals Elimination Initiative 


An investigation is underway to determine the feasibility 
and desirability of mounting a collaborative MIT-Industrial 
consortium effort to pursue mutually desirable avenues of 
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research with the objective of reducing or eliminating the 
use of precious metals in standard electrical contact and 
connector applications. Activities which duplicate lines of 
inquiry currently being pursued within industrial labora- 
tories or which have been found unproductive in the past 
would have low priority in comparison to activities relating 
to more exotic, longer range, innovative approaches, e.g,, 
intercalated graphite, integrated composites, surface 
treatments, preferred orientation processing, etc. The 
immediate task is to evaluate the degree of industry 
interest in such an effort and to identify industry can- 
didates ready to develop a mutually desirable research 
agenda and format. 

Sponsors : Department of Energy, West German Science 

Foundation, National Science Foundation, 

Honeywell Co. , Inc. 

Faculty ? J. Bostock, M.L.A. MacVicar 

Staff : A. Callegari, R. Gelinas, P. Stohr, P. Tedrow 

Graduate Students : S. Chen, Y. Epps, L. Salamanca-Riba 

Undergraduate Students ; K. Allen, S. Arney, J. Bellingham, 
Z. Bern, A. Cohen, J. Izatt, R. Jaimes, D. Kieda, 

E, Leiser, K. Sidikman 
Publications : 

1. M.L.A. MacVicar and J. Bostock, "Influences on 

2 

Electron-Coupled Phonon Spectra a F Obtained From 
Tunneling Data Deconvolutions," Physics of Transition 
Metals II-Univ, of Leeds, England, 8/80, p. 527. 

2. M.L.A. MacVicar, K.R. Milkove, W.N. Cheung and J. 
Bostock, "Tunelling Characteristics as a Diagnostic of 
Superconducting Electrode Disposition," 158th 
Electrochemical Society Meeting, Hollywood, Florida, 
10/80. Invited paper. 

3. M.L.A. MacVicar, J.F. DeBroux , V. Diadiuk and J. 

Bostock, "Effects of Nb Substrate Orientation on Nb Sn 

3 

Diffusion Layer Properties," 158th Electrochemical 
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Society Meeting# Hollywood, Florida, 10/80, Invited 
paper. 

4, M.L.A. MaeVicar and J. Boahock, "a P from Tunelling; 
Considerations and Influences," To be published, 

5, M.L.A. MaeVicar and P.C. Wyatt, "Carburetor," The World 
nook Encyclopedia , 10 02 Edition; Invited author. 

6, M.L.A. MaeVicar and P.C. Wyatt, "Drakes," The World Hook 
Encyclopedia , 1982 Edition; Invited author. 

7, M.L.A, MaeVicar, Encyclopedia of Materials Science and 
Engineering , subject editor and overview author for 
superconductivity, Pergamon Press Ltd. Encyclopedia to 
be published (19B2). 

8, P.L. Stohr, "High Resolution Imaging of Magnetic 
Structures and Inhomogeneities of a Superconductor Using 
Scanning Electron Microscopy," 16th International 
Conference on Low Temperature Physics, Los Angeles, 
California, 8/81; proceedings published in Physica , 

10713 , 1901, pp. 441. 


7. PROCESSING/PROPERTY RELATIONSHIPS IN ZIRCALOY AND NICKEL 
-DAS E ALLOYS 

Environmentally Induced Cracking Under Cyclic Loading of 
Nickel- Base Alloys Used., in Li g h t Water Reactors t . - 
Mo d i f i ca t Ion o f Ag ing 

This program is oriented toward study of the intergranular 
fracture of Inconel 600 under cyclic loading and electroche- 
mical control in PWR environments. During the past year 
particular emphasis has been placed on the synergistic 
effects of thermal treatment and cathodic polarization on 
the fatigue crack growth behavior of Inconel 600 at ambient 
temperatures and pressures. The essence of this work is 
that aging treatments of longer duration than those typi- 
cally used in the industry at present lead to both desen- 
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oitization (i»e., an increase in the Cr concentration in the 
grain boundary proximity following an initial sensitizing 
thermal treatment) and phosphorus desorption from the grain 
boundaries, the consequence of which is reduced suscep- 
tibility to intergranular cracking at cathodic potentials. 

Sponsor ; Electric Power Research Institute 
Faculty ; R.M. hatanision, R.M.N. Pelloux 
Staff; R.G. Dellinger 
Graduate Student; N. Moshier 


Crystallographic Tonturo Control in Zircaloy Tubes 

Mechanical properties, stress corrosion cracking resistance, 
and irradiation creep resistance of zircaloy cladding tubes 
can all be improved by taking advantage of crystallographic 
texture effects. The fundamental objective of this program 
is to conduct a systematic investigation of the manufac- 
turing parameters which control the formation of 
crystallographic textures in Zircaloy-2 and -4 tubing. 
Precise control of texture development during tube produc- 
tion must be based on; (1) close observation and detailed 
analysis of the process, and (2) a thorough understanding of 
the deformation behavior of the material. In order to 
accomplish these tasks, samples of Zircaloy from inter- 
mediate stages in the tube rocking process are being ana- 
lyzed dimensionally, micros true turally, and 
crystallographically . In addition, mechanical tests simu- 
lating the deformations incurred during tube production, arc 
being performed with similar analyses. The texture shar- 
pening and rotation during the intermediate recrystalliza- 
tion steps is also being evaluated. A computer model based 
on the micromechanics of deformation in Zircaloy is used to 
simulate the tube rocking process and to predict the 
resulting crystallographic textures. The main goal of this 
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research program Is to achieve process control in the for- 
mation of crystallographic textures. 

Sponsor! Exxon Nuclear Company, Inc. 

Faculty! E. M. Pelloux 
Graduato Student! J. Shewbridge 

mauMsam: •«&*=- - 


8 . PHOTOVOLTAIC DEVICES 
Summary 

Metal-Insulator-Semiconductor (MIS) junctions show promise 
for low cost production of efficient photovoltaic cello and 
other devices. An MIS junction consists of a insulator- 
metal sandwich in which the energy barrier height depends on 
the relative work functions of the metal and the semiconduc- 
tor and is strongly affected by the properties of the inter- 
facial layer. 

Metal-Insulator-Semiconductor Solar Cells (MISSC) are being 
fabricated using p-type silicon (boron doped, 0.5-ICh-cm 
resistivity, (100) and (111) orientations) and a low work 
function metal (Al, Co, Hf , Mn, Y). The wafer is first che- 
mically cleaned and 0.5 m thick Al film is evaporated on 

the back surface. This back contact is then heat-treated in 

— 5 o 

a vacuum better than 2 x 10 torr to 400-500 C to form an 
alloyed ohmic contact. When the equilibrium temperature has 
been reached, ultrapure dry oxygen is introduced in the 
chamber to grow a thin oxide on the front surface (10-15 A 
thick). This oxidation is performed after obtaining well- 
defined starting surface conditions, e.g. a chemical 
cleaning procedure has been developed to give reproducibly a 
6 A thick ultrathin native-silicon oxide. Our first diodes 
were fabricated using aluminum as a contacting metal. 

Device fabrication is completed by evaporating a copper film 
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(. «* 'r h'? 

on both sides to lower contact resistance and to attach gold 
wires with silver epoxy. A C94 mv barrier height for a 
completed device has been determined from X-V and C-V 
electrical measurements, 


MXSSC behavior and performance will be optimized in terms of 
silicon resistivity and orientation , silicon oxide thickness 
and growth conditions, and metal work function. Surface 
studies performed on silicon before and after oxidation 
(Auger, XPS) will be correlated with device characteristics 
and performance. 


Faculty ; R. Kaplow* 

Graduate Student ; M, Baddi 
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*We regret the death of a highly valued member of our 
faculty, Professor Roy Kaplow. 
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L. TECHNOLOGY TRANSFER AND INTERNATIONAL DEVELOPMENT: 


MATERIALS 

AND MANUFACTURING 

TECHNOLOGY 



Faculty: R. 

F. Baddour 

Staff: 

C. T. 

Hill 

J. 

H. Hollomon 


K. N. 

Rao 

J. 

Oliveira 


J. A. 

Hansen 


SUMMARY 

The Center for Policy Alternatives (CPA) under the 
directorship of J. Herbert Hollomon has conducted, over the 
past eight years, a series of policy-oriented studies on 
technological development in several relatively advanced 
developing countries. These have been directed by Dr. 
Nagaraja Rao, Senior Research Associate at CPA. Priority 
sectors defined in terms of technological sophistication, 
capital intensity, value added, and export potential have 
been studied in the cases of Brazil, Venezuela, Israel, and 
Korea. These projects have been supported directly by. the 
countries themselves or through loans and grants made to 
them by international agencies such as the World Bank and 
the U.S. Agency for International Development. Research 
teams include graduate students from M.I.T. Schools of 
Engineering and Management, and senior researchers from the 
countries sponsoring the study.- 

Although the objectives of these studies are generally to 
develop technological policy alternatives for the sponsoring 
country, much emphasis is placed on understanding the 
dynamics of the sectors through structured interviews with a 
large sample of firms in the leading manufacturing and 
materials processing sectors. 
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Technological Status and Dcvelrntnont of Porti 

fcass a s^^ ssajfcjBagss^gsjq^ lft t ^^ - ?^^ 


This Is a study at the Center for Policy Alternate, voo by the 
National Laboratory for Engineering and Industrial 
Technology { LNETI ) of Portugal and partially financed by the 
Government of Portugal through a project loan from the World 
Bank. 


The objective of the study is to assess the present status 
of teennoiogy in the .major industrial sectors of Portugal 
and examine alternative approaches to the upgrading of 
technology in industry and to strengthen the relationship of 
industry and public research and development institutes , and 
promote better utilization ol at andardi station f engineering * 
and quality control services. The Portuguese Government is 
actively promoting industrial development through a variety 
of incentive programs to enable industry to become 
competitive with those of the members ol the European 
Economic Community* which Portugal hopes to join in the 
mid-eighties * 

Plans for the improvement of product and process technology 
in a sample of one hundred Portuguese firms are assessed 
through structured interviews with entrepreneurs. 

Literature studies are conducted and expert opinions are 
sought on the potentials for diversification of Portuguese 
manufacturing and industries based on locally available 
mineral- agricultural* and ocean resources. Those economic 
and technical assessments will form the basis for LNETI to 
develop a mid-term technology development plan and to 
identify projects for funding by the Government or 
international financial agencies. 


Sponsor : 


Laboratorio Pa clonal de Kugonharia e Tee no log in 
Industrial - LNETI 
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Faculty ; R. Baddour, J. Ilollomon, Jao Oliveira 
Staff; K. Rao, T. Hill, J . Hansen 
Graduate Student : S. Bengali 

Research Collaborators from LNET .lt Campos Rodrigues, 
Director of the Research Project and his 
Associates 

Publications : 

1. K. Negara j a Rao, R. F. Baddour and C. T. Hill, The 

Ch a n g i ng World Chemical Indus t r y and the N e ed for New 
Strategies , "Strategic Aspects of Chemical Industry 
Development iri the Rapidly Industrializing Nations," 2nd 
World Congress of Chemical Engineering, IX Interamorican 
Congress of Chemical Engineering, Montreal, October 
1981. 

2. C. Rodrigues, "Report from Abroad," Policy Choices , 
CPA/MIT, Fall, 1981. 




